1 Plot the covariances

pops<-c("PWS","TB","SS")
covs<-data.frame()
Variance<-data.frame()
winsize<-c("50k","75k","100k","250k")

for (w in 1: length(winsize)){
    covs<-data.frame()
    for (p in 1: length(pops)){
        #covariance output file
        cov<-read.csv(paste0("~/Projects/Pacherring_Vincent/MD7000/3Pops_maf05_temp_cov_matrix_",pops[p],"_",winsize[w],".csv"))
        cov<-cov[,-1]
        
        ci<-read.csv(paste0("~/Projects/Pacherring_Vincent/MD7000/3Pops_maf05_",pops[p],"_Cov_CIs_bootstrap5000_",winsize[w],"window.csv"))
        ci<-ci[,-1]
        
        #reshape the matrix
        mat1<-cov[1:3,]
        mat2<-cov[4:6,]
        
        covdf<-data.frame()
        k=1
        for (i in 1:nrow(mat1)){
            for (j in 1:ncol(mat1)){
                covdf[k,1]<-mat2[i,j]
                covdf[k,2]<-mat1[i,j]
                k=k+1
            }
        }
        colnames(covdf)<-c("label","value")
        covdf$value<-as.numeric(covdf$value)
        covar<-covdf[grep("cov",covdf$label),]
        vars<-covdf[grep("var",covdf$label),]
        
        #remove the redundant values
        if (pops[p]!="SS") covar<-covar[!duplicated(covar[, 2]),] 
        if (pops[p]=="SS") covar<-covar[c(1,2,4),]
        
        #assign the starting time period and covering period values
        covar$year<-c(1,2,2)
        covar$series<-c("1991","1991","1996")
        
        vars$year<-c(1,2,2)
        vars$series<-c("1991","1991","1996")
        
        #assign population name
        covar$location<-pops[p]
        vars$location<-pops[p]
        
        #attach ci info
        covar$ci_l<-c(ci[1,2], ci[1,3],ci[2,3])
        covar$ci_u<-c(ci[4,2], ci[4,3],ci[5,3])
        
        #combine in to one matrix
        covs<-rbind(covs, covar)
        Variance<-rbind(Variance, vars)
    }
    
    covs$ci_l<-as.numeric(covs$ci_l)
    covs$ci_u<-as.numeric(covs$ci_u)
    ggplot(data=covs, aes(x=year, y=value, color=location, shape=series, group=interaction(location, series)))+
        geom_point(size=3, position=position_dodge(width = 0.1,preserve ="total"))+
        #geom_errorbar(data=covs, aes(x=year, y=value, ymin=ci_l, ymax=ci_u), width=.2, size=.2, position=position_dodge(width = 0.1,preserve ="total"))+
        geom_line(data=covs, aes(x=year, y=value,color=location, group=interaction(location, series)), position=position_dodge(width = 0.1,preserve ="total"))+
        ylab("Covariance")+xlab('')+theme_classic()+
        theme(axis.text.x = element_blank(),legend.title = element_blank())+
        geom_hline(yintercept = 0,color="gray70", size=0.3)+
        geom_errorbar(aes(ymin=ci_l, ymax=ci_u), width=.2, size=.2, position=position_dodge(width = 0.1,preserve ="total"))+
        scale_shape_manual(values=c(16,17),labels=c("1991-","1996-"))+
        scale_x_continuous(breaks = c(1,2))
    ggsave(paste0("../Output/COV/3Pops_Cov_overtime_CI_",winsize[w],".window.png"),width = 4.7, height = 3, dpi=300)
}    

2 Find regions with high covariances in each population

  • From Temporal Covariance analysis -output covariances for each time period

2.1 Plot the covariances across the genome

#Find the regions with a high temporal covariance 
pops<-c("PWS","TB","SS")
winsize<-"100k"
evens<-paste0("chr",seq(2,26, by=2))
cov.list<-list()
covs_all<-list()
k=1
for (p in 1: length(pops)){
    pop<-pops[p]
    iv<-read.csv(paste0("~/Projects/Pacherring_Vincent/MD7000/3pops_intervals_",winsize,"window.csv"), row.names = 1)
    if (p==3) {
        cov23<-read.csv(paste0("~/Projects/Pacherring_Vincent/MD7000/",pop,"_cov23_2017-2006_2006-1996_3Pops_",winsize,"window.csv"), header = F)
        covs<-cbind(iv, cov23)
        colnames(covs)[4]<-c("cov23")
        covs$index=1:nrow(covs)
        covs$color<-"col1"
        covs$color[covs$chrom %in% evens]<-"col2"

        covs[sapply(covs, is.infinite)] <- NA
        covs[sapply(covs, is.nan)] <- NA
        
        cov.list[[k]]<-covs
        names(cov.list)[k]<-paste0(pop,"_",winsize)    
        k=k+1
            
        y<-min(covs$cov23, na.rm=T)
        ymin<-ifelse (y<=-0.1,-0.1, y) 
        ymax<-max(covs$cov23, na.rm=T)
        ggplot(covs, aes(x=index, y=cov23, color=color))+
            geom_point(size=1, alpha=0.5)+
            theme_classic()+
            ylim(ymin,ymax)+
            scale_color_manual(values=c("gray70","steelblue"), guide="none")+
            ylab("Covariance")+xlab('Chromosome')+
            theme(axis.text.x = element_blank())+
            ggtitle(paste0(pop," ", winsize," window"))
        #ggsave(paste0("../Output/COV/3Pops.",pop,"_tempCovs_acrossGenome_",winsize[i], "Window.png"), width = 8, height = 2.7, dpi=300) 
        }
    else {
        cov12<-read.csv(paste0("~/Projects/Pacherring_Vincent/MD7000/",pop,"_cov12_1996-1991_2006-1996_3Pops_",winsize,"window.csv"), header = F)
        cov23<-read.csv(paste0("~/Projects/Pacherring_Vincent/MD7000/",pop,"_cov23_2017-2006_2006-1996_3Pops_",winsize,"window.csv"), header = F)
        cov13<-read.csv(paste0("~/Projects/Pacherring_Vincent/MD7000/",pop,"_cov13_2017-2006_1996-1991_3Pops_",winsize,"window.csv"), header = F)
        covs<-cbind(iv, cov12, cov23,cov13)
        colnames(covs)[4:6]<-c("cov12","cov23","cov13")
        covs$index=1:nrow(covs)
    
        covs$color<-"col1"
        covs$color[covs$chrom %in% evens]<-"col2"
    
        covs[sapply(covs, is.infinite)] <- NA
        covs[sapply(covs, is.nan)] <- NA
        
        cov.list[[k]]<-covs
        names(cov.list)[k]<-paste0(pop,"_",winsize)    
        k=k+1
        covsm<-melt(covs[,c("index","color","cov12","cov23","cov13")], id.vars = c("index", "color"))
        ymax<-max(covsm$value, na.rm=T)
        y<-min(covsm$value, na.rm=T)
        ymin<-ifelse (y<=-0.1,-0.1, y) 
        ggplot(covsm, aes(x=index, y=value, color=color))+
            facet_wrap(~variable, nrow=3)+
            geom_point(size=1, alpha=0.5)+
            theme_classic()+
            ylim(ymin,ymax)+
            scale_color_manual(values=c("gray70","steelblue"), guide="none")+
            ylab("Covariance")+xlab('Chromosome')+
            theme(axis.text.x = element_blank())+
            ggtitle(paste0(pop," ", winsize," window"))
        #ggsave(paste0("../Output/COV/3Pops.",pop,"_tempCovs_acrossGenome_",winsize, "Window.png"), width = 8, height = 8, dpi=300)    
    }
}

2.2 Find the outlier regions for each time period

#find how outliers overlap between different windows
cov12<-data.frame()
cov23<-data.frame()
cov13<-data.frame()

for (i in 1:length(cov.list)){
 if (grepl("PWS",names(cov.list)[i])|grepl("TB",names(cov.list)[i])){
    covs<-cov.list[[i]]
    
    pop<-gsub("_.+",'', names(cov.list)[i])
    
    covs<-covs[order(covs$cov12, decreasing=T),]
    n<-ceiling(nrow(covs)*0.01) #top1% region
    covs12_top<-covs[1:n,c(1:4)]
    covs12_top<-covs12_top[order(covs12_top$chrom, covs12_top$start),]
    covs12_top$window<-'100k'
    covs12_top$pop<-pop
    cov12<-rbind(cov12, covs12_top)
    
    covs<-covs[order(covs$cov13, decreasing=T),]
    n<-ceiling(nrow(covs)*0.01) #top1% region
    covs13_top<-covs[1:n,c(1:3,6)]
    covs13_top<-covs13_top[order(covs13_top$chrom, covs13_top$start),]
    covs13_top$window<-'100k'
    covs13_top$pop<-pop
    cov13<-rbind(cov13, covs13_top)
    
    covs<-covs[order(covs$cov23, decreasing=T),]
    n<-ceiling(nrow(covs)*0.01) #top1% region
    covs23_top<-covs[1:n,c(1:3,5)]
    covs23_top<-covs23_top[order(covs23_top$chrom, covs23_top$start),]
    covs23_top$window<-'100k'
    covs23_top$pop<-pop
    cov23<-rbind(cov23, covs23_top)
 }
 if (grepl("SS",names(cov.list)[i])){
    covs<-cov.list[[i]]
    
    pop<-gsub("_.+",'', names(cov.list)[i])
     
    covs<-covs[order(covs$cov23, decreasing=T),]
    n<-ceiling(nrow(covs)*0.01) #top1% region
    covs23_top<-covs[1:n,c(1:4)]
    covs23_top<-covs23_top[order(covs23_top$chrom, covs23_top$start),]
    covs23_top$window<-'100k'
    covs23_top$pop<-pop
    cov23<-rbind(cov23, covs23_top)
    
    }
}

write.csv(cov12, "../Output/COV/3pops_top1percent_outlier_regions.cov12_100k.csv",row.names = F)
write.csv(cov23, "../Output/COV/3pops_top1percent_outlier_regions.cov23_100k.csv",row.names = F)
write.csv(cov13, "../Output/COV/3pops_top1percent_outlier_regions.cov13_100k.csv",row.names = F)
#Create plots with different colors for outliers
#for COV12 and COV13 for TB and PWS (100K)
cv<-c("cov12","cov13","cov23")
winsize<-"100k"

for (i in 1:length(cv)){
    if (i==1|i==2){
            #PWS
            df1<-cov.list[[paste0("PWS_", winsize)]]
            df1<-df1[order(df1[,cv[i]], decreasing=T),]
            n<-ceiling(nrow(df1)*0.01) #top1% region
            df1$top1<-"N"
            df1$top1[1:n]<-"PWS"
            
            #tb
            df2<-cov.list[[paste0("TB_", winsize)]]
            df2<-df2[order(df2[,cv[i]], decreasing=T),]
            df2$top1<-"N"
            df2$top1[1:n]<-"TB"
            
            co<-rbind(df1, df2)
    
            co$chrom<-factor(co$chrom, levels=paste0("chr", 1:26))
            co$top1<-factor(co$top1, levels=c("PWS","TB","N"))
            colnames(co)[which(colnames(co)==cv[i])]<-"cov"
    
            ymax<-max(co$cov, na.rm=T)
            ggplot(co, aes(x=start/1000000, y=cov, color=top1))+
                geom_point(size=0.5)+
                facet_wrap(~chrom, ncol=4)+
                theme_classic()+ylim(-0.1,ymax)+
                scale_color_manual(values=c("deeppink","orange" ,"#ADD8E680"), labels=c("PWS", "TB", ""))+
                ylab("Covariance")+xlab('Postion (Mb)')+
                ggtitle(paste0(winsize[w]," window ",cv[i]))+
                scale_x_continuous(labels = comma)+
                guides(color = guide_legend(override.aes = list(color=c("deeppink","orange","white")), title=element_text("Top 1%")))
   
                ggsave(paste0("../Output/COV/COVscan_3pop/3Pops.",cv[i],"_perChrom_",winsize, "Window_Outliers.png"), width = 10, height = 8, dpi=300)
       
    }
   
    if (i==3){
        
        #pws
        df1<-cov.list[[paste0("PWS_", winsize[w])]]
        df1<-df1[,c("chrom","start","end","cov23")]
        df1<-df1[order(df1$cov23, decreasing=T),]
        n<-ceiling(nrow(df1)*0.01) #top1% region
        df1$top1<-"N"
        df1$top1[1:n]<-"PWS"
    
        #tb
        df2<-cov.list[[paste0("TB_", winsize)]]
        df2<-df2[,c("chrom","start","end","cov23")]
        df2<-df2[order(df2$cov23, decreasing=T),]
        df2$top1<-"N"
        df2$top1[1:n]<-"TB"
    
        #ss
        df3<-cov.list[[paste0("SS_", winsize)]
        df3<-df3[,c("chrom","start","end","cov23")]
        df3<-df3[order(df3$cov23, decreasing=T),]
        df3$top1<-"N"
        df3$top1[1:n]<-"SS"

        co<-rbind(df1,df2,df3)

        co$chrom<-factor(co$chrom, levels=paste0("chr", 1:26))
        co$top1<-factor(co$top1, levels=c("PWS","TB","SS","N"))
        ymax<-max(co$cov23, na.rm=T)
        ggplot(co, aes(x=start/1000000, y=cov23, color=top1))+
            geom_point(size=0.5)+
            facet_wrap(~chrom, ncol=4)+
            theme_classic()+ylim(-0.1,ymax)+
            ylab("Covariance")+xlab('Postion (Mb)')+
            ggtitle(paste0(winsize[w]," window ",cv[i]))+
            scale_x_continuous(labels = comma)+
            #scale_color_discrete(breaks=c("PWS","SS","TB"))+
            scale_color_manual(values=c("deeppink","orange",gre,"#ADD8E666"), labels=c("PWS","TB","SS", ""))+
            guides(color = guide_legend(override.aes = list(color=c("deeppink","orange",gre,         "white")),title=element_text("Top 1% outliers"))) 
        ggsave(paste0("../Output/COV/COVscan_3pop/3Pops.cov23_3Pops_perChrom_",winsize, "Window_Outliers.png"), width = 10, height = 9, dpi=300)
        }
    }
}

2.2.1 Whole genome plots without facet wrap

cv<-c("cov12","cov13","cov23")
for (i in 1:length(cv)){
    if (i==1|i==2){
        df1<-cov.list[["PWS_100k"]]
        df1<-df1[order(df1[,cv[i]], decreasing=T),]
        n<-ceiling(nrow(df1)*0.01) #top1% region
        df1$top1<-"N"
        df1$top1[1:n]<-"PWS"
        
        df2<-cov.list[["TB_100k"]]
        df2<-df2[order(df2[,cv[i]], decreasing=T),]
        df2$top1<-"N"
        df2$top1[1:n]<-"TB"

        co<-rbind(df1, df2)
        co$chrom<-factor(co$chrom, levels=paste0("chr", 1:26))
        colnames(co)[which(colnames(co)==cv[i])]<-"cov"
        #assgin colors
        co$top1<-apply(co, 1, function(x) {ifelse (x['top1']=="N", x['color'], x['top1'])} )
        co$top1<-factor(co$top1, levels=c("PWS","TB","col1","col2"))
        
        #count the number of sites per chromosomes
        poss<-data.frame(chr=paste0("chr",1:26))
        k=1
        for (j in 1:26){
            df<-df1[df1$chr==paste0("chr",j),]
            poss$start[j]<-k
            poss$end[j]<-k+nrow(df)-1
            k=k+nrow(df)
        }
        poss$x<-poss$start+(poss$end-poss$start)/2
        ymax<-max(co$cov, na.rm=T)
        ggplot(co, aes(x=index, y=cov, color=top1))+
            geom_point(size=0.5)+
            theme_classic()+ylim(-0.1,ymax)+
            scale_color_manual(values=c("#FF3293B3","#0433FFB3" ,"#ADD8E680","#C0C0C088"), labels=c("PWS", "TB", "",""))+
                ylab("Covariance")+
                ggtitle(paste0(" 100k window ",cv[i]))+
                guides(color = guide_legend(override.aes = list(color=c("deeppink","#0433FF","white","white"), size=2), title=element_text("Outlier (1%)", size=10)))+
            scale_x_continuous(name="Chromosome", breaks=poss$x, labels=1:26)
        ggsave(paste0("../Output/COV/COVscan_3pop/3Pops.",cv[i],"_100k_Window_Outliers.png"), width = 10, height = 3.5, dpi=300)
        }}
    if (i==3){
        #pws
        df1<-cov.list[["PWS_100k"]]
        df1<-df1[,c("chrom","start","end","cov23","index","color")]
        df1<-df1[order(df1$cov23, decreasing=T),]
        n<-ceiling(nrow(df1)*0.01) #top1% region
        df1$top1<-df1$color
        df1$top1[1:n]<-"PWS"
        #tb
        df2<-cov.list[["TB_100k"]]
        df2<-df2[,c("chrom","start","end","cov23","index","color")]
        df2<-df2[order(df2$cov23, decreasing=T),]
        df2$top1<-df2$color
        df2$top1[1:n]<-"TB"
        #ss
        df3<-cov.list[["SS_100k"]]
        df3<-df3[,c("chrom","start","end","cov23","index","color")]
        df3<-df3[order(df3$cov23, decreasing=T),]
        df3$top1<-df3$color
        df3$top1[1:n]<-"SS"

        co<-rbind(df1,df2,df3)
        co$top1<-factor(co$top1, levels=c("PWS","TB","SS","col1","col2"))
        
        #count the number of sites per chromosomes
        poss<-data.frame(chr=paste0("chr",1:26))
        k=1
        for (j in 1:26){
            df<-df1[df1$chr==paste0("chr",j),]
            poss$start[j]<-k
            poss$end[j]<-k+nrow(df)-1
            k=k+nrow(df)
        }
        poss$x<-poss$start+(poss$end-poss$start)/2
        ymax<-max(co$cov, na.rm=T)
        ggplot(co, aes(x=index, y=cov23, color=top1))+
            geom_point(size=0.5)+
            theme_classic()+ylim(-0.1,ymax)+
            scale_color_manual(values=c("#FF3293B3","#0433FFB3","#008F00B3" ,"#ADD8E680","#C0C0C088"), labels=c("PWS", "TB","SS", "",""))+
                ylab("Covariance")+
                ggtitle(paste0(" 100k window ",cv[i]))+
                guides(color = guide_legend(override.aes = list(color=c("deeppink","#0433FF","#008F00","white","white"), size=2), title=element_text("Outlier (1%)")))+
            scale_x_continuous(name="Chromosome", breaks=poss$x, labels=1:26)+
            theme(legend.title = element_text(size=10))
        ggsave(paste0("../Output/COV/COVscan_3pop/3Pops.",cv[i],"_100k_Window_Outliers.png"), width = 10, height = 3.5, dpi=300)
        }
    }
}

## Plot 3 time periods together for PWS and TB
Cov<-data.frame()
for (i in 1:length(cv)){
    df1<-cov.list[["PWS_100k"]]
    df1<-df1[order(df1[,cv[i]], decreasing=T),]
    n<-ceiling(nrow(df1)*0.01) #top1% region
    df1$top1<-"N"
    df1$top1[1:n]<-"PWS"
    
    df2<-cov.list[["TB_100k"]]
    df2<-df2[order(df2[,cv[i]], decreasing=T),]
    df2$top1<-"N"
    df2$top1[1:n]<-"TB"
    co<-rbind(df1, df2)
    co$chrom<-factor(co$chrom, levels=paste0("chr", 1:26))
    colnames(co)[which(colnames(co)==cv[i])]<-"cov"
    #assgin colors
    co$top1<-apply(co, 1, function(x) {ifelse (x['top1']=="N", x['color'], x['top1'])} )
    co$top1<-factor(co$top1, levels=c("PWS","TB","col1","col2"))
    co$time<-cv[i]
    
    Cov<-rbind(Cov, co[,c("index", "cov","top1","time")])
}

#count the number of sites per chromosomes
df1<-cov.list[["PWS_100k"]]
poss<-data.frame(chr=paste0("chr",1:26))
k=1
for (j in 1:26){
        df<-df1[df1$chr==paste0("chr",j),]
        poss$start[j]<-k
        poss$end[j]<-k+nrow(df)-1
        k=k+nrow(df)
}
poss$x<-poss$start+(poss$end-poss$start)/2
ymax<-max(co$cov, na.rm=T)
ggplot(Cov, aes(x=index, y=cov, color=top1))+
    facet_wrap(~time, ncol=1)+
    geom_point(size=0.5)+
    theme_classic()+ylim(-0.1,ymax)+
    scale_color_manual(values=c("#FF3293B3","#0433FFB3" ,"#ADD8E680","#C0C0C088"), labels=c("PWS", "TB", "",""))+
    ylab("Covariance")+
    guides(color = guide_legend(override.aes = list(color=c("deeppink","#0433FF","white","white"), size=2), title=element_text("Outlier (1%)", size=10)))+
    scale_x_continuous(name="Chromosome", breaks=poss$x, labels=1:26)

ggsave(paste0("../Output/COV/COVscan_3pop/PWS_TB_100k_Window_Outliers.png"), width = 11, height = 5, dpi=300)
        }}

2.3 Overlapping outlier regions between different populations

#100k
cv<-c("cov12","cov13","cov23")
pairs<-t(combn(pops, 2))
pairs<-data.frame(pairs)
colnames(pairs)<-paste0("pop",1:2)
Ov_direct<-data.frame(cov=c(cv[1:2],"cov23-PT","cov23-PS","cov23-ST" ,"cov23-3"))
Ov_300<-data.frame(cov=c(cv[1:2],"cov23-PT","cov23-PS","cov23-ST" ,"cov23-3"))
for (i in 1:length(cv)){
    df<-read.csv(paste0("../Output/COV/COVscan_3pop/3pops_top1percent_outlier_regions.", cv[i], "_new.csv"))
    df<-df[df$window=="100k",]
    df$id<-paste0(df$chrom,"_",df$start)
    
    if (i!=3){
        #exact overlaps
        isec<-intersect(df$id[df$pop=="PWS"], df$id[df$pop=="TB"]) 
        Ov_direct$count[i]<-length(isec)
        
        #### Check chromosome region overlap +-200,000 bases
        pop1<-df[df$pop=="PWS",]
        pop2<-df[df$pop=="TB",]
        overlps<-data.frame()
        overlps2<-data.frame()
        for (n in 1: nrow(pop1)){
            re<-pop2[pop2$chrom==pop1$chrom[n],]
            if (nrow(re)>=1){
                for (s in 1: nrow(re)){
                    if (re$start[s]<=pop1$start[n]+300000 & re$start[s]>=pop1$start[n]-300000){
                        overlps<-rbind(overlps, re[s,])
                        overlps2<-rbind(overlps2,pop1[n,])}
                }
            }
        }
        # Merge two tables into one summary overlap table:
        ov<-data.frame(id=overlps$id)
        for (n in 1: nrow(overlps)){
            if (overlps$start[n]<overlps2$start[n]) {ov$start[n]<-overlps$start[n]; ov$end[n]<-overlps2$end[n]}
            if (overlps$start[n]>=overlps2$start[n]) {ov$start[n]<-overlps2$start[n];ov$end[n]<-overlps$end[n]}
        }
        ov[,"cov.PWS"]<-overlps[,4]
        ov[,"cov.TB"]<-overlps2[,4]
        write.csv(ov, paste0("../Output/COV/COVscan_3pop/Overlap_regions_",cv[i],"_plusminus300k.csv"), row.names = F)
        Ov_300$count[i]<-nrow(ov)
        }
        
    if (i==3){
        isec<-intersect(df$id[df$pop=="PWS"], df$id[df$pop=="TB"]) 
        isec2<-intersect(df$id[df$pop=="PWS"], df$id[df$pop=="SS"]) 
        isec3<-intersect(df$id[df$pop=="SS"], df$id[df$pop=="TB"]) 
        Ov_direct$count[i]<-length(isec)
        Ov_direct$count[i+1]<-length(isec2)
        Ov_direct$count[i+2]<-length(isec)
        Ov_direct$count[i+3]<-length(intersect(df$id[df$pop=="SS"], intersect(df$id[df$pop=="PWS"], df$id[df$pop=="TB"])))
        
        for(j in 1:nrow(pairs)){
        #### Check chromosome region overlap +-200,000 bases
            pop1<-df[df$pop==pairs[j,1],]
            pop2<-df[df$pop==pairs[j,2],]
            overlps<-data.frame()
            overlps2<-data.frame()
            for (n in 1: nrow(pop1)){
                re<-pop2[pop2$chrom==pop1$chrom[n],]
                if (nrow(re)>=1){
                    for (s in 1: nrow(re)){
                        if (re$start[s]<=pop1$start[n]+300000 & re$start[s]>=pop1$start[n]-300000){
                            overlps<-rbind(overlps, re[s,])
                            overlps2<-rbind(overlps2,pop1[n,])}
                    }
                }
            }
        # Merge two tables into one summary overlap table:
            ov<-data.frame(id=overlps$id)
            for (n in 1: nrow(overlps)){
                if (overlps$start[n]<overlps2$start[n]) {ov$start[n]<-overlps$start[n]; ov$end[n]<-overlps2$end[n]}
                if (overlps$start[n]>=overlps2$start[n]) {ov$start[n]<-overlps2$start[n];ov$end[n]<-overlps$end[n]}
            }
        
            ov[,paste0("cov.",pairs[j,1])]<-overlps[,4]
            ov[,paste0("cov.",pairs[j,2])]<-overlps2[,4]
            ov<-ov[!duplicated(ov),]
            write.csv(ov, paste0("../Output/COV/COVscan_3pop/Overlap_regions_",cv[i],"_",pairs[j,1],".", pairs[j,2],"_plusminus300k.csv"), row.names = F)
            Ov_300$count[i+j-1]<-nrow(ov)
    }
    }
}
write.csv(Ov_direct, paste0("../Output/COV/COVscan_3pop/DDirect_Overlapping_regions_counts_3pop_summary.csv"))
Ov_300$count[6]<-NA
write.csv(Ov_300, paste0("../Output/COV/COVscan_3pop/Overlapping_regions_counts_3pop_plusMinus300k.csv"))

3 Run the snpEff pipeline to find annotation in the outlier regions (100k-window+-100k)

3.1 Create a script to run SnpEff

Create VCF files with selected regions & run snpEff

#Create bed files
cv<-c("cov12","cov13","cov23")
#Prevent scientific notation in bed files
options(scipen=999)

#The first line of bed files is often not red by vcftools
for (i in 1:3){
    df<-read.csv(paste0("../Output/COV/COVscan_3pop/3pops_top1percent_outlier_regions.",cv[i],"_new.csv"))
    #select 100k
    df<-df[df$window=="100k",]
    #add 100k
    df$start<-df$start-100000
    df$end<-df$end+100000
    dfp<-df[df$pop=="PWS",1:3]
    colnames(dfp)<-c('track type=bedGraph', '1','1')
    write.table(dfp, paste0("../Output/COV/COVscan_3pop/PWS_outliers_",cv[i],"_new.bed"),quote = F, row.names = F, col.names = T,sep = "\t")
    dft<-df[df$pop=="TB",1:3]
    colnames(dft)<-c('track type=bedGraph', '1','1')
    write.table(dft, paste0("../Output/COV/COVscan_3pop/TB_outliers_",cv[i],"_new.bed"),quote = F, row.names = F, col.names = F,sep = "\t")
    
    if (i==3){
        dfs<-df[df$pop=="SS",1:3]
        colnames(dfs)<-c('track type=bedGraph', '1','1')
        write.table(dfs, paste0("../Output/COV/COVscan_3pop/SS_outliers_",cv[i],"_new.bed"),quote = F, row.names = F, col.names = F,sep = "\t")
    }
}

# Create a bash script to run snpEff
bedfiles<-list.files("../Output/COV/COVscan_3pop/", pattern="*_new.bed")

sink("../COVscan_createVCFs_3Pops2.sh")
cat("#!/bin/bash \n\n")
for (i in 1:length(bedfiles)){
    fname<-gsub(".bed",'', bedfiles[i])
    cat(paste0("vcftools --gzvcf Data/new_vcf/3pop/3pops.MD7000_NS0.5_maf05.vcf.gz --bed Output/COV/COVscan_3pop/", bedfiles[i], " --out Output/COV/COVscan_3pop/", fname," --recode --keep-INFO-all \n"))
}
sink(NULL)  

#create a bash script to run snpEff
vfiles<-list.files("../Output/COV/COVscan_3pop/", pattern=".recode.vcf")

sink("~/programs/snpEff/runsnpEff_cov_3pop.sh")
cat("#!/bin/bash \n\n")
for (i in 1:length(vfiles)){
    fname<-gsub("_new.recode.vcf","",vfiles[i])
    cat(paste0("java -Xmx8g -jar snpEff.jar Ch_v2.0.2.99 ~/Projects/PacHerring/Output/COV/COVscan_3pop/",vfiles[i], " -stats ~/Projects/PacHerring/Output/COV/COVscan_3pop/",fname,".html >  ~/Projects/PacHerring/Output/COV/COVscan_3pop/Anno.",fname,".vcf \n"))
    
    #extract the annotation information
    cat(paste0("bcftools query -f '%CHROM %POS %INFO/AF %INFO/ANN\\n' ~/Projects/PacHerring/Output/COV/COVscan_3pop/Anno.",fname,".vcf > ~/Projects/PacHerring/Output/COV/COVscan_3pop/",fname,"_annotation \n\n"))

}
sink(NULL)  

3.2 Create summary gene files from snpEff and check overlapping genes.

## Create summary files of snpEff results (gene annotations in the regions of interest) and reformat as a ShinyGo input 

#create gene list 
gfiles<-list.files("../Output/COV/COVscan_3pop/", pattern="genes.txt")

for (i in 1:length(gfiles)){
    df<-read.table(paste0("../Output/COV/COVscan_3pop/",gfiles[i]), sep="\t")
    df<-df[,1:7]
    colnames(df)<-c("GeneName","GeneId","TranscriptId","BioType","variants_impact_HIGH","variants_impact_LOW",  "variants_impact_MODERATE")
    
    fname<-gsub(".genes.txt","",gfiles[i])
    genes<-unique(df$GeneId)
    sink(paste0("../Output/COV/COVscan_3pop/geneIDlist_",fname,".txt"))
    cat(paste0(genes,"; "))
    sink(NULL)
}

#Annotation infor from SnpEff
cv<-c("cov12","cov13","cov23")
for (c in 1:3){
    if (c!=3){
    for (p in 1:2){
        ano<-read.table(paste0("../Output/COV/COVscan_3pop/",pops[p],"_outliers_",cv[c],"_annotation"), header = F)
        annotations<-data.frame()
        for (i in 1: nrow(ano)){
            anns<-unlist(strsplit(ano$V4[i], "\\,|\\|"))
            annm<-data.frame(matrix(anns,ncol = 16, byrow = TRUE))
            annm<-annm[,c(2,3,4,5,8)]
            colnames(annm)<-c("Effect","Putative_impact","Gene_name","Gene_ID","Feature type")
            annm<-annm[!duplicated(annm), ]
            annm$chr<-ano$V1[i]
            annm$pos<-ano$V2[i]
            annm$AF<- ano$V3[i]
            annotations<-rbind(annotations, annm)
        }     
        annotations<-annotations[,c(6:8,1:5)]
        annotations<-annotations[!duplicated(annotations[,1:2]),]
        write.csv(annotations, paste0("../Output/COV/COVscan_3pop/Genes_",pops[p],"_outliers_100k_",cv[c],".csv"), row.names = F)
    }
    }
    if (c==3){
        for (p in 1:3){
        ano<-read.table(paste0("../Output/COV/COVscan_3pop/",pops[p],"_outliers_",cv[c],"_annotation"), header = F)
        annotations<-data.frame()
        for (i in 1: nrow(ano)){
            anns<-unlist(strsplit(ano$V4[i], "\\,|\\|"))
            annm<-data.frame(matrix(anns,ncol = 16, byrow = TRUE))
            annm<-annm[,c(2,3,4,5,8)]
            colnames(annm)<-c("Effect","Putative_impact","Gene_name","Gene_ID","Feature type")
            annm<-annm[!duplicated(annm), ]
            annm$chr<-ano$V1[i]
            annm$chr<-ano$V1[i]
            annm$pos<-ano$V2[i]
            annm$AF<- ano$V3[i]
            annotations<-rbind(annotations, annm)
        }     
        annotations<-annotations[,c(6:8,1:5)]
        annotations<-annotations[!duplicated(annotations[,1:2]),]
        write.csv(annotations, paste0("../Output/COV/COVscan_3pop/Genes_",pops[p],"_outliers_100k_",cv[c],".csv"), row.names = F)
    }
}
  
}

3.3 Find the overlapping gene names

gnamesfiles<-list.files("../Output/COV/COVscan_3pop/", pattern="Genes_.+outliers_100k.+\\d.csv$")

for (i in 1:length(gnamesfiles)){
    df<-read.csv(paste0("../Output/COV/COVscan_3pop/",gnamesfiles[i]))
    df<-df[,c(1,6:7)]
    df<-df[!duplicated(df),]
    
    fname<-gsub(".csv","", gnamesfiles[i])
    fname<-gsub("Genes_","", fname)
    
    
    #add gene names for front and back of intergenic regions
    df2<-df[grep("-", df$Gene_ID),]
    k=1
    df_div<-data.frame()
    oddnames<-data.frame()
    for (j in 1:nrow(df2)){
        names<-unlist(strsplit(df2$Gene_name[j], "-"))
        ids<-unlist(strsplit(df2$Gene_ID[j], "-"))
        
        if (length(names)==2){
            df_div<-rbind(df_div, c(df2$chr[j],names[1],ids[1]))
            k=k+1
            df_div<-rbind(df_div, c(df2$chr[j],names[2],ids[2]))
            k=k+1
        }
       
        if (length(names)!=2){
            n<-grep("si:", names)
            if (length(n)>0){
                if (n==1) newnames<-c(paste0(names[1],"-",names[2]), names[3])
                if (n==2) newnames<-c(names[1],paste0(names[2],"-",names[3]))
                df_div<-rbind(df_div, c(df2$chr[j],newnames[1],ids[1]))
                k=k+1
                df_div<-rbind(df_div, c(df2$chr[j],newnames[2],ids[2]))
                k=k+1
            }
            
            if (length(n)==0) {
                oddnames<-rbind(oddnames, df2[j,])
            }
        }
    }
    df_div<-df_div[!duplicated(df_div),]
    df_div<-df_div[df_div$Gene_ID!="CHR_END",]
    df_div<-df_div[df_div$Gene_ID!="CHR_START",]
    
    remove<-grep("-", df$Gene_ID)
    df<-df[-remove,]
    df<-rbind(df, df_div)
    df<-df[!duplicated(df),]
    
    if (nrow(oddnames)!=0){
        write.csv(df, paste0("../Output/COV/COVscan_3pop/",fname,"GeneList_withIntergenicGenes.csv" ), row.names = F)
        write.csv(oddnames, paste0("../Output/COV/COVscan_3pop/Oddnames_", fname,".csv"))
    }
    if (nrow(oddnames)==0){
        write.csv(df, paste0("../Output/COV/COVscan_3pop/",fname,"GeneList_withIntergenicGenes_new.csv" ), row.names = F)
     }
}
   

## !! ##
## Manually change the oddnames and add them eteo the GeneList files 
#(updated file names has "_new" at the end)

#aggregate all gene names
gnew<-list.files("../Output/COV/COVscan_3pop/", pattern="GeneList_withIntergenicGenes_new.csv$")
Genes<-data.frame()
GeneList<-list()
for (i in 1:length(gnew)){
    df<-read.csv(paste0("../Output/COV/COVscan_3pop/", gnew[i]))
    GeneList[[i]]<-df
    fname<-gsub("GeneList_withIntergenicGenes_new.csv",'',gnew[i])
    names(GeneList)[i]<-fname
    dup<-df[duplicated(df),]
    df<-df[!duplicated(df),]
    Genes<-rbind(Genes, df)
    Genes<-Genes[!duplicated(Genes),]
    
}


#1. Between populations
times<-c("cov12","cov13","cov23")
common<-list()
common_summary<-data.frame(time=times)
for (i in 1:3){
    tlist<-GeneList[grep(times[i], names(GeneList))]
    if (i !=3){
        common_genes<-intersect(tlist[[1]]["Gene_name"], tlist[[2]]["Gene_name"])
        common[[i]]<-common_genes
        names(common)[[i]]<-times[i]
        common_summary$PWS[i]<-nrow(tlist[[grep("PWS", names(tlist))]])
        common_summary$TB[i]<-nrow(tlist[[grep("TB", names(tlist))]])
        common_summary$SS[i]<-NA
        common_summary$common_PWS.TB[i]<-nrow(common_genes)
        
        pws<-tlist[[1]]["Gene_name"]
        tb<-tlist[[2]]["Gene_name"]
        x<-list(PWS=pws$Gene_name,TB=tb$Gene_name)
        ggvenn(x, fill_color = cols[c(2,1)], stroke_size = 0.5, set_name_size = 4,text_size=3)+ggtitle(times[i])
        ggsave(paste0("../Output/COV/COVscan_3pop/Venn_",times[i],".png"), width = 3, height=3, dpi=300)
    }
    if (i==3){
        common_summary$PWS[i]<-nrow(tlist[[grep("PWS", names(tlist))]])
        common_summary$TB[i]<- nrow(tlist[[grep("TB", names(tlist))]])
        common_summary$SS[i]<- nrow(tlist[[grep("SS", names(tlist))]])
        
        genes1<-intersect(tlist[[1]]["Gene_name"], tlist[[3]]["Gene_name"])
        genes2<-intersect(tlist[[1]]["Gene_name"], tlist[[2]]["Gene_name"])
        genes3<-intersect(tlist[[2]]["Gene_name"], tlist[[3]]["Gene_name"])
        genes4<-intersect(tlist[[1]]["Gene_name"],intersect(tlist[[2]]["Gene_name"], tlist[[3]]["Gene_name"]))
        common_summary$common_PWS.TB[i]<-nrow(genes1)
        common_summary$common_PWS.SS[i]<-nrow(genes2)
        common_summary$common_SS.TB[i]<-nrow(genes3)
        common_summary$common3[i]<-nrow(genes4)
        k=i
        common[[k]]<-genes2
        names(common)[[k]]<-paste0(times[i],"_PWS.SS")
        k=k+1
        common[[k]]<-genes1
        names(common)[[k]]<-paste0(times[i],"_PWS.TB")
        k=k+1
        common[[k]]<-genes3
        names(common)[[k]]<-paste0(times[i],"_SS.TB")
        k=k+1
        common[[k]]<-genes4
        names(common)[[k]]<-paste0(times[i],"_3pops")
        
        pws<-tlist[[1]]["Gene_name"]
        tb<-tlist[[3]]["Gene_name"]
        ss<-tlist[[2]]["Gene_name"]
        x<-list(PWS=pws$Gene_name,TB=tb$Gene_name, SS=ss$Gene_name)
        ggvenn(x, fill_color = cols[c(2,1,3)], stroke_size = 0.5, set_name_size = 4,text_size=3)+ggtitle(times[i])
        ggsave(paste0("../Output/COV/COVscan_3pop/Venn_",times[i],".png"), width = 4, height=4, dpi=300)
        
         x1<-list(PWS=pws$Gene_name,TB=tb$Gene_name)
        ggvenn(x1, fill_color = cols[c(2,1)], stroke_size = 0.5, set_name_size = 4,text_size=3)+ggtitle(times[i])
        ggsave(paste0("../Output/COV/COVscan_3pop/Venn_PWS_TB_",times[i],".png"), width = 3, height=3, dpi=300)
         x2<-list(PWS=pws$Gene_name,SS=ss$Gene_name)
        ggvenn(x2, fill_color = cols[c(2,3)], stroke_size = 0.5, set_name_size = 4,text_size=3)+ggtitle(times[i])
        ggsave(paste0("../Output/COV/COVscan_3pop/Venn_PWS_SS_",times[i],".png"), width = 3, height=3, dpi=300)
          x3<-list(SS=ss$Gene_name, TB=tb$Gene_name)
        ggvenn(x3, fill_color = cols[c(3,1)], stroke_size = 0.5, set_name_size = 4,text_size=3)+ggtitle(times[i])
        ggsave(paste0("../Output/COV/COVscan_3pop/Venn_SS_TB_",times[i],".png"), width = 3, height=3, dpi=300)
        
        
        }
}
write.csv(common_summary, "../Output/COV/COVscan_3pop/Common_genes_withIntergenes_3pops.csv")


#What are the overlapping gene names between populations
common_times<-list()
for (i in 1: length(common)){
    gids<-common[[i]]
    df<-data.frame(Gene_name=gids)
    
    df2<-merge(df, Genes, by="Gene_name")
    write.csv(df2, paste0("../Output/COV/COVscan_3pop/Common_genes_", names(common)[i],".csv"), row.names = F)
    common_times[[i]]<-df2
    names(common_times)[i]<- names(common)[i]
}

#non_overlapping genes COV23

tlist<-GeneList[grep(times[3], names(GeneList))]
pws23<-tlist[[1]]["Gene_ID"]
ss23<-tlist[[2]]["Gene_ID"]
tb23<-tlist[[3]]["Gene_ID"]
genes1<-intersect(tlist[[1]]["Gene_ID"], tlist[[3]]["Gene_ID"])
genes2<-intersect(tlist[[1]]["Gene_ID"], tlist[[2]]["Gene_ID"])
genes3<-intersect(tlist[[2]]["Gene_ID"], tlist[[3]]["Gene_ID"])
genes4<-intersect(tlist[[1]]["Gene_ID"],intersect(tlist[[2]]["Gene_ID"], tlist[[3]]["Gene_ID"]))

overlaps<-rbind(genes1, genes2, genes3, genes4)
overlaps<-overlaps[!duplicated(overlaps$Gene_ID),]
pws23_only<-data.frame(Gene_ID=pws23[!(pws23$Gene_ID %in% overlaps), ])
g1<-unique(pws23_only$Gene_ID)
sink("../Output/COV/COVscan_3pop/genes/GeneID_PWSonly_COV23.txt")
cat(paste0(g1, ";"))
sink(NULL)
        
tb23_only<-tb23[!(tb23$Gene_ID %in% overlaps), ]
g1<-unique(tb23_only)
sink("../Output/COV/COVscan_3pop/genes/GeneID_TBonly_COV23.txt")
cat(paste0(g1, ";"))
sink(NULL)

ss23_only<-ss23[!(ss23$Gene_ID %in% overlaps), ]
g1<-unique(ss23_only)
sink("../Output/COV/COVscan_3pop/genes/GeneID_SSonly_COV23.txt")
cat(paste0(g1, ";"))
sink(NULL)




#2. Between Time-Points within a population

times<-c("cov12","cov13","cov23")
pops<-c("PWS","TB")
common2<-list()
common_summary2<-data.frame(pop=rep(pops[1:2], each=4))
for (i in 1:length(pops)){
    plist<-GeneList[grep(pops[i], names(GeneList))]
    k=4*i-3
    #common genes between COV12 and COV13
    common_genes1<-intersect(plist[[1]]["Gene_name"], plist[[2]]["Gene_name"])
    common2[[k]]<-common_genes1
    names(common2)[[k]]<-paste0(pops[i],".", times[1],"_",times[2])
    common_summary2$Time[k]<-paste0(times[1],"_",times[2])
    common_summary2$no.of.genes[k]<-nrow(common_genes1) 
    
    c12<-plist[[1]]["Gene_name"]
    c13<-plist[[2]]["Gene_name"]
    c23<-plist[[3]]["Gene_name"]
    x<-list(COV12=c12$Gene_name,COV13=c13$Gene_name, COV23=c23$Gene_name)
    ggvenn(x, fill_color = cols[c(1,5,7)], stroke_size = 0.5, set_name_size = 4,text_size=3)+ggtitle(pops[i])
    ggsave(paste0("../Output/COV/COVscan_3pop/Venn_",pops[i],".png"), width = 4, height=4, dpi=300)

    
    k=k+1
    #common genes between COV12 and COV23
    common_genes2<-intersect(plist[[1]]["Gene_name"], plist[[3]]["Gene_name"])
    common2[[k]]<-common_genes2
    names(common2)[[k]]<-paste0(pops[i],".", times[1],"_",times[3])
    common_summary2$Time[k]<-paste0(times[1],"_",times[3])
    common_summary2$no.of.genes[k]<-nrow(common_genes2) 
 
    k=k+1
    #common genes between COV13 and COV23
    common_genes3<-intersect(plist[[2]]["Gene_name"], plist[[3]]["Gene_name"])
    common2[[k]]<-common_genes3
    names(common2)[[k]]<-paste0(pops[i],".", times[2],"_",times[3])
    common_summary2$Time[k]<-paste0(times[2],"_",times[3])
    common_summary2$no.of.genes[k]<-nrow(common_genes3) 
 
    k=k+1
    #common genes among all time periods
    common_genes4<-intersect(plist[[1]]["Gene_name"], (intersect(plist[[2]]["Gene_name"], plist[[3]]["Gene_name"])))
    common2[[k]]<-common_genes4
    names(common2)[[k]]<-paste0(pops[i],".all")
    common_summary2$Time[k]<-"All"
    common_summary2$no.of.genes[k]<-nrow(common_genes4) 
}
write.csv(common_summary2, "../Output/COV/COVscan_3pop/Common_genes_betweenTimePoints.csv")


#Common gene names between time points

for (i in 1:2){
    CommonGenes<-data.frame()
    glist<-common2[grep(pops[i], names(common2))]
    for(j in 1:length(glist)){
        gids<-glist[[j]]
        df<-data.frame(Gene_name=gids)
        df2<-merge(df, Genes, by="Gene_name", all.x=T)
        write.csv(df2, paste0("../Output/COV/COVscan_3pop/Common_genes_", names(glist)[j],".csv"), row.names = F)
        df2$Time<-names(glist)[j]
        CommonGenes<-rbind(CommonGenes, df2)
    }
    write.csv(CommonGenes, paste0("../Output/COV/COVscan_3pop/Common_genes_",pops[i] ,".csv"), row.names = F)
}

3.3.0.1 Overlapping gene numbers

time PWS TB SS common_PWS.TB common_PWS.SS common_SS.TB common3
cov12 568 635 NA 89 NA NA NA
cov13 586 643 NA 80 NA NA NA
cov23 495 628 417 64 80 30 10

3.4 What are the genes overlapping across different time points between populations?

#1. between PWS and TB
pws.tb<-common_times[c(1,2,4)]

genes1213<-intersect(pws.tb[[1]]["Gene_name"], pws.tb[[2]]["Gene_name"])
genes1213<-merge(genes1213, Genes, by="Gene_name")
write.csv(genes1223, "../Output/COV/COVscan_3pop/Common_genes_PWS.TB.cov12-cov23.csv")
#           Gene_name   chr            Gene_ID
#1 ENSCHAG00000001687 chr13 ENSCHAG00000001687
#2             ndst2a chr13 ENSCHAG00000002649
#3             zswim8 chr13 ENSCHAG00000005956

## Common genes between population across time points (COV12 - COV13)
p1213<-pws.tb[[1]]["Gene_name"]
t1213<-pws.tb[[2]]["Gene_name"]
x<-list(PWS=p1213$Gene_name,TB=t1213$Gene_name)
ggvenn(x, fill_color = cols[c(2,1)], stroke_size = 0.5, set_name_size = 4,text_size=3)+ggtitle("COV12-COV13 in PWS & TB")
ggsave(paste0("../Output/COV/COVscan_3pop/Venn_PWS_TB_COV12-COV13.png"), width = 3, height=3, dpi=300)
        
genes1223<-intersect(pws.tb[[1]]["Gene_name"], pws.tb[[3]]["Gene_name"])
genes1223<-merge(genes1223, Genes, by="Gene_name")
write.csv(genes1223, "../Output/COV/COVscan_3pop/Common_genes_PWS.TB.cov12-cov13.csv")
#           Gene_name   chr            Gene_ID
#1 ENSCHAG00000001687 chr13 ENSCHAG00000001687
#2 ENSCHAG00000022709 chr20 ENSCHAG00000022709
#3 ENSCHAG00000022815 chr20 ENSCHAG00000022815
#4             ndst2a chr13 ENSCHAG00000002649

p1223<-pws.tb[[1]]["Gene_name"]
t1223<-pws.tb[[3]]["Gene_name"]
x<-list(PWS=p1223$Gene_name,TB=t1223$Gene_name)
ggvenn(x, fill_color = cols[c(2,1)], stroke_size = 0.5, set_name_size = 4,text_size=3)+ggtitle("COV12-COV23 in PWS & TB")
ggsave(paste0("../Output/COV/COVscan_3pop/Venn_PWS_TB_COV12-COV23.png"), width = 3, height=3, dpi=300)


genes1323<-intersect(pws.tb[[2]]["Gene_name"], pws.tb[[3]]["Gene_name"])
genes1323<-merge(genes1323, Genes, by="Gene_name")
write.csv(genes1323, "../Output/COV/COVscan_3pop/Common_genes_PWS.TB.cov13-cov23.csv")
#           Gene_name   chr            Gene_ID
#1 ENSCHAG00000001687 chr13 ENSCHAG00000001687
#2             ndst2a chr13 ENSCHAG00000002649

p1323<-pws.tb[[2]]["Gene_name"]
t1323<-pws.tb[[3]]["Gene_name"]
x<-list(PWS=p1323$Gene_name,TB=t1323$Gene_name)
ggvenn(x, fill_color = cols[c(2,1)], stroke_size = 0.5, set_name_size = 4,text_size=3)+ggtitle("COV13-COV23 in PWS & TB")
ggsave(paste0("../Output/COV/COVscan_3pop/Venn_PWS_TB_COV13-COV23.png"), width = 3, height=3, dpi=300)

4 Compare results from PH_MD7000 (all pops together) and 3Pops-NS0.5 (filter for 3 pops) VCF files

pws1<-read.csv("../Output/COV/COVscan_3pop/3pops_top1percent_outlier_regions.cov12_new.csv")
pws1<-pws1[pws1$pop=="PWS" & pws1$window=="100k",]
pws2<-read.csv("../Output/COV/3pops_top1percent_outlier_regions.cov12.csv")
pws2<-pws2[pws2$pop=="PWS"&pws2$window=="100k",]

pws1$vcf<-"3Pops"
pws2$vcf<-"PH"
pws<-rbind(pws1[,c("chrom","start","end","cov12","vcf")],pws2[,c("chrom","start","end","cov12","vcf")])
pws$chrom<-factor(pws$chrom, levels=paste0("chr", 1:26))
ggplot(data=pws,aes(x=start/1000000, y=cov12, color=vcf, fill=vcf))+
    facet_wrap(~chrom)+
    geom_point(position=position_dodge(width = 0.7), alpha=0.5)+xlab("Position (Mb)")+
    ggtitle("PWS COV12")
ggsave("../Output/COV/COVscan_3pop/PWS_PH.vs.3Pops_outlier_overlap_cov12.png", width = 10, height = 8, dpi=300)

pws1<-read.csv("../Output/COV/COVscan_3pop/3pops_top1percent_outlier_regions.cov23_new.csv")
pws1<-pws1[pws1$pop=="PWS" & pws1$window=="100k",]
pws2<-read.csv("../Output/COV/3pops_top1percent_outlier_regions.cov23.csv")
pws2<-pws2[pws2$pop=="PWS"&pws2$window=="100k",]

pws1$vcf<-"PH"
pws2$vcf<-"3Pops"
pws<-rbind(pws1[,c("chrom","start","end","cov23","vcf")],pws2[,c("chrom","start","end","cov23","vcf")])
pws$chrom<-factor(pws$chrom, levels=paste0("chr", 1:26))
ggplot(data=pws,aes(x=start/1000000, y=cov23, color=vcf, fill=vcf))+
    facet_wrap(~chrom)+
    geom_point(position=position_dodge(width = 0.7), alpha=0.5)+xlab("Position (Mb)")+
    ggtitle("PWS COV23")
ggsave("../Output/COV/COVscan_3pop/PWS_PH.vs.3Pops_outlier_overlap_cov23.png", width = 10, height = 8, dpi=300)

pws1<-read.csv("../Output/COV/COVscan_3pop/3pops_top1percent_outlier_regions.cov13_new.csv")
pws1<-pws1[pws1$pop=="PWS" & pws1$window=="100k",]
pws2<-read.csv("../Output/COV/3pops_top1percent_outlier_regions.cov13.csv")
pws2<-pws2[pws2$pop=="PWS"&pws2$window=="100k",]

pws1$vcf<-"PH"
pws2$vcf<-"3Pops"
pws<-rbind(pws1[,c("chrom","start","end","cov13","vcf")],pws2[,c("chrom","start","end","cov13","vcf")])
pws$chrom<-factor(pws$chrom, levels=paste0("chr", 1:26))
ggplot(data=pws,aes(x=start/1000000, y=cov13, color=vcf, fill=vcf))+
    facet_wrap(~chrom)+
    geom_point(position=position_dodge(width = 0.7), alpha=0.5)+xlab("Position (Mb)")+
    ggtitle("PWS COV13")
ggsave("../Output/COV/COVscan_3pop/PWS_PH.vs.3Pops_outlier_overlap_cov13.png", width = 10, height = 8, dpi=300)

5 Interpopulation comparison per time period

### Interpopulation comparisons
#decode the samples to create the right matrix
cv<-read.csv("~/Projects/Pacherring_Vincent/MD7000/GW_Covs_interPopulations_100k_3pops.csv", header = F)
labs<-read.csv("~/Projects/Pacherring_Vincent/MD7000/GW_Covs_interPopulations_100k_labels_3pops.csv" )
labs<-labs[,-1]
cvm<-data.frame(label=as.vector(t(labs)), cov=as.vector(t(cv)))

#rearrange based on comparions: covariance between populations within the same period
#PopYr Symbols
# PH 1 'PWS', 1991
# PH 2 'PWS', 1996
# PH 3 'PWS', 2006
# PH 4 'PWS', 2017
# PH 5 'SS',  1991
# PH 6 'SS',  1996
# PH 7 'SS',  2006
# PH 8 'SS',  2017
# PH 9 'TB',  1991
# PH 10'TB',  1996
# PH 11'TB',  2006
# PH 12'TB',  2017

Covs<-data.frame(pops=rep(c("PWS.vs.SS", "PWS.vs.TB",  "SS.vs.TB"), times=6),
                 period=c(rep("1991-1996", times=3),rep("1996-2006", times=3), rep("2006-2017", times=3)))

Covs$cov<-c(NA, cvm$cov[cvm$label=="cov(PH: 2-1, PH: 10-9)"],NA,
            cvm$cov[cvm$label=="cov(PH: 3-2, PH: 7-6)"],cvm$cov[cvm$label=="cov(PH: 3-2, PH: 11-10)"], 
            cvm$cov[cvm$label=="cov(PH: 7-6, PH: 11-10)"],
            cvm$cov[cvm$label=="cov(PH: 4-3, PH: 8-7)"],cvm$cov[cvm$label=="cov(PH: 4-3, PH: 12-11)"],cvm$cov[cvm$label=="cov(PH: 8-7, PH: 12-11)"])


#C.I.
cis<-read.csv("~/Projects/Pacherring_Vincent/MD7000/GW_COV_Interpop_comparison_CIs.csv")
cis<-cis[,-1]
cim<-data.frame(label=as.vector(t(labs)), ci_l=as.vector(t(cis[1:11,])))
cim$ci_h<-as.vector(t(cis[12:22,]))

Covs$ci_l<-as.numeric(c(NA,cim$ci_l[cim$label=="cov(PH: 2-1, PH: 10-9)"],NA,
                      cim$ci_l[cim$label=="cov(PH: 3-2, PH: 7-6)"],cim$ci_l[cim$label=="cov(PH: 3-2, PH: 11-10)"], cim$ci_l[cim$label=="cov(PH: 7-6, PH: 11-10)"],
                      cim$ci_l[cim$label=="cov(PH: 4-3, PH: 8-7)"],cim$ci_l[cim$label=="cov(PH: 4-3, PH: 12-11)"], cim$ci_l[cim$label=="cov(PH: 8-7, PH: 12-11)"]))

Covs$ci_h<-as.numeric(c(NA, cim$ci_h[cim$label=="cov(PH: 2-1, PH: 10-9)"],NA,
                      cim$ci_h[cim$label=="cov(PH: 3-2, PH: 7-6)"],cim$ci_h[cim$label=="cov(PH: 3-2, PH: 11-10)"], cim$ci_h[cim$label=="cov(PH: 7-6, PH: 11-10)"],
                      cim$ci_h[cim$label=="cov(PH: 4-3, PH: 8-7)"],cim$ci_h[cim$label=="cov(PH: 4-3, PH: 12-11)"], cim$ci_h[cim$label=="cov(PH: 8-7, PH: 12-11)"]))

#Barplot
ggplot(Covs, aes(x=period, y=cov, fill=pops))+
    geom_bar(stat="identity",position=position_dodge(width = 0.7), width=0.8)+
    ylab("Covariance")+xlab('')+theme_classic()+
    geom_hline(yintercept = 0,color="gray70", size=0.3)+
    scale_fill_manual(values=cols[c(2,1,3)])+
    theme(legend.title = element_blank())+
    scale_y_continuous(labels = comma)+
    ylim(-0.0013, 0.002)+
    geom_errorbar(aes(ymin=ci_l, ymax=ci_h), width=.2, size=.2, position=position_dodge(width = 0.7))
ggsave("../Output/COV/Interpop_cov_comparison_3Pops_new.png",width = 4.7, height = 3, dpi=300)

#Point plot
ggplot(Covs, aes(x=period, y=cov, color=pops))+
    geom_point(position=position_dodge(width = 0.7), size=4)+
    ylab("Covariance")+xlab('')+theme_classic()+
    geom_hline(yintercept = 0,color="gray70", size=0.3)+
    scale_color_manual(values=cols[c(2,1,3)])+
    theme(legend.title = element_blank())+
    scale_y_continuous(labels = comma)+
    geom_errorbar(aes(ymin=ci_l, ymax=ci_h), width=.2, size=.2, position=position_dodge(width = 0.7))+
    ylim(-0.0013, 0.002)+
     geom_vline(xintercept = c(1.5,2.5), color="gray", size=0.3)
ggsave("../Output/COV/Interpop_cov_comparison_3Pops_new_PointPlot.png",width = 4.7, height = 3, dpi=300)

#line plot
Covs$time<-1
Covs$time[Covs$period=="1996-2006"]<-2
Covs$time[Covs$period=="2006-2017"]<-3
Covs<-Covs[order(Covs$time),]
ggplot(Covs, aes(x=time, y=cov, color=pops, group=pops))+
    geom_point(position=position_dodge(width = 0.7), size=4)+
    geom_path(position=position_dodge(width = 0.7))+
    ylab("Covariance")+xlab('')+theme_classic()+
    geom_hline(yintercept = 0,color="gray70", size=0.3)+
    scale_color_manual(values=cols[c(2,1,3)])+
    theme(legend.title = element_blank())+
    scale_y_continuous(labels = comma)+
    geom_errorbar(aes(ymin=ci_l, ymax=ci_h), width=.2, size=.2, position=position_dodge(width = 0.7))+
    ylim(-0.0013, 0.002)+
     geom_vline(xintercept = c(1.5,2.5), color="gray", size=0.3)+
    scale_x_continuous(breaks=c(1,2,3), labels = c("1991-1996","1996-2006","2006-2017"))
ggsave("../Output/COV/Interpop_cov_comparison_3Pops_new_LinePlot.png",width = 4.7, height = 3, dpi=300)

5.1 Longer time period

## Longer time-period
Covs2<-data.frame(pops=rep(c("PWS.vs.SS", "PWS.vs.TB",  "SS.vs.TB"), times=3),
                 period=c(rep("1991-2006", times=3),rep("1991-2017", times=3),rep("1996-2017", times=3)))

cv1<-read.csv("~/Projects/Pacherring_Vincent/MD7000/GW_Covs_interPopulations_100k_1991-2006.csv", header = F)
labs1<-read.csv("~/Projects/Pacherring_Vincent/MD7000/GW_Covs_interPopulations_100k_labels_1991-2006.csv" )
labs1<-labs1[,-1]
cvm1<-data.frame(label=as.vector(t(labs1)), cov=as.vector(t(cv1)))

cv2<-read.csv("~/Projects/Pacherring_Vincent/MD7000/GW_Covs_interPopulations_100k_1991-2017.csv", header = F)
labs2<-read.csv("~/Projects/Pacherring_Vincent/MD7000/GW_Covs_interPopulations_100k_labels_1991-2017.csv" )
labs2<-labs2[,-1]
cvm2<-data.frame(label=as.vector(t(labs2)), cov=as.vector(t(cv2)))

cv3<-read.csv("~/Projects/Pacherring_Vincent/MD7000/GW_Covs_interPopulations_100k_1996-2017.csv", header = F)
labs3<-read.csv("~/Projects/Pacherring_Vincent/MD7000/GW_Covs_interPopulations_100k_labels_1996-2017.csv" )
labs3<-labs3[,-1]
cvm3<-data.frame(label=as.vector(t(labs3)), cov=as.vector(t(cv3)))

Covs2$cov<-c(NA, cvm1$cov[cvm1$label=="cov(PH: 2-1, PH: 4-3)"], NA,
             NA, cvm2$cov[cvm2$label=="cov(PH: 2-1, PH: 4-3)"], NA,
             cvm3$cov[cvm3$label=="cov(PH: 2-1, PH: 4-3)"], cvm3$cov[cvm3$label=="cov(PH: 2-1, PH: 6-5)"], cvm3$cov[cvm3$label=="cov(PH: 4-3, PH: 6-5)"])

#C.I.
cis1<-read.csv("~/Projects/Pacherring_Vincent/MD7000/GW_COV_Interpop_comparison_CIs_1991-2006.csv")
cis1<-cis1[,-1]
cis2<-read.csv("~/Projects/Pacherring_Vincent/MD7000/GW_COV_Interpop_comparison_CIs_1991-2017.csv")
cis2<-cis2[,-1]
cis3<-read.csv("~/Projects/Pacherring_Vincent/MD7000/GW_COV_Interpop_comparison_CIs_1996-2017.csv")
cis3<-cis3[,-1]

#cim<-data.frame(label=as.vector(t(labs)), ci_l=as.vector(t(cis1[1:4,])))
#cim$ci_h<-as.vector(t(cis[12:22,]))

Covs2$ci_l<-as.numeric(c(NA,cis1[1,3],NA,
                        NA,cis2[1,3],NA,
                      cis3[1,3],cis3[1,5],cis3[3,5]))

Covs2$ci_h<-as.numeric(c(NA,cis1[4,3],NA,
                        NA,cis2[4,3],NA,
                      cis3[6,3],cis3[6,5],cis3[8,5]))


ggplot(Covs2, aes(x=period, y=cov, fill=pops))+
    geom_bar(stat="identity",position=position_dodge(width = 0.7), width=0.8)+
    ylab("Covariance")+xlab('')+theme_classic()+
    geom_hline(yintercept = 0,color="gray70", size=0.3)+
    scale_fill_manual(values=cols[c(2,1,3)])+
    theme(legend.title = element_blank())+scale_y_continuous(labels = comma)+
    ylim(-0.0013, 0.002)+
    geom_errorbar(aes(ymin=ci_l, ymax=ci_h), width=.2, size=.2, position=position_dodge(width = 0.7))
ggsave("../Output/COV/Interpop_cov_comparison_3Pops_LonogerPeriod.png",width = 4.7, height = 3, dpi=300)

ggplot(Covs2, aes(x=period, y=cov, color=pops))+
    geom_point(position=position_dodge(width = 0.7), size=4)+
    ylab("Covariance")+xlab('')+theme_classic()+
    geom_hline(yintercept = 0,color="gray70", size=0.3)+
    scale_color_manual(values=cols[c(2,1,3)])+
    theme(legend.title = element_blank())+
    scale_y_continuous(labels = comma)+
     ylim(-0.0013, 0.002)+
    geom_errorbar(aes(ymin=ci_l, ymax=ci_h), width=.2, size=.2, position=position_dodge(width = 0.7))+
    geom_vline(xintercept = c(1.5,2.5), color="gray", size=0.3)
ggsave("../Output/COV/Interpop_cov_comparison_3PopsLonogerPeriod_PointPlot.png",width = 4.7, height = 3, dpi=300)

6 Focused freq analysis

6.1 rel gene (chr13: 23070000 - 23080000)

pops<-c("PWS91","PWS96","PWS07","PWS17")
yr<-c(1991,1996,2007,2017)
maf<-data.frame()
for (i in 1:4){
    af<-read.table(paste0("../Data/new_vcf/AF/",pops[i],".mafs"),sep="\t", header = T)
    af<-af[af$chromo=="chr13"&af$position>=23070000&af$position<=23080000,]
    af$year<-yr[i]
    maf<-rbind(maf,af)
}
#write.csv(maf,"../Output/COV/COVscan_3pop/AF_maf_chr13_23Mb.csv")

ggplot(maf, aes(x=year, y=knownEM, color=factor(position)))+
    geom_point()+
    geom_path()+ggtitle("MAF (ANGSD)")+
    ylab("maf")+
    theme(legend.title=element_blank())
ggsave("../Output/COV/COVscan/AF_ch13_23.07-23.08_angsd.png", width = 5, height=3, dpi=300)


AF<-data.frame()
for (i in 1:4){
    af<-read.table(paste0("../Data/new_vcf/AF/",pops[i],"_maf05_af.frq"),header = FALSE, skip=1, col.names = c("chr","pos","n_allele","n_sample","MajorAF","MAF"))
    af<-af[af$chr=="chr13"&af$pos>=23070000&af$pos<=23080000,]
    af$year<-yr[i]
    af$maf<-substr(af$MAF, 3,10)
    af$maf<-as.numeric(af$maf)
    af<-af[,c("chr","pos","maf","year")]
    AF<-rbind(AF,af)
}
ggplot(AF, aes(x=year, y=maf, color=factor(pos)))+
    geom_point()+
    geom_path()+ggtitle("MAF (vcftools)")+
    theme(legend.title=element_blank())
ggsave("../Output/COV/COVscan/AF_ch13_23.07-23.08.png", width = 5, height=3, dpi=300)


###TB
pops<-c("TB91","TB96","TB06","TB17","PWS91","PWS96","PWS07","PWS17","SS96","SS06","SS17")
yr<-c(1991,1996,2006,2017,1991,1996,2007,2017,1996,2006,2017)
maf<-data.frame()
for (i in 1:length(pops)){
    af<-read.table(paste0("../Data/new_vcf/AF/",pops[i],".mafs"),sep="\t", header = T)
    af<-af[af$chromo=="chr13"&af$position>=23070000&af$position<=23080000,]
    af$year<-yr[i]
    af$pop<-sub("\\d\\d","", pops[i])
    maf<-rbind(maf,af)
}
#write.csv(maf,"../Output/COV/COVscan_3pop/AF_maf_chr13_23Mb.csv")
ggplot(maf, aes(x=year, y=knownEM, color=pop))+
    facet_wrap(~factor(position))+
    geom_point()+
    geom_path()+ggtitle("Chr13 (rel gene)")+
    ylab("maf")+
    theme(legend.title=element_blank())
ggsave("../Output/COV/COVscan/AF_ch13_23.07-23.08_angsd.png", width = 5, height=3, dpi=300)

LS0tCnRpdGxlOiAiQ09WIHNjYW4gMyBwb3B1bGF0aW9ucyIKb3V0cHV0OgogIGh0bWxfbm90ZWJvb2s6CiAgICAgIHRvYzogdHJ1ZSAKICAgICAgdG9jX2Zsb2F0OiB0cnVlCiAgICAgIG51bWJlcl9zZWN0aW9uczogdHJ1ZQogICAgICB0aGVtZTogbHVtZW4KICAgICAgaGlnaGxpZ2h0OiB0YW5nbwogICAgICBjb2RlX2ZvbGRpbmc6IGhpZGUKICAgICAgZGZfcHJpbnQ6IHBhZ2VkCi0tLQoKYGBge3IgZXZhbD1GQUxTRSwgbWVzc2FnZT1GQUxTRSwgd2FybmluZz1GQUxTRSwgaW5jbHVkZT1GQUxTRX0Kc291cmNlKCIuLi9Sc2NyaXB0cy9CYXNlU2NyaXB0cy5SIikKbGlicmFyeSh0aWR5dmVyc2UpCmxpYnJhcnkoZHBseXIpCmxpYnJhcnkoY293cGxvdCkKbGlicmFyeShzY2FsZXMpCmxpYnJhcnkoZ2d2ZW5uKQpgYGAKCiMgUGxvdCB0aGUgY292YXJpYW5jZXMKCmBgYHtyIGV2YWw9RkFMU0UsIG1lc3NhZ2U9RkFMU0UsIHdhcm5pbmc9RkFMU0V9CnBvcHM8LWMoIlBXUyIsIlRCIiwiU1MiKQpjb3ZzPC1kYXRhLmZyYW1lKCkKVmFyaWFuY2U8LWRhdGEuZnJhbWUoKQp3aW5zaXplPC1jKCI1MGsiLCI3NWsiLCIxMDBrIiwiMjUwayIpCgpmb3IgKHcgaW4gMTogbGVuZ3RoKHdpbnNpemUpKXsKICAgIGNvdnM8LWRhdGEuZnJhbWUoKQogICAgZm9yIChwIGluIDE6IGxlbmd0aChwb3BzKSl7CiAgICAgICAgI2NvdmFyaWFuY2Ugb3V0cHV0IGZpbGUKICAgICAgICBjb3Y8LXJlYWQuY3N2KHBhc3RlMCgifi9Qcm9qZWN0cy9QYWNoZXJyaW5nX1ZpbmNlbnQvTUQ3MDAwLzNQb3BzX21hZjA1X3RlbXBfY292X21hdHJpeF8iLHBvcHNbcF0sIl8iLHdpbnNpemVbd10sIi5jc3YiKSkKICAgICAgICBjb3Y8LWNvdlssLTFdCiAgICAgICAgCiAgICAgICAgY2k8LXJlYWQuY3N2KHBhc3RlMCgifi9Qcm9qZWN0cy9QYWNoZXJyaW5nX1ZpbmNlbnQvTUQ3MDAwLzNQb3BzX21hZjA1XyIscG9wc1twXSwiX0Nvdl9DSXNfYm9vdHN0cmFwNTAwMF8iLHdpbnNpemVbd10sIndpbmRvdy5jc3YiKSkKICAgICAgICBjaTwtY2lbLC0xXQogICAgICAgIAogICAgICAgICNyZXNoYXBlIHRoZSBtYXRyaXgKICAgICAgICBtYXQxPC1jb3ZbMTozLF0KICAgICAgICBtYXQyPC1jb3ZbNDo2LF0KICAgICAgICAKICAgICAgICBjb3ZkZjwtZGF0YS5mcmFtZSgpCiAgICAgICAgaz0xCiAgICAgICAgZm9yIChpIGluIDE6bnJvdyhtYXQxKSl7CiAgICAgICAgICAgIGZvciAoaiBpbiAxOm5jb2wobWF0MSkpewogICAgICAgICAgICAgICAgY292ZGZbaywxXTwtbWF0MltpLGpdCiAgICAgICAgICAgICAgICBjb3ZkZltrLDJdPC1tYXQxW2ksal0KICAgICAgICAgICAgICAgIGs9aysxCiAgICAgICAgICAgIH0KICAgICAgICB9CiAgICAgICAgY29sbmFtZXMoY292ZGYpPC1jKCJsYWJlbCIsInZhbHVlIikKICAgICAgICBjb3ZkZiR2YWx1ZTwtYXMubnVtZXJpYyhjb3ZkZiR2YWx1ZSkKICAgICAgICBjb3ZhcjwtY292ZGZbZ3JlcCgiY292Iixjb3ZkZiRsYWJlbCksXQogICAgICAgIHZhcnM8LWNvdmRmW2dyZXAoInZhciIsY292ZGYkbGFiZWwpLF0KICAgICAgICAKICAgICAgICAjcmVtb3ZlIHRoZSByZWR1bmRhbnQgdmFsdWVzCiAgICAgICAgaWYgKHBvcHNbcF0hPSJTUyIpIGNvdmFyPC1jb3ZhclshZHVwbGljYXRlZChjb3ZhclssIDJdKSxdIAogICAgICAgIGlmIChwb3BzW3BdPT0iU1MiKSBjb3ZhcjwtY292YXJbYygxLDIsNCksXQogICAgICAgIAogICAgICAgICNhc3NpZ24gdGhlIHN0YXJ0aW5nIHRpbWUgcGVyaW9kIGFuZCBjb3ZlcmluZyBwZXJpb2QgdmFsdWVzCiAgICAgICAgY292YXIkeWVhcjwtYygxLDIsMikKICAgICAgICBjb3ZhciRzZXJpZXM8LWMoIjE5OTEiLCIxOTkxIiwiMTk5NiIpCiAgICAgICAgCiAgICAgICAgdmFycyR5ZWFyPC1jKDEsMiwyKQogICAgICAgIHZhcnMkc2VyaWVzPC1jKCIxOTkxIiwiMTk5MSIsIjE5OTYiKQogICAgICAgIAogICAgICAgICNhc3NpZ24gcG9wdWxhdGlvbiBuYW1lCiAgICAgICAgY292YXIkbG9jYXRpb248LXBvcHNbcF0KICAgICAgICB2YXJzJGxvY2F0aW9uPC1wb3BzW3BdCiAgICAgICAgCiAgICAgICAgI2F0dGFjaCBjaSBpbmZvCiAgICAgICAgY292YXIkY2lfbDwtYyhjaVsxLDJdLCBjaVsxLDNdLGNpWzIsM10pCiAgICAgICAgY292YXIkY2lfdTwtYyhjaVs0LDJdLCBjaVs0LDNdLGNpWzUsM10pCiAgICAgICAgCiAgICAgICAgI2NvbWJpbmUgaW4gdG8gb25lIG1hdHJpeAogICAgICAgIGNvdnM8LXJiaW5kKGNvdnMsIGNvdmFyKQogICAgICAgIFZhcmlhbmNlPC1yYmluZChWYXJpYW5jZSwgdmFycykKICAgIH0KICAgIAogICAgY292cyRjaV9sPC1hcy5udW1lcmljKGNvdnMkY2lfbCkKICAgIGNvdnMkY2lfdTwtYXMubnVtZXJpYyhjb3ZzJGNpX3UpCiAgICBnZ3Bsb3QoZGF0YT1jb3ZzLCBhZXMoeD15ZWFyLCB5PXZhbHVlLCBjb2xvcj1sb2NhdGlvbiwgc2hhcGU9c2VyaWVzLCBncm91cD1pbnRlcmFjdGlvbihsb2NhdGlvbiwgc2VyaWVzKSkpKwogICAgICAgIGdlb21fcG9pbnQoc2l6ZT0zLCBwb3NpdGlvbj1wb3NpdGlvbl9kb2RnZSh3aWR0aCA9IDAuMSxwcmVzZXJ2ZSA9InRvdGFsIikpKwogICAgICAgICNnZW9tX2Vycm9yYmFyKGRhdGE9Y292cywgYWVzKHg9eWVhciwgeT12YWx1ZSwgeW1pbj1jaV9sLCB5bWF4PWNpX3UpLCB3aWR0aD0uMiwgc2l6ZT0uMiwgcG9zaXRpb249cG9zaXRpb25fZG9kZ2Uod2lkdGggPSAwLjEscHJlc2VydmUgPSJ0b3RhbCIpKSsKICAgICAgICBnZW9tX2xpbmUoZGF0YT1jb3ZzLCBhZXMoeD15ZWFyLCB5PXZhbHVlLGNvbG9yPWxvY2F0aW9uLCBncm91cD1pbnRlcmFjdGlvbihsb2NhdGlvbiwgc2VyaWVzKSksIHBvc2l0aW9uPXBvc2l0aW9uX2RvZGdlKHdpZHRoID0gMC4xLHByZXNlcnZlID0idG90YWwiKSkrCiAgICAgICAgeWxhYigiQ292YXJpYW5jZSIpK3hsYWIoJycpK3RoZW1lX2NsYXNzaWMoKSsKICAgICAgICB0aGVtZShheGlzLnRleHQueCA9IGVsZW1lbnRfYmxhbmsoKSxsZWdlbmQudGl0bGUgPSBlbGVtZW50X2JsYW5rKCkpKwogICAgICAgIGdlb21faGxpbmUoeWludGVyY2VwdCA9IDAsY29sb3I9ImdyYXk3MCIsIHNpemU9MC4zKSsKICAgICAgICBnZW9tX2Vycm9yYmFyKGFlcyh5bWluPWNpX2wsIHltYXg9Y2lfdSksIHdpZHRoPS4yLCBzaXplPS4yLCBwb3NpdGlvbj1wb3NpdGlvbl9kb2RnZSh3aWR0aCA9IDAuMSxwcmVzZXJ2ZSA9InRvdGFsIikpKwogICAgICAgIHNjYWxlX3NoYXBlX21hbnVhbCh2YWx1ZXM9YygxNiwxNyksbGFiZWxzPWMoIjE5OTEtIiwiMTk5Ni0iKSkrCiAgICAgICAgc2NhbGVfeF9jb250aW51b3VzKGJyZWFrcyA9IGMoMSwyKSkKICAgIGdnc2F2ZShwYXN0ZTAoIi4uL091dHB1dC9DT1YvM1BvcHNfQ292X292ZXJ0aW1lX0NJXyIsd2luc2l6ZVt3XSwiLndpbmRvdy5wbmciKSx3aWR0aCA9IDQuNywgaGVpZ2h0ID0gMywgZHBpPTMwMCkKfSAgICAKCmBgYAohW10oLi4vT3V0cHV0L0NPVi8zUG9wc19Db3Zfb3ZlcnRpbWVfQ0lfMTAway53aW5kb3cucG5nKQoKCiMgRmluZCByZWdpb25zIHdpdGggaGlnaCBjb3ZhcmlhbmNlcyBpbiBlYWNoIHBvcHVsYXRpb24KKiBGcm9tIFRlbXBvcmFsIENvdmFyaWFuY2UgYW5hbHlzaXMgIC1vdXRwdXQgY292YXJpYW5jZXMgZm9yIGVhY2ggdGltZSBwZXJpb2QKCiMjIFBsb3QgdGhlIGNvdmFyaWFuY2VzIGFjcm9zcyB0aGUgZ2Vub21lICAKCmBgYHtyIGV2YWw9RkFMU0UsIG1lc3NhZ2U9RkFMU0UsIHdhcm5pbmc9RkFMU0V9CgojRmluZCB0aGUgcmVnaW9ucyB3aXRoIGEgaGlnaCB0ZW1wb3JhbCBjb3ZhcmlhbmNlIApwb3BzPC1jKCJQV1MiLCJUQiIsIlNTIikKd2luc2l6ZTwtIjEwMGsiCmV2ZW5zPC1wYXN0ZTAoImNociIsc2VxKDIsMjYsIGJ5PTIpKQpjb3YubGlzdDwtbGlzdCgpCmNvdnNfYWxsPC1saXN0KCkKaz0xCmZvciAocCBpbiAxOiBsZW5ndGgocG9wcykpewogICAgcG9wPC1wb3BzW3BdCiAgICBpdjwtcmVhZC5jc3YocGFzdGUwKCJ+L1Byb2plY3RzL1BhY2hlcnJpbmdfVmluY2VudC9NRDcwMDAvM3BvcHNfaW50ZXJ2YWxzXyIsd2luc2l6ZSwid2luZG93LmNzdiIpLCByb3cubmFtZXMgPSAxKQogICAgaWYgKHA9PTMpIHsKICAgICAgICBjb3YyMzwtcmVhZC5jc3YocGFzdGUwKCJ+L1Byb2plY3RzL1BhY2hlcnJpbmdfVmluY2VudC9NRDcwMDAvIixwb3AsIl9jb3YyM18yMDE3LTIwMDZfMjAwNi0xOTk2XzNQb3BzXyIsd2luc2l6ZSwid2luZG93LmNzdiIpLCBoZWFkZXIgPSBGKQogICAgICAgIGNvdnM8LWNiaW5kKGl2LCBjb3YyMykKICAgICAgICBjb2xuYW1lcyhjb3ZzKVs0XTwtYygiY292MjMiKQogICAgICAgIGNvdnMkaW5kZXg9MTpucm93KGNvdnMpCiAgICAgICAgY292cyRjb2xvcjwtImNvbDEiCiAgICAgICAgY292cyRjb2xvcltjb3ZzJGNocm9tICVpbiUgZXZlbnNdPC0iY29sMiIKCiAgICAgICAgY292c1tzYXBwbHkoY292cywgaXMuaW5maW5pdGUpXSA8LSBOQQogICAgICAgIGNvdnNbc2FwcGx5KGNvdnMsIGlzLm5hbildIDwtIE5BCiAgICAgICAgCiAgICAgICAgY292Lmxpc3RbW2tdXTwtY292cwogICAgICAgIG5hbWVzKGNvdi5saXN0KVtrXTwtcGFzdGUwKHBvcCwiXyIsd2luc2l6ZSkgICAgCiAgICAgICAgaz1rKzEKICAgICAgICAgICAgCiAgICAgICAgeTwtbWluKGNvdnMkY292MjMsIG5hLnJtPVQpCiAgICAgICAgeW1pbjwtaWZlbHNlICh5PD0tMC4xLC0wLjEsIHkpIAogICAgICAgIHltYXg8LW1heChjb3ZzJGNvdjIzLCBuYS5ybT1UKQogICAgICAgIGdncGxvdChjb3ZzLCBhZXMoeD1pbmRleCwgeT1jb3YyMywgY29sb3I9Y29sb3IpKSsKICAgICAgICAgICAgZ2VvbV9wb2ludChzaXplPTEsIGFscGhhPTAuNSkrCiAgICAgICAgICAgIHRoZW1lX2NsYXNzaWMoKSsKICAgICAgICAgICAgeWxpbSh5bWluLHltYXgpKwogICAgICAgICAgICBzY2FsZV9jb2xvcl9tYW51YWwodmFsdWVzPWMoImdyYXk3MCIsInN0ZWVsYmx1ZSIpLCBndWlkZT0ibm9uZSIpKwogICAgICAgICAgICB5bGFiKCJDb3ZhcmlhbmNlIikreGxhYignQ2hyb21vc29tZScpKwogICAgICAgICAgICB0aGVtZShheGlzLnRleHQueCA9IGVsZW1lbnRfYmxhbmsoKSkrCiAgICAgICAgICAgIGdndGl0bGUocGFzdGUwKHBvcCwiICIsIHdpbnNpemUsIiB3aW5kb3ciKSkKICAgICAgICAjZ2dzYXZlKHBhc3RlMCgiLi4vT3V0cHV0L0NPVi8zUG9wcy4iLHBvcCwiX3RlbXBDb3ZzX2Fjcm9zc0dlbm9tZV8iLHdpbnNpemVbaV0sICJXaW5kb3cucG5nIiksIHdpZHRoID0gOCwgaGVpZ2h0ID0gMi43LCBkcGk9MzAwKSAKICAgICAgICB9CiAgICBlbHNlIHsKICAgICAgICBjb3YxMjwtcmVhZC5jc3YocGFzdGUwKCJ+L1Byb2plY3RzL1BhY2hlcnJpbmdfVmluY2VudC9NRDcwMDAvIixwb3AsIl9jb3YxMl8xOTk2LTE5OTFfMjAwNi0xOTk2XzNQb3BzXyIsd2luc2l6ZSwid2luZG93LmNzdiIpLCBoZWFkZXIgPSBGKQogICAgICAgIGNvdjIzPC1yZWFkLmNzdihwYXN0ZTAoIn4vUHJvamVjdHMvUGFjaGVycmluZ19WaW5jZW50L01ENzAwMC8iLHBvcCwiX2NvdjIzXzIwMTctMjAwNl8yMDA2LTE5OTZfM1BvcHNfIix3aW5zaXplLCJ3aW5kb3cuY3N2IiksIGhlYWRlciA9IEYpCiAgICAgICAgY292MTM8LXJlYWQuY3N2KHBhc3RlMCgifi9Qcm9qZWN0cy9QYWNoZXJyaW5nX1ZpbmNlbnQvTUQ3MDAwLyIscG9wLCJfY292MTNfMjAxNy0yMDA2XzE5OTYtMTk5MV8zUG9wc18iLHdpbnNpemUsIndpbmRvdy5jc3YiKSwgaGVhZGVyID0gRikKICAgICAgICBjb3ZzPC1jYmluZChpdiwgY292MTIsIGNvdjIzLGNvdjEzKQogICAgICAgIGNvbG5hbWVzKGNvdnMpWzQ6Nl08LWMoImNvdjEyIiwiY292MjMiLCJjb3YxMyIpCiAgICAgICAgY292cyRpbmRleD0xOm5yb3coY292cykKICAgIAogICAgICAgIGNvdnMkY29sb3I8LSJjb2wxIgogICAgICAgIGNvdnMkY29sb3JbY292cyRjaHJvbSAlaW4lIGV2ZW5zXTwtImNvbDIiCiAgICAKICAgICAgICBjb3ZzW3NhcHBseShjb3ZzLCBpcy5pbmZpbml0ZSldIDwtIE5BCiAgICAgICAgY292c1tzYXBwbHkoY292cywgaXMubmFuKV0gPC0gTkEKICAgICAgICAKICAgICAgICBjb3YubGlzdFtba11dPC1jb3ZzCiAgICAgICAgbmFtZXMoY292Lmxpc3QpW2tdPC1wYXN0ZTAocG9wLCJfIix3aW5zaXplKSAgICAKICAgICAgICBrPWsrMQogICAgICAgIGNvdnNtPC1tZWx0KGNvdnNbLGMoImluZGV4IiwiY29sb3IiLCJjb3YxMiIsImNvdjIzIiwiY292MTMiKV0sIGlkLnZhcnMgPSBjKCJpbmRleCIsICJjb2xvciIpKQogICAgICAgIHltYXg8LW1heChjb3ZzbSR2YWx1ZSwgbmEucm09VCkKICAgICAgICB5PC1taW4oY292c20kdmFsdWUsIG5hLnJtPVQpCiAgICAgICAgeW1pbjwtaWZlbHNlICh5PD0tMC4xLC0wLjEsIHkpIAogICAgICAgIGdncGxvdChjb3ZzbSwgYWVzKHg9aW5kZXgsIHk9dmFsdWUsIGNvbG9yPWNvbG9yKSkrCiAgICAgICAgICAgIGZhY2V0X3dyYXAofnZhcmlhYmxlLCBucm93PTMpKwogICAgICAgICAgICBnZW9tX3BvaW50KHNpemU9MSwgYWxwaGE9MC41KSsKICAgICAgICAgICAgdGhlbWVfY2xhc3NpYygpKwogICAgICAgICAgICB5bGltKHltaW4seW1heCkrCiAgICAgICAgICAgIHNjYWxlX2NvbG9yX21hbnVhbCh2YWx1ZXM9YygiZ3JheTcwIiwic3RlZWxibHVlIiksIGd1aWRlPSJub25lIikrCiAgICAgICAgICAgIHlsYWIoIkNvdmFyaWFuY2UiKSt4bGFiKCdDaHJvbW9zb21lJykrCiAgICAgICAgICAgIHRoZW1lKGF4aXMudGV4dC54ID0gZWxlbWVudF9ibGFuaygpKSsKICAgICAgICAgICAgZ2d0aXRsZShwYXN0ZTAocG9wLCIgIiwgd2luc2l6ZSwiIHdpbmRvdyIpKQogICAgICAgICNnZ3NhdmUocGFzdGUwKCIuLi9PdXRwdXQvQ09WLzNQb3BzLiIscG9wLCJfdGVtcENvdnNfYWNyb3NzR2Vub21lXyIsd2luc2l6ZSwgIldpbmRvdy5wbmciKSwgd2lkdGggPSA4LCBoZWlnaHQgPSA4LCBkcGk9MzAwKSAgICAKICAgIH0KfQpgYGAKCiFbXSguLi9PdXRwdXQvQ09WLzNQb3BzLlBXU190ZW1wQ292c19hY3Jvc3NHZW5vbWVfMTAwa1dpbmRvdy5wbmcpe3dpZHRoPTY1JX0KCiFbXSguLi9PdXRwdXQvQ09WLzNQb3BzLlRCX3RlbXBDb3ZzX2Fjcm9zc0dlbm9tZV8xMDBrV2luZG93LnBuZyl7d2lkdGg9NjUlfSAgIAoKIVtdKC4uL091dHB1dC9DT1YvM1BvcHMuU1NfdGVtcENvdnNfYWNyb3NzR2Vub21lXzEwMGtXaW5kb3cucG5nKXt3aWR0aD02NSV9ICAKCgojIyBGaW5kIHRoZSBvdXRsaWVyIHJlZ2lvbnMgZm9yIGVhY2ggdGltZSBwZXJpb2QgIAoKYGBge3IgZXZhbD1GQUxTRSwgbWVzc2FnZT1GQUxTRSwgd2FybmluZz1GQUxTRX0KCiNmaW5kIGhvdyBvdXRsaWVycyBvdmVybGFwIGJldHdlZW4gZGlmZmVyZW50IHdpbmRvd3MKY292MTI8LWRhdGEuZnJhbWUoKQpjb3YyMzwtZGF0YS5mcmFtZSgpCmNvdjEzPC1kYXRhLmZyYW1lKCkKCmZvciAoaSBpbiAxOmxlbmd0aChjb3YubGlzdCkpewogaWYgKGdyZXBsKCJQV1MiLG5hbWVzKGNvdi5saXN0KVtpXSl8Z3JlcGwoIlRCIixuYW1lcyhjb3YubGlzdClbaV0pKXsKICAgIGNvdnM8LWNvdi5saXN0W1tpXV0KICAgIAogICAgcG9wPC1nc3ViKCJfLisiLCcnLCBuYW1lcyhjb3YubGlzdClbaV0pCiAgICAKICAgIGNvdnM8LWNvdnNbb3JkZXIoY292cyRjb3YxMiwgZGVjcmVhc2luZz1UKSxdCiAgICBuPC1jZWlsaW5nKG5yb3coY292cykqMC4wMSkgI3RvcDElIHJlZ2lvbgogICAgY292czEyX3RvcDwtY292c1sxOm4sYygxOjQpXQogICAgY292czEyX3RvcDwtY292czEyX3RvcFtvcmRlcihjb3ZzMTJfdG9wJGNocm9tLCBjb3ZzMTJfdG9wJHN0YXJ0KSxdCiAgICBjb3ZzMTJfdG9wJHdpbmRvdzwtJzEwMGsnCiAgICBjb3ZzMTJfdG9wJHBvcDwtcG9wCiAgICBjb3YxMjwtcmJpbmQoY292MTIsIGNvdnMxMl90b3ApCiAgICAKICAgIGNvdnM8LWNvdnNbb3JkZXIoY292cyRjb3YxMywgZGVjcmVhc2luZz1UKSxdCiAgICBuPC1jZWlsaW5nKG5yb3coY292cykqMC4wMSkgI3RvcDElIHJlZ2lvbgogICAgY292czEzX3RvcDwtY292c1sxOm4sYygxOjMsNildCiAgICBjb3ZzMTNfdG9wPC1jb3ZzMTNfdG9wW29yZGVyKGNvdnMxM190b3AkY2hyb20sIGNvdnMxM190b3Akc3RhcnQpLF0KICAgIGNvdnMxM190b3Akd2luZG93PC0nMTAwaycKICAgIGNvdnMxM190b3AkcG9wPC1wb3AKICAgIGNvdjEzPC1yYmluZChjb3YxMywgY292czEzX3RvcCkKICAgIAogICAgY292czwtY292c1tvcmRlcihjb3ZzJGNvdjIzLCBkZWNyZWFzaW5nPVQpLF0KICAgIG48LWNlaWxpbmcobnJvdyhjb3ZzKSowLjAxKSAjdG9wMSUgcmVnaW9uCiAgICBjb3ZzMjNfdG9wPC1jb3ZzWzE6bixjKDE6Myw1KV0KICAgIGNvdnMyM190b3A8LWNvdnMyM190b3Bbb3JkZXIoY292czIzX3RvcCRjaHJvbSwgY292czIzX3RvcCRzdGFydCksXQogICAgY292czIzX3RvcCR3aW5kb3c8LScxMDBrJwogICAgY292czIzX3RvcCRwb3A8LXBvcAogICAgY292MjM8LXJiaW5kKGNvdjIzLCBjb3ZzMjNfdG9wKQogfQogaWYgKGdyZXBsKCJTUyIsbmFtZXMoY292Lmxpc3QpW2ldKSl7CiAgICBjb3ZzPC1jb3YubGlzdFtbaV1dCiAgICAKICAgIHBvcDwtZ3N1YigiXy4rIiwnJywgbmFtZXMoY292Lmxpc3QpW2ldKQogICAgIAogICAgY292czwtY292c1tvcmRlcihjb3ZzJGNvdjIzLCBkZWNyZWFzaW5nPVQpLF0KICAgIG48LWNlaWxpbmcobnJvdyhjb3ZzKSowLjAxKSAjdG9wMSUgcmVnaW9uCiAgICBjb3ZzMjNfdG9wPC1jb3ZzWzE6bixjKDE6NCldCiAgICBjb3ZzMjNfdG9wPC1jb3ZzMjNfdG9wW29yZGVyKGNvdnMyM190b3AkY2hyb20sIGNvdnMyM190b3Akc3RhcnQpLF0KICAgIGNvdnMyM190b3Akd2luZG93PC0nMTAwaycKICAgIGNvdnMyM190b3AkcG9wPC1wb3AKICAgIGNvdjIzPC1yYmluZChjb3YyMywgY292czIzX3RvcCkKICAgIAogICAgfQp9Cgp3cml0ZS5jc3YoY292MTIsICIuLi9PdXRwdXQvQ09WLzNwb3BzX3RvcDFwZXJjZW50X291dGxpZXJfcmVnaW9ucy5jb3YxMl8xMDBrLmNzdiIscm93Lm5hbWVzID0gRikKd3JpdGUuY3N2KGNvdjIzLCAiLi4vT3V0cHV0L0NPVi8zcG9wc190b3AxcGVyY2VudF9vdXRsaWVyX3JlZ2lvbnMuY292MjNfMTAway5jc3YiLHJvdy5uYW1lcyA9IEYpCndyaXRlLmNzdihjb3YxMywgIi4uL091dHB1dC9DT1YvM3BvcHNfdG9wMXBlcmNlbnRfb3V0bGllcl9yZWdpb25zLmNvdjEzXzEwMGsuY3N2Iixyb3cubmFtZXMgPSBGKQoKYGBgCgoKCmBgYHtyIGV2YWw9RkFMU0UsIG1lc3NhZ2U9RkFMU0UsIHdhcm5pbmc9RkFMU0V9CiNDcmVhdGUgcGxvdHMgd2l0aCBkaWZmZXJlbnQgY29sb3JzIGZvciBvdXRsaWVycwojZm9yIENPVjEyIGFuZCBDT1YxMyBmb3IgVEIgYW5kIFBXUyAoMTAwSykKY3Y8LWMoImNvdjEyIiwiY292MTMiLCJjb3YyMyIpCndpbnNpemU8LSIxMDBrIgoKZm9yIChpIGluIDE6bGVuZ3RoKGN2KSl7CiAgICBpZiAoaT09MXxpPT0yKXsKICAgICAgICAgICAgI1BXUwogICAgICAgICAgICBkZjE8LWNvdi5saXN0W1twYXN0ZTAoIlBXU18iLCB3aW5zaXplKV1dCiAgICAgICAgICAgIGRmMTwtZGYxW29yZGVyKGRmMVssY3ZbaV1dLCBkZWNyZWFzaW5nPVQpLF0KICAgICAgICAgICAgbjwtY2VpbGluZyhucm93KGRmMSkqMC4wMSkgI3RvcDElIHJlZ2lvbgogICAgICAgICAgICBkZjEkdG9wMTwtIk4iCiAgICAgICAgICAgIGRmMSR0b3AxWzE6bl08LSJQV1MiCiAgICAgICAgICAgIAogICAgICAgICAgICAjdGIKICAgICAgICAgICAgZGYyPC1jb3YubGlzdFtbcGFzdGUwKCJUQl8iLCB3aW5zaXplKV1dCiAgICAgICAgICAgIGRmMjwtZGYyW29yZGVyKGRmMlssY3ZbaV1dLCBkZWNyZWFzaW5nPVQpLF0KICAgICAgICAgICAgZGYyJHRvcDE8LSJOIgogICAgICAgICAgICBkZjIkdG9wMVsxOm5dPC0iVEIiCiAgICAgICAgICAgIAogICAgICAgICAgICBjbzwtcmJpbmQoZGYxLCBkZjIpCiAgICAKICAgICAgICAgICAgY28kY2hyb208LWZhY3RvcihjbyRjaHJvbSwgbGV2ZWxzPXBhc3RlMCgiY2hyIiwgMToyNikpCiAgICAgICAgICAgIGNvJHRvcDE8LWZhY3RvcihjbyR0b3AxLCBsZXZlbHM9YygiUFdTIiwiVEIiLCJOIikpCiAgICAgICAgICAgIGNvbG5hbWVzKGNvKVt3aGljaChjb2xuYW1lcyhjbyk9PWN2W2ldKV08LSJjb3YiCiAgICAKICAgICAgICAgICAgeW1heDwtbWF4KGNvJGNvdiwgbmEucm09VCkKICAgICAgICAgICAgZ2dwbG90KGNvLCBhZXMoeD1zdGFydC8xMDAwMDAwLCB5PWNvdiwgY29sb3I9dG9wMSkpKwogICAgICAgICAgICAgICAgZ2VvbV9wb2ludChzaXplPTAuNSkrCiAgICAgICAgICAgICAgICBmYWNldF93cmFwKH5jaHJvbSwgbmNvbD00KSsKICAgICAgICAgICAgICAgIHRoZW1lX2NsYXNzaWMoKSt5bGltKC0wLjEseW1heCkrCiAgICAgICAgICAgICAgICBzY2FsZV9jb2xvcl9tYW51YWwodmFsdWVzPWMoImRlZXBwaW5rIiwib3JhbmdlIiAsIiNBREQ4RTY4MCIpLCBsYWJlbHM9YygiUFdTIiwgIlRCIiwgIiIpKSsKICAgICAgICAgICAgICAgIHlsYWIoIkNvdmFyaWFuY2UiKSt4bGFiKCdQb3N0aW9uIChNYiknKSsKICAgICAgICAgICAgICAgIGdndGl0bGUocGFzdGUwKHdpbnNpemVbd10sIiB3aW5kb3cgIixjdltpXSkpKwogICAgICAgICAgICAgICAgc2NhbGVfeF9jb250aW51b3VzKGxhYmVscyA9IGNvbW1hKSsKICAgICAgICAgICAgICAgIGd1aWRlcyhjb2xvciA9IGd1aWRlX2xlZ2VuZChvdmVycmlkZS5hZXMgPSBsaXN0KGNvbG9yPWMoImRlZXBwaW5rIiwib3JhbmdlIiwid2hpdGUiKSksIHRpdGxlPWVsZW1lbnRfdGV4dCgiVG9wIDElIikpKQogICAKICAgICAgICAgICAgICAgIGdnc2F2ZShwYXN0ZTAoIi4uL091dHB1dC9DT1YvQ09Wc2Nhbl8zcG9wLzNQb3BzLiIsY3ZbaV0sIl9wZXJDaHJvbV8iLHdpbnNpemUsICJXaW5kb3dfT3V0bGllcnMucG5nIiksIHdpZHRoID0gMTAsIGhlaWdodCA9IDgsIGRwaT0zMDApCiAgICAgICAKICAgIH0KICAgCiAgICBpZiAoaT09Myl7CiAgICAgICAgCiAgICAgICAgI3B3cwogICAgICAgIGRmMTwtY292Lmxpc3RbW3Bhc3RlMCgiUFdTXyIsIHdpbnNpemVbd10pXV0KICAgICAgICBkZjE8LWRmMVssYygiY2hyb20iLCJzdGFydCIsImVuZCIsImNvdjIzIildCiAgICAgICAgZGYxPC1kZjFbb3JkZXIoZGYxJGNvdjIzLCBkZWNyZWFzaW5nPVQpLF0KICAgICAgICBuPC1jZWlsaW5nKG5yb3coZGYxKSowLjAxKSAjdG9wMSUgcmVnaW9uCiAgICAgICAgZGYxJHRvcDE8LSJOIgogICAgICAgIGRmMSR0b3AxWzE6bl08LSJQV1MiCiAgICAKICAgICAgICAjdGIKICAgICAgICBkZjI8LWNvdi5saXN0W1twYXN0ZTAoIlRCXyIsIHdpbnNpemUpXV0KICAgICAgICBkZjI8LWRmMlssYygiY2hyb20iLCJzdGFydCIsImVuZCIsImNvdjIzIildCiAgICAgICAgZGYyPC1kZjJbb3JkZXIoZGYyJGNvdjIzLCBkZWNyZWFzaW5nPVQpLF0KICAgICAgICBkZjIkdG9wMTwtIk4iCiAgICAgICAgZGYyJHRvcDFbMTpuXTwtIlRCIgogICAgCiAgICAgICAgI3NzCiAgICAgICAgZGYzPC1jb3YubGlzdFtbcGFzdGUwKCJTU18iLCB3aW5zaXplKV0KICAgICAgICBkZjM8LWRmM1ssYygiY2hyb20iLCJzdGFydCIsImVuZCIsImNvdjIzIildCiAgICAgICAgZGYzPC1kZjNbb3JkZXIoZGYzJGNvdjIzLCBkZWNyZWFzaW5nPVQpLF0KICAgICAgICBkZjMkdG9wMTwtIk4iCiAgICAgICAgZGYzJHRvcDFbMTpuXTwtIlNTIgoKICAgICAgICBjbzwtcmJpbmQoZGYxLGRmMixkZjMpCgogICAgICAgIGNvJGNocm9tPC1mYWN0b3IoY28kY2hyb20sIGxldmVscz1wYXN0ZTAoImNociIsIDE6MjYpKQogICAgICAgIGNvJHRvcDE8LWZhY3RvcihjbyR0b3AxLCBsZXZlbHM9YygiUFdTIiwiVEIiLCJTUyIsIk4iKSkKICAgICAgICB5bWF4PC1tYXgoY28kY292MjMsIG5hLnJtPVQpCiAgICAgICAgZ2dwbG90KGNvLCBhZXMoeD1zdGFydC8xMDAwMDAwLCB5PWNvdjIzLCBjb2xvcj10b3AxKSkrCiAgICAgICAgICAgIGdlb21fcG9pbnQoc2l6ZT0wLjUpKwogICAgICAgICAgICBmYWNldF93cmFwKH5jaHJvbSwgbmNvbD00KSsKICAgICAgICAgICAgdGhlbWVfY2xhc3NpYygpK3lsaW0oLTAuMSx5bWF4KSsKICAgICAgICAgICAgeWxhYigiQ292YXJpYW5jZSIpK3hsYWIoJ1Bvc3Rpb24gKE1iKScpKwogICAgICAgICAgICBnZ3RpdGxlKHBhc3RlMCh3aW5zaXplW3ddLCIgd2luZG93ICIsY3ZbaV0pKSsKICAgICAgICAgICAgc2NhbGVfeF9jb250aW51b3VzKGxhYmVscyA9IGNvbW1hKSsKICAgICAgICAgICAgI3NjYWxlX2NvbG9yX2Rpc2NyZXRlKGJyZWFrcz1jKCJQV1MiLCJTUyIsIlRCIikpKwogICAgICAgICAgICBzY2FsZV9jb2xvcl9tYW51YWwodmFsdWVzPWMoImRlZXBwaW5rIiwib3JhbmdlIixncmUsIiNBREQ4RTY2NiIpLCBsYWJlbHM9YygiUFdTIiwiVEIiLCJTUyIsICIiKSkrCiAgICAgICAgICAgIGd1aWRlcyhjb2xvciA9IGd1aWRlX2xlZ2VuZChvdmVycmlkZS5hZXMgPSBsaXN0KGNvbG9yPWMoImRlZXBwaW5rIiwib3JhbmdlIixncmUsICAgICAgICAgIndoaXRlIikpLHRpdGxlPWVsZW1lbnRfdGV4dCgiVG9wIDElIG91dGxpZXJzIikpKSAKICAgICAgICBnZ3NhdmUocGFzdGUwKCIuLi9PdXRwdXQvQ09WL0NPVnNjYW5fM3BvcC8zUG9wcy5jb3YyM18zUG9wc19wZXJDaHJvbV8iLHdpbnNpemUsICJXaW5kb3dfT3V0bGllcnMucG5nIiksIHdpZHRoID0gMTAsIGhlaWdodCA9IDksIGRwaT0zMDApCiAgICAgICAgfQogICAgfQp9CmBgYAoKIVtdKC4uL091dHB1dC9DT1YvQ09Wc2Nhbl8zcG9wLzNQb3BzLmNvdjEyX3BlckNocm9tXzEwMGtXaW5kb3dfT3V0bGllcnMucG5nKQohW10oLi4vT3V0cHV0L0NPVi9DT1ZzY2FuXzNwb3AvM1BvcHMuY292MTNfcGVyQ2hyb21fMTAwa1dpbmRvd19PdXRsaWVycy5wbmcpCgoKIVtdKC4uL091dHB1dC9DT1YvQ09Wc2Nhbl8zcG9wLzNQb3BzLmNvdjIzXzNQb3BzX3BlckNocm9tXzEwMGtXaW5kb3dfT3V0bGllcnMucG5nKSAgCgojIyMgV2hvbGUgZ2Vub21lIHBsb3RzIHdpdGhvdXQgZmFjZXQgd3JhcAoKYGBge3IgZXZhbD1GQUxTRSwgbWVzc2FnZT1GQUxTRSwgd2FybmluZz1GQUxTRX0KY3Y8LWMoImNvdjEyIiwiY292MTMiLCJjb3YyMyIpCmZvciAoaSBpbiAxOmxlbmd0aChjdikpewogICAgaWYgKGk9PTF8aT09Mil7CiAgICAgICAgZGYxPC1jb3YubGlzdFtbIlBXU18xMDBrIl1dCiAgICAgICAgZGYxPC1kZjFbb3JkZXIoZGYxWyxjdltpXV0sIGRlY3JlYXNpbmc9VCksXQogICAgICAgIG48LWNlaWxpbmcobnJvdyhkZjEpKjAuMDEpICN0b3AxJSByZWdpb24KICAgICAgICBkZjEkdG9wMTwtIk4iCiAgICAgICAgZGYxJHRvcDFbMTpuXTwtIlBXUyIKICAgICAgICAKICAgICAgICBkZjI8LWNvdi5saXN0W1siVEJfMTAwayJdXQogICAgICAgIGRmMjwtZGYyW29yZGVyKGRmMlssY3ZbaV1dLCBkZWNyZWFzaW5nPVQpLF0KICAgICAgICBkZjIkdG9wMTwtIk4iCiAgICAgICAgZGYyJHRvcDFbMTpuXTwtIlRCIgoKICAgICAgICBjbzwtcmJpbmQoZGYxLCBkZjIpCiAgICAgICAgY28kY2hyb208LWZhY3RvcihjbyRjaHJvbSwgbGV2ZWxzPXBhc3RlMCgiY2hyIiwgMToyNikpCiAgICAgICAgY29sbmFtZXMoY28pW3doaWNoKGNvbG5hbWVzKGNvKT09Y3ZbaV0pXTwtImNvdiIKICAgICAgICAjYXNzZ2luIGNvbG9ycwogICAgICAgIGNvJHRvcDE8LWFwcGx5KGNvLCAxLCBmdW5jdGlvbih4KSB7aWZlbHNlICh4Wyd0b3AxJ109PSJOIiwgeFsnY29sb3InXSwgeFsndG9wMSddKX0gKQogICAgICAgIGNvJHRvcDE8LWZhY3RvcihjbyR0b3AxLCBsZXZlbHM9YygiUFdTIiwiVEIiLCJjb2wxIiwiY29sMiIpKQogICAgICAgIAogICAgICAgICNjb3VudCB0aGUgbnVtYmVyIG9mIHNpdGVzIHBlciBjaHJvbW9zb21lcwogICAgICAgIHBvc3M8LWRhdGEuZnJhbWUoY2hyPXBhc3RlMCgiY2hyIiwxOjI2KSkKICAgICAgICBrPTEKICAgICAgICBmb3IgKGogaW4gMToyNil7CiAgICAgICAgICAgIGRmPC1kZjFbZGYxJGNocj09cGFzdGUwKCJjaHIiLGopLF0KICAgICAgICAgICAgcG9zcyRzdGFydFtqXTwtawogICAgICAgICAgICBwb3NzJGVuZFtqXTwtaytucm93KGRmKS0xCiAgICAgICAgICAgIGs9aytucm93KGRmKQogICAgICAgIH0KICAgICAgICBwb3NzJHg8LXBvc3Mkc3RhcnQrKHBvc3MkZW5kLXBvc3Mkc3RhcnQpLzIKICAgICAgICB5bWF4PC1tYXgoY28kY292LCBuYS5ybT1UKQogICAgICAgIGdncGxvdChjbywgYWVzKHg9aW5kZXgsIHk9Y292LCBjb2xvcj10b3AxKSkrCiAgICAgICAgICAgIGdlb21fcG9pbnQoc2l6ZT0wLjUpKwogICAgICAgICAgICB0aGVtZV9jbGFzc2ljKCkreWxpbSgtMC4xLHltYXgpKwogICAgICAgICAgICBzY2FsZV9jb2xvcl9tYW51YWwodmFsdWVzPWMoIiNGRjMyOTNCMyIsIiMwNDMzRkZCMyIgLCIjQUREOEU2ODAiLCIjQzBDMEMwODgiKSwgbGFiZWxzPWMoIlBXUyIsICJUQiIsICIiLCIiKSkrCiAgICAgICAgICAgICAgICB5bGFiKCJDb3ZhcmlhbmNlIikrCiAgICAgICAgICAgICAgICBnZ3RpdGxlKHBhc3RlMCgiIDEwMGsgd2luZG93ICIsY3ZbaV0pKSsKICAgICAgICAgICAgICAgIGd1aWRlcyhjb2xvciA9IGd1aWRlX2xlZ2VuZChvdmVycmlkZS5hZXMgPSBsaXN0KGNvbG9yPWMoImRlZXBwaW5rIiwiIzA0MzNGRiIsIndoaXRlIiwid2hpdGUiKSwgc2l6ZT0yKSwgdGl0bGU9ZWxlbWVudF90ZXh0KCJPdXRsaWVyICgxJSkiLCBzaXplPTEwKSkpKwogICAgICAgICAgICBzY2FsZV94X2NvbnRpbnVvdXMobmFtZT0iQ2hyb21vc29tZSIsIGJyZWFrcz1wb3NzJHgsIGxhYmVscz0xOjI2KQogICAgICAgIGdnc2F2ZShwYXN0ZTAoIi4uL091dHB1dC9DT1YvQ09Wc2Nhbl8zcG9wLzNQb3BzLiIsY3ZbaV0sIl8xMDBrX1dpbmRvd19PdXRsaWVycy5wbmciKSwgd2lkdGggPSAxMCwgaGVpZ2h0ID0gMy41LCBkcGk9MzAwKQogICAgICAgIH19CiAgICBpZiAoaT09Myl7CiAgICAgICAgI3B3cwogICAgICAgIGRmMTwtY292Lmxpc3RbWyJQV1NfMTAwayJdXQogICAgICAgIGRmMTwtZGYxWyxjKCJjaHJvbSIsInN0YXJ0IiwiZW5kIiwiY292MjMiLCJpbmRleCIsImNvbG9yIildCiAgICAgICAgZGYxPC1kZjFbb3JkZXIoZGYxJGNvdjIzLCBkZWNyZWFzaW5nPVQpLF0KICAgICAgICBuPC1jZWlsaW5nKG5yb3coZGYxKSowLjAxKSAjdG9wMSUgcmVnaW9uCiAgICAgICAgZGYxJHRvcDE8LWRmMSRjb2xvcgogICAgICAgIGRmMSR0b3AxWzE6bl08LSJQV1MiCiAgICAgICAgI3RiCiAgICAgICAgZGYyPC1jb3YubGlzdFtbIlRCXzEwMGsiXV0KICAgICAgICBkZjI8LWRmMlssYygiY2hyb20iLCJzdGFydCIsImVuZCIsImNvdjIzIiwiaW5kZXgiLCJjb2xvciIpXQogICAgICAgIGRmMjwtZGYyW29yZGVyKGRmMiRjb3YyMywgZGVjcmVhc2luZz1UKSxdCiAgICAgICAgZGYyJHRvcDE8LWRmMiRjb2xvcgogICAgICAgIGRmMiR0b3AxWzE6bl08LSJUQiIKICAgICAgICAjc3MKICAgICAgICBkZjM8LWNvdi5saXN0W1siU1NfMTAwayJdXQogICAgICAgIGRmMzwtZGYzWyxjKCJjaHJvbSIsInN0YXJ0IiwiZW5kIiwiY292MjMiLCJpbmRleCIsImNvbG9yIildCiAgICAgICAgZGYzPC1kZjNbb3JkZXIoZGYzJGNvdjIzLCBkZWNyZWFzaW5nPVQpLF0KICAgICAgICBkZjMkdG9wMTwtZGYzJGNvbG9yCiAgICAgICAgZGYzJHRvcDFbMTpuXTwtIlNTIgoKICAgICAgICBjbzwtcmJpbmQoZGYxLGRmMixkZjMpCiAgICAgICAgY28kdG9wMTwtZmFjdG9yKGNvJHRvcDEsIGxldmVscz1jKCJQV1MiLCJUQiIsIlNTIiwiY29sMSIsImNvbDIiKSkKICAgICAgICAKICAgICAgICAjY291bnQgdGhlIG51bWJlciBvZiBzaXRlcyBwZXIgY2hyb21vc29tZXMKICAgICAgICBwb3NzPC1kYXRhLmZyYW1lKGNocj1wYXN0ZTAoImNociIsMToyNikpCiAgICAgICAgaz0xCiAgICAgICAgZm9yIChqIGluIDE6MjYpewogICAgICAgICAgICBkZjwtZGYxW2RmMSRjaHI9PXBhc3RlMCgiY2hyIixqKSxdCiAgICAgICAgICAgIHBvc3Mkc3RhcnRbal08LWsKICAgICAgICAgICAgcG9zcyRlbmRbal08LWsrbnJvdyhkZiktMQogICAgICAgICAgICBrPWsrbnJvdyhkZikKICAgICAgICB9CiAgICAgICAgcG9zcyR4PC1wb3NzJHN0YXJ0Kyhwb3NzJGVuZC1wb3NzJHN0YXJ0KS8yCiAgICAgICAgeW1heDwtbWF4KGNvJGNvdiwgbmEucm09VCkKICAgICAgICBnZ3Bsb3QoY28sIGFlcyh4PWluZGV4LCB5PWNvdjIzLCBjb2xvcj10b3AxKSkrCiAgICAgICAgICAgIGdlb21fcG9pbnQoc2l6ZT0wLjUpKwogICAgICAgICAgICB0aGVtZV9jbGFzc2ljKCkreWxpbSgtMC4xLHltYXgpKwogICAgICAgICAgICBzY2FsZV9jb2xvcl9tYW51YWwodmFsdWVzPWMoIiNGRjMyOTNCMyIsIiMwNDMzRkZCMyIsIiMwMDhGMDBCMyIgLCIjQUREOEU2ODAiLCIjQzBDMEMwODgiKSwgbGFiZWxzPWMoIlBXUyIsICJUQiIsIlNTIiwgIiIsIiIpKSsKICAgICAgICAgICAgICAgIHlsYWIoIkNvdmFyaWFuY2UiKSsKICAgICAgICAgICAgICAgIGdndGl0bGUocGFzdGUwKCIgMTAwayB3aW5kb3cgIixjdltpXSkpKwogICAgICAgICAgICAgICAgZ3VpZGVzKGNvbG9yID0gZ3VpZGVfbGVnZW5kKG92ZXJyaWRlLmFlcyA9IGxpc3QoY29sb3I9YygiZGVlcHBpbmsiLCIjMDQzM0ZGIiwiIzAwOEYwMCIsIndoaXRlIiwid2hpdGUiKSwgc2l6ZT0yKSwgdGl0bGU9ZWxlbWVudF90ZXh0KCJPdXRsaWVyICgxJSkiKSkpKwogICAgICAgICAgICBzY2FsZV94X2NvbnRpbnVvdXMobmFtZT0iQ2hyb21vc29tZSIsIGJyZWFrcz1wb3NzJHgsIGxhYmVscz0xOjI2KSsKICAgICAgICAgICAgdGhlbWUobGVnZW5kLnRpdGxlID0gZWxlbWVudF90ZXh0KHNpemU9MTApKQogICAgICAgIGdnc2F2ZShwYXN0ZTAoIi4uL091dHB1dC9DT1YvQ09Wc2Nhbl8zcG9wLzNQb3BzLiIsY3ZbaV0sIl8xMDBrX1dpbmRvd19PdXRsaWVycy5wbmciKSwgd2lkdGggPSAxMCwgaGVpZ2h0ID0gMy41LCBkcGk9MzAwKQogICAgICAgIH0KICAgIH0KfQoKIyMgUGxvdCAzIHRpbWUgcGVyaW9kcyB0b2dldGhlciBmb3IgUFdTIGFuZCBUQgpDb3Y8LWRhdGEuZnJhbWUoKQpmb3IgKGkgaW4gMTpsZW5ndGgoY3YpKXsKICAgIGRmMTwtY292Lmxpc3RbWyJQV1NfMTAwayJdXQogICAgZGYxPC1kZjFbb3JkZXIoZGYxWyxjdltpXV0sIGRlY3JlYXNpbmc9VCksXQogICAgbjwtY2VpbGluZyhucm93KGRmMSkqMC4wMSkgI3RvcDElIHJlZ2lvbgogICAgZGYxJHRvcDE8LSJOIgogICAgZGYxJHRvcDFbMTpuXTwtIlBXUyIKICAgIAogICAgZGYyPC1jb3YubGlzdFtbIlRCXzEwMGsiXV0KICAgIGRmMjwtZGYyW29yZGVyKGRmMlssY3ZbaV1dLCBkZWNyZWFzaW5nPVQpLF0KICAgIGRmMiR0b3AxPC0iTiIKICAgIGRmMiR0b3AxWzE6bl08LSJUQiIKICAgIGNvPC1yYmluZChkZjEsIGRmMikKICAgIGNvJGNocm9tPC1mYWN0b3IoY28kY2hyb20sIGxldmVscz1wYXN0ZTAoImNociIsIDE6MjYpKQogICAgY29sbmFtZXMoY28pW3doaWNoKGNvbG5hbWVzKGNvKT09Y3ZbaV0pXTwtImNvdiIKICAgICNhc3NnaW4gY29sb3JzCiAgICBjbyR0b3AxPC1hcHBseShjbywgMSwgZnVuY3Rpb24oeCkge2lmZWxzZSAoeFsndG9wMSddPT0iTiIsIHhbJ2NvbG9yJ10sIHhbJ3RvcDEnXSl9ICkKICAgIGNvJHRvcDE8LWZhY3RvcihjbyR0b3AxLCBsZXZlbHM9YygiUFdTIiwiVEIiLCJjb2wxIiwiY29sMiIpKQogICAgY28kdGltZTwtY3ZbaV0KICAgIAogICAgQ292PC1yYmluZChDb3YsIGNvWyxjKCJpbmRleCIsICJjb3YiLCJ0b3AxIiwidGltZSIpXSkKfQoKI2NvdW50IHRoZSBudW1iZXIgb2Ygc2l0ZXMgcGVyIGNocm9tb3NvbWVzCmRmMTwtY292Lmxpc3RbWyJQV1NfMTAwayJdXQpwb3NzPC1kYXRhLmZyYW1lKGNocj1wYXN0ZTAoImNociIsMToyNikpCms9MQpmb3IgKGogaW4gMToyNil7CiAgICAgICAgZGY8LWRmMVtkZjEkY2hyPT1wYXN0ZTAoImNociIsaiksXQogICAgICAgIHBvc3Mkc3RhcnRbal08LWsKICAgICAgICBwb3NzJGVuZFtqXTwtaytucm93KGRmKS0xCiAgICAgICAgaz1rK25yb3coZGYpCn0KcG9zcyR4PC1wb3NzJHN0YXJ0Kyhwb3NzJGVuZC1wb3NzJHN0YXJ0KS8yCnltYXg8LW1heChjbyRjb3YsIG5hLnJtPVQpCmdncGxvdChDb3YsIGFlcyh4PWluZGV4LCB5PWNvdiwgY29sb3I9dG9wMSkpKwogICAgZmFjZXRfd3JhcCh+dGltZSwgbmNvbD0xKSsKICAgIGdlb21fcG9pbnQoc2l6ZT0wLjUpKwogICAgdGhlbWVfY2xhc3NpYygpK3lsaW0oLTAuMSx5bWF4KSsKICAgIHNjYWxlX2NvbG9yX21hbnVhbCh2YWx1ZXM9YygiI0ZGMzI5M0IzIiwiIzA0MzNGRkIzIiAsIiNBREQ4RTY4MCIsIiNDMEMwQzA4OCIpLCBsYWJlbHM9YygiUFdTIiwgIlRCIiwgIiIsIiIpKSsKICAgIHlsYWIoIkNvdmFyaWFuY2UiKSsKICAgIGd1aWRlcyhjb2xvciA9IGd1aWRlX2xlZ2VuZChvdmVycmlkZS5hZXMgPSBsaXN0KGNvbG9yPWMoImRlZXBwaW5rIiwiIzA0MzNGRiIsIndoaXRlIiwid2hpdGUiKSwgc2l6ZT0yKSwgdGl0bGU9ZWxlbWVudF90ZXh0KCJPdXRsaWVyICgxJSkiLCBzaXplPTEwKSkpKwogICAgc2NhbGVfeF9jb250aW51b3VzKG5hbWU9IkNocm9tb3NvbWUiLCBicmVha3M9cG9zcyR4LCBsYWJlbHM9MToyNikKCmdnc2F2ZShwYXN0ZTAoIi4uL091dHB1dC9DT1YvQ09Wc2Nhbl8zcG9wL1BXU19UQl8xMDBrX1dpbmRvd19PdXRsaWVycy5wbmciKSwgd2lkdGggPSAxMSwgaGVpZ2h0ID0gNSwgZHBpPTMwMCkKICAgICAgICB9fQoKCmBgYAoKIVtdKC4uL091dHB1dC9DT1YvQ09Wc2Nhbl8zcG9wL1BXU19UQl8xMDBrX1dpbmRvd19PdXRsaWVycy5wbmcpCgoKCiMjIE92ZXJsYXBwaW5nIG91dGxpZXIgcmVnaW9ucyBiZXR3ZWVuIGRpZmZlcmVudCBwb3B1bGF0aW9ucyAKYGBge3IgZXZhbD1GQUxTRSwgbWVzc2FnZT1GQUxTRSwgd2FybmluZz1GQUxTRX0KIzEwMGsKY3Y8LWMoImNvdjEyIiwiY292MTMiLCJjb3YyMyIpCnBhaXJzPC10KGNvbWJuKHBvcHMsIDIpKQpwYWlyczwtZGF0YS5mcmFtZShwYWlycykKY29sbmFtZXMocGFpcnMpPC1wYXN0ZTAoInBvcCIsMToyKQpPdl9kaXJlY3Q8LWRhdGEuZnJhbWUoY292PWMoY3ZbMToyXSwiY292MjMtUFQiLCJjb3YyMy1QUyIsImNvdjIzLVNUIiAsImNvdjIzLTMiKSkKT3ZfMzAwPC1kYXRhLmZyYW1lKGNvdj1jKGN2WzE6Ml0sImNvdjIzLVBUIiwiY292MjMtUFMiLCJjb3YyMy1TVCIgLCJjb3YyMy0zIikpCmZvciAoaSBpbiAxOmxlbmd0aChjdikpewogICAgZGY8LXJlYWQuY3N2KHBhc3RlMCgiLi4vT3V0cHV0L0NPVi9DT1ZzY2FuXzNwb3AvM3BvcHNfdG9wMXBlcmNlbnRfb3V0bGllcl9yZWdpb25zLiIsIGN2W2ldLCAiX25ldy5jc3YiKSkKICAgIGRmPC1kZltkZiR3aW5kb3c9PSIxMDBrIixdCiAgICBkZiRpZDwtcGFzdGUwKGRmJGNocm9tLCJfIixkZiRzdGFydCkKICAgIAogICAgaWYgKGkhPTMpewogICAgICAgICNleGFjdCBvdmVybGFwcwogICAgICAgIGlzZWM8LWludGVyc2VjdChkZiRpZFtkZiRwb3A9PSJQV1MiXSwgZGYkaWRbZGYkcG9wPT0iVEIiXSkgCiAgICAgICAgT3ZfZGlyZWN0JGNvdW50W2ldPC1sZW5ndGgoaXNlYykKICAgICAgICAKICAgICAgICAjIyMjIENoZWNrIGNocm9tb3NvbWUgcmVnaW9uIG92ZXJsYXAgKy0yMDAsMDAwIGJhc2VzCiAgICAgICAgcG9wMTwtZGZbZGYkcG9wPT0iUFdTIixdCiAgICAgICAgcG9wMjwtZGZbZGYkcG9wPT0iVEIiLF0KICAgICAgICBvdmVybHBzPC1kYXRhLmZyYW1lKCkKICAgICAgICBvdmVybHBzMjwtZGF0YS5mcmFtZSgpCiAgICAgICAgZm9yIChuIGluIDE6IG5yb3cocG9wMSkpewogICAgICAgICAgICByZTwtcG9wMltwb3AyJGNocm9tPT1wb3AxJGNocm9tW25dLF0KICAgICAgICAgICAgaWYgKG5yb3cocmUpPj0xKXsKICAgICAgICAgICAgICAgIGZvciAocyBpbiAxOiBucm93KHJlKSl7CiAgICAgICAgICAgICAgICAgICAgaWYgKHJlJHN0YXJ0W3NdPD1wb3AxJHN0YXJ0W25dKzMwMDAwMCAmIHJlJHN0YXJ0W3NdPj1wb3AxJHN0YXJ0W25dLTMwMDAwMCl7CiAgICAgICAgICAgICAgICAgICAgICAgIG92ZXJscHM8LXJiaW5kKG92ZXJscHMsIHJlW3MsXSkKICAgICAgICAgICAgICAgICAgICAgICAgb3ZlcmxwczI8LXJiaW5kKG92ZXJscHMyLHBvcDFbbixdKX0KICAgICAgICAgICAgICAgIH0KICAgICAgICAgICAgfQogICAgICAgIH0KICAgICAgICAjIE1lcmdlIHR3byB0YWJsZXMgaW50byBvbmUgc3VtbWFyeSBvdmVybGFwIHRhYmxlOgogICAgICAgIG92PC1kYXRhLmZyYW1lKGlkPW92ZXJscHMkaWQpCiAgICAgICAgZm9yIChuIGluIDE6IG5yb3cob3ZlcmxwcykpewogICAgICAgICAgICBpZiAob3ZlcmxwcyRzdGFydFtuXTxvdmVybHBzMiRzdGFydFtuXSkge292JHN0YXJ0W25dPC1vdmVybHBzJHN0YXJ0W25dOyBvdiRlbmRbbl08LW92ZXJscHMyJGVuZFtuXX0KICAgICAgICAgICAgaWYgKG92ZXJscHMkc3RhcnRbbl0+PW92ZXJscHMyJHN0YXJ0W25dKSB7b3Ykc3RhcnRbbl08LW92ZXJscHMyJHN0YXJ0W25dO292JGVuZFtuXTwtb3ZlcmxwcyRlbmRbbl19CiAgICAgICAgfQogICAgICAgIG92WywiY292LlBXUyJdPC1vdmVybHBzWyw0XQogICAgICAgIG92WywiY292LlRCIl08LW92ZXJscHMyWyw0XQogICAgICAgIHdyaXRlLmNzdihvdiwgcGFzdGUwKCIuLi9PdXRwdXQvQ09WL0NPVnNjYW5fM3BvcC9PdmVybGFwX3JlZ2lvbnNfIixjdltpXSwiX3BsdXNtaW51czMwMGsuY3N2IiksIHJvdy5uYW1lcyA9IEYpCiAgICAgICAgT3ZfMzAwJGNvdW50W2ldPC1ucm93KG92KQogICAgICAgIH0KICAgICAgICAKICAgIGlmIChpPT0zKXsKICAgICAgICBpc2VjPC1pbnRlcnNlY3QoZGYkaWRbZGYkcG9wPT0iUFdTIl0sIGRmJGlkW2RmJHBvcD09IlRCIl0pIAogICAgICAgIGlzZWMyPC1pbnRlcnNlY3QoZGYkaWRbZGYkcG9wPT0iUFdTIl0sIGRmJGlkW2RmJHBvcD09IlNTIl0pIAogICAgICAgIGlzZWMzPC1pbnRlcnNlY3QoZGYkaWRbZGYkcG9wPT0iU1MiXSwgZGYkaWRbZGYkcG9wPT0iVEIiXSkgCiAgICAgICAgT3ZfZGlyZWN0JGNvdW50W2ldPC1sZW5ndGgoaXNlYykKICAgICAgICBPdl9kaXJlY3QkY291bnRbaSsxXTwtbGVuZ3RoKGlzZWMyKQogICAgICAgIE92X2RpcmVjdCRjb3VudFtpKzJdPC1sZW5ndGgoaXNlYykKICAgICAgICBPdl9kaXJlY3QkY291bnRbaSszXTwtbGVuZ3RoKGludGVyc2VjdChkZiRpZFtkZiRwb3A9PSJTUyJdLCBpbnRlcnNlY3QoZGYkaWRbZGYkcG9wPT0iUFdTIl0sIGRmJGlkW2RmJHBvcD09IlRCIl0pKSkKICAgICAgICAKICAgICAgICBmb3IoaiBpbiAxOm5yb3cocGFpcnMpKXsKICAgICAgICAjIyMjIENoZWNrIGNocm9tb3NvbWUgcmVnaW9uIG92ZXJsYXAgKy0yMDAsMDAwIGJhc2VzCiAgICAgICAgICAgIHBvcDE8LWRmW2RmJHBvcD09cGFpcnNbaiwxXSxdCiAgICAgICAgICAgIHBvcDI8LWRmW2RmJHBvcD09cGFpcnNbaiwyXSxdCiAgICAgICAgICAgIG92ZXJscHM8LWRhdGEuZnJhbWUoKQogICAgICAgICAgICBvdmVybHBzMjwtZGF0YS5mcmFtZSgpCiAgICAgICAgICAgIGZvciAobiBpbiAxOiBucm93KHBvcDEpKXsKICAgICAgICAgICAgICAgIHJlPC1wb3AyW3BvcDIkY2hyb209PXBvcDEkY2hyb21bbl0sXQogICAgICAgICAgICAgICAgaWYgKG5yb3cocmUpPj0xKXsKICAgICAgICAgICAgICAgICAgICBmb3IgKHMgaW4gMTogbnJvdyhyZSkpewogICAgICAgICAgICAgICAgICAgICAgICBpZiAocmUkc3RhcnRbc108PXBvcDEkc3RhcnRbbl0rMzAwMDAwICYgcmUkc3RhcnRbc10+PXBvcDEkc3RhcnRbbl0tMzAwMDAwKXsKICAgICAgICAgICAgICAgICAgICAgICAgICAgIG92ZXJscHM8LXJiaW5kKG92ZXJscHMsIHJlW3MsXSkKICAgICAgICAgICAgICAgICAgICAgICAgICAgIG92ZXJscHMyPC1yYmluZChvdmVybHBzMixwb3AxW24sXSl9CiAgICAgICAgICAgICAgICAgICAgfQogICAgICAgICAgICAgICAgfQogICAgICAgICAgICB9CiAgICAgICAgIyBNZXJnZSB0d28gdGFibGVzIGludG8gb25lIHN1bW1hcnkgb3ZlcmxhcCB0YWJsZToKICAgICAgICAgICAgb3Y8LWRhdGEuZnJhbWUoaWQ9b3ZlcmxwcyRpZCkKICAgICAgICAgICAgZm9yIChuIGluIDE6IG5yb3cob3ZlcmxwcykpewogICAgICAgICAgICAgICAgaWYgKG92ZXJscHMkc3RhcnRbbl08b3ZlcmxwczIkc3RhcnRbbl0pIHtvdiRzdGFydFtuXTwtb3ZlcmxwcyRzdGFydFtuXTsgb3YkZW5kW25dPC1vdmVybHBzMiRlbmRbbl19CiAgICAgICAgICAgICAgICBpZiAob3ZlcmxwcyRzdGFydFtuXT49b3ZlcmxwczIkc3RhcnRbbl0pIHtvdiRzdGFydFtuXTwtb3ZlcmxwczIkc3RhcnRbbl07b3YkZW5kW25dPC1vdmVybHBzJGVuZFtuXX0KICAgICAgICAgICAgfQogICAgICAgIAogICAgICAgICAgICBvdlsscGFzdGUwKCJjb3YuIixwYWlyc1tqLDFdKV08LW92ZXJscHNbLDRdCiAgICAgICAgICAgIG92WyxwYXN0ZTAoImNvdi4iLHBhaXJzW2osMl0pXTwtb3ZlcmxwczJbLDRdCiAgICAgICAgICAgIG92PC1vdlshZHVwbGljYXRlZChvdiksXQogICAgICAgICAgICB3cml0ZS5jc3Yob3YsIHBhc3RlMCgiLi4vT3V0cHV0L0NPVi9DT1ZzY2FuXzNwb3AvT3ZlcmxhcF9yZWdpb25zXyIsY3ZbaV0sIl8iLHBhaXJzW2osMV0sIi4iLCBwYWlyc1tqLDJdLCJfcGx1c21pbnVzMzAway5jc3YiKSwgcm93Lm5hbWVzID0gRikKICAgICAgICAgICAgT3ZfMzAwJGNvdW50W2krai0xXTwtbnJvdyhvdikKICAgIH0KICAgIH0KfQp3cml0ZS5jc3YoT3ZfZGlyZWN0LCBwYXN0ZTAoIi4uL091dHB1dC9DT1YvQ09Wc2Nhbl8zcG9wL0REaXJlY3RfT3ZlcmxhcHBpbmdfcmVnaW9uc19jb3VudHNfM3BvcF9zdW1tYXJ5LmNzdiIpKQpPdl8zMDAkY291bnRbNl08LU5BCndyaXRlLmNzdihPdl8zMDAsIHBhc3RlMCgiLi4vT3V0cHV0L0NPVi9DT1ZzY2FuXzNwb3AvT3ZlcmxhcHBpbmdfcmVnaW9uc19jb3VudHNfM3BvcF9wbHVzTWludXMzMDBrLmNzdiIpKQoKYGBgCgoKIyBSdW4gdGhlIHNucEVmZiBwaXBlbGluZSB0byBmaW5kIGFubm90YXRpb24gaW4gdGhlIG91dGxpZXIgcmVnaW9ucyAoMTAway13aW5kb3crLTEwMGspICAKCiMjIENyZWF0ZSBhIHNjcmlwdCB0byBydW4gU25wRWZmIAoKQ3JlYXRlIFZDRiBmaWxlcyB3aXRoIHNlbGVjdGVkIHJlZ2lvbnMgJiBydW4gc25wRWZmICAKYGBge3IgZXZhbD1GQUxTRSwgbWVzc2FnZT1GQUxTRSwgd2FybmluZz1GQUxTRX0KI0NyZWF0ZSBiZWQgZmlsZXMKY3Y8LWMoImNvdjEyIiwiY292MTMiLCJjb3YyMyIpCiNQcmV2ZW50IHNjaWVudGlmaWMgbm90YXRpb24gaW4gYmVkIGZpbGVzCm9wdGlvbnMoc2NpcGVuPTk5OSkKCiNUaGUgZmlyc3QgbGluZSBvZiBiZWQgZmlsZXMgaXMgb2Z0ZW4gbm90IHJlZCBieSB2Y2Z0b29scwpmb3IgKGkgaW4gMTozKXsKICAgIGRmPC1yZWFkLmNzdihwYXN0ZTAoIi4uL091dHB1dC9DT1YvQ09Wc2Nhbl8zcG9wLzNwb3BzX3RvcDFwZXJjZW50X291dGxpZXJfcmVnaW9ucy4iLGN2W2ldLCJfbmV3LmNzdiIpKQogICAgI3NlbGVjdCAxMDBrCiAgICBkZjwtZGZbZGYkd2luZG93PT0iMTAwayIsXQogICAgI2FkZCAxMDBrCiAgICBkZiRzdGFydDwtZGYkc3RhcnQtMTAwMDAwCiAgICBkZiRlbmQ8LWRmJGVuZCsxMDAwMDAKICAgIGRmcDwtZGZbZGYkcG9wPT0iUFdTIiwxOjNdCiAgICBjb2xuYW1lcyhkZnApPC1jKCd0cmFjayB0eXBlPWJlZEdyYXBoJywgJzEnLCcxJykKICAgIHdyaXRlLnRhYmxlKGRmcCwgcGFzdGUwKCIuLi9PdXRwdXQvQ09WL0NPVnNjYW5fM3BvcC9QV1Nfb3V0bGllcnNfIixjdltpXSwiX25ldy5iZWQiKSxxdW90ZSA9IEYsIHJvdy5uYW1lcyA9IEYsIGNvbC5uYW1lcyA9IFQsc2VwID0gIlx0IikKICAgIGRmdDwtZGZbZGYkcG9wPT0iVEIiLDE6M10KICAgIGNvbG5hbWVzKGRmdCk8LWMoJ3RyYWNrIHR5cGU9YmVkR3JhcGgnLCAnMScsJzEnKQogICAgd3JpdGUudGFibGUoZGZ0LCBwYXN0ZTAoIi4uL091dHB1dC9DT1YvQ09Wc2Nhbl8zcG9wL1RCX291dGxpZXJzXyIsY3ZbaV0sIl9uZXcuYmVkIikscXVvdGUgPSBGLCByb3cubmFtZXMgPSBGLCBjb2wubmFtZXMgPSBGLHNlcCA9ICJcdCIpCiAgICAKICAgIGlmIChpPT0zKXsKICAgICAgICBkZnM8LWRmW2RmJHBvcD09IlNTIiwxOjNdCiAgICAgICAgY29sbmFtZXMoZGZzKTwtYygndHJhY2sgdHlwZT1iZWRHcmFwaCcsICcxJywnMScpCiAgICAgICAgd3JpdGUudGFibGUoZGZzLCBwYXN0ZTAoIi4uL091dHB1dC9DT1YvQ09Wc2Nhbl8zcG9wL1NTX291dGxpZXJzXyIsY3ZbaV0sIl9uZXcuYmVkIikscXVvdGUgPSBGLCByb3cubmFtZXMgPSBGLCBjb2wubmFtZXMgPSBGLHNlcCA9ICJcdCIpCiAgICB9Cn0KCiMgQ3JlYXRlIGEgYmFzaCBzY3JpcHQgdG8gcnVuIHNucEVmZgpiZWRmaWxlczwtbGlzdC5maWxlcygiLi4vT3V0cHV0L0NPVi9DT1ZzY2FuXzNwb3AvIiwgcGF0dGVybj0iKl9uZXcuYmVkIikKCnNpbmsoIi4uL0NPVnNjYW5fY3JlYXRlVkNGc18zUG9wczIuc2giKQpjYXQoIiMhL2Jpbi9iYXNoIFxuXG4iKQpmb3IgKGkgaW4gMTpsZW5ndGgoYmVkZmlsZXMpKXsKICAgIGZuYW1lPC1nc3ViKCIuYmVkIiwnJywgYmVkZmlsZXNbaV0pCiAgICBjYXQocGFzdGUwKCJ2Y2Z0b29scyAtLWd6dmNmIERhdGEvbmV3X3ZjZi8zcG9wLzNwb3BzLk1ENzAwMF9OUzAuNV9tYWYwNS52Y2YuZ3ogLS1iZWQgT3V0cHV0L0NPVi9DT1ZzY2FuXzNwb3AvIiwgYmVkZmlsZXNbaV0sICIgLS1vdXQgT3V0cHV0L0NPVi9DT1ZzY2FuXzNwb3AvIiwgZm5hbWUsIiAtLXJlY29kZSAtLWtlZXAtSU5GTy1hbGwgXG4iKSkKfQpzaW5rKE5VTEwpICAKCiNjcmVhdGUgYSBiYXNoIHNjcmlwdCB0byBydW4gc25wRWZmCnZmaWxlczwtbGlzdC5maWxlcygiLi4vT3V0cHV0L0NPVi9DT1ZzY2FuXzNwb3AvIiwgcGF0dGVybj0iLnJlY29kZS52Y2YiKQoKc2luaygifi9wcm9ncmFtcy9zbnBFZmYvcnVuc25wRWZmX2Nvdl8zcG9wLnNoIikKY2F0KCIjIS9iaW4vYmFzaCBcblxuIikKZm9yIChpIGluIDE6bGVuZ3RoKHZmaWxlcykpewogICAgZm5hbWU8LWdzdWIoIl9uZXcucmVjb2RlLnZjZiIsIiIsdmZpbGVzW2ldKQogICAgY2F0KHBhc3RlMCgiamF2YSAtWG14OGcgLWphciBzbnBFZmYuamFyIENoX3YyLjAuMi45OSB+L1Byb2plY3RzL1BhY0hlcnJpbmcvT3V0cHV0L0NPVi9DT1ZzY2FuXzNwb3AvIix2ZmlsZXNbaV0sICIgLXN0YXRzIH4vUHJvamVjdHMvUGFjSGVycmluZy9PdXRwdXQvQ09WL0NPVnNjYW5fM3BvcC8iLGZuYW1lLCIuaHRtbCA+ICB+L1Byb2plY3RzL1BhY0hlcnJpbmcvT3V0cHV0L0NPVi9DT1ZzY2FuXzNwb3AvQW5uby4iLGZuYW1lLCIudmNmIFxuIikpCiAgICAKICAgICNleHRyYWN0IHRoZSBhbm5vdGF0aW9uIGluZm9ybWF0aW9uCiAgICBjYXQocGFzdGUwKCJiY2Z0b29scyBxdWVyeSAtZiAnJUNIUk9NICVQT1MgJUlORk8vQUYgJUlORk8vQU5OXFxuJyB+L1Byb2plY3RzL1BhY0hlcnJpbmcvT3V0cHV0L0NPVi9DT1ZzY2FuXzNwb3AvQW5uby4iLGZuYW1lLCIudmNmID4gfi9Qcm9qZWN0cy9QYWNIZXJyaW5nL091dHB1dC9DT1YvQ09Wc2Nhbl8zcG9wLyIsZm5hbWUsIl9hbm5vdGF0aW9uIFxuXG4iKSkKCn0Kc2luayhOVUxMKSAgCgoKYGBgCgpgYGB7YmFzaCBldmFsPUZBTFNFLCBpbmNsdWRlPUZBTFNFfQpjZCB+L1Byb2plY3RzL1BhY0hlcnJpbmcKYmFzaCBDT1ZzY2FuX2NyZWF0ZVZDRnNfM3BvcHMuc2gKCmNkIH4vcHJvZ3JhbXMvc25wRWZmCmJhc2ggcnVuc25wRWZmX2Nvdl8zcG9wLnNoCmBgYAoKCgojIyBDcmVhdGUgc3VtbWFyeSBnZW5lIGZpbGVzIGZyb20gc25wRWZmIGFuZCBjaGVjayBvdmVybGFwcGluZyBnZW5lcy4KCmBgYHtyIGV2YWw9RkFMU0UsIG1lc3NhZ2U9RkFMU0UsIHdhcm5pbmc9RkFMU0V9CiMjIENyZWF0ZSBzdW1tYXJ5IGZpbGVzIG9mIHNucEVmZiByZXN1bHRzIChnZW5lIGFubm90YXRpb25zIGluIHRoZSByZWdpb25zIG9mIGludGVyZXN0KSBhbmQgcmVmb3JtYXQgYXMgYSBTaGlueUdvIGlucHV0IAoKI2NyZWF0ZSBnZW5lIGxpc3QgCmdmaWxlczwtbGlzdC5maWxlcygiLi4vT3V0cHV0L0NPVi9DT1ZzY2FuXzNwb3AvIiwgcGF0dGVybj0iZ2VuZXMudHh0IikKCmZvciAoaSBpbiAxOmxlbmd0aChnZmlsZXMpKXsKICAgIGRmPC1yZWFkLnRhYmxlKHBhc3RlMCgiLi4vT3V0cHV0L0NPVi9DT1ZzY2FuXzNwb3AvIixnZmlsZXNbaV0pLCBzZXA9Ilx0IikKICAgIGRmPC1kZlssMTo3XQogICAgY29sbmFtZXMoZGYpPC1jKCJHZW5lTmFtZSIsIkdlbmVJZCIsIlRyYW5zY3JpcHRJZCIsIkJpb1R5cGUiLCJ2YXJpYW50c19pbXBhY3RfSElHSCIsInZhcmlhbnRzX2ltcGFjdF9MT1ciLAkidmFyaWFudHNfaW1wYWN0X01PREVSQVRFIikKICAgIAogICAgZm5hbWU8LWdzdWIoIi5nZW5lcy50eHQiLCIiLGdmaWxlc1tpXSkKICAgIGdlbmVzPC11bmlxdWUoZGYkR2VuZUlkKQogICAgc2luayhwYXN0ZTAoIi4uL091dHB1dC9DT1YvQ09Wc2Nhbl8zcG9wL2dlbmVJRGxpc3RfIixmbmFtZSwiLnR4dCIpKQogICAgY2F0KHBhc3RlMChnZW5lcywiOyAiKSkKICAgIHNpbmsoTlVMTCkKfQoKI0Fubm90YXRpb24gaW5mb3IgZnJvbSBTbnBFZmYKY3Y8LWMoImNvdjEyIiwiY292MTMiLCJjb3YyMyIpCmZvciAoYyBpbiAxOjMpewogICAgaWYgKGMhPTMpewogICAgZm9yIChwIGluIDE6Mil7CiAgICAgICAgYW5vPC1yZWFkLnRhYmxlKHBhc3RlMCgiLi4vT3V0cHV0L0NPVi9DT1ZzY2FuXzNwb3AvIixwb3BzW3BdLCJfb3V0bGllcnNfIixjdltjXSwiX2Fubm90YXRpb24iKSwgaGVhZGVyID0gRikKICAgICAgICBhbm5vdGF0aW9uczwtZGF0YS5mcmFtZSgpCiAgICAgICAgZm9yIChpIGluIDE6IG5yb3coYW5vKSl7CiAgICAgICAgICAgIGFubnM8LXVubGlzdChzdHJzcGxpdChhbm8kVjRbaV0sICJcXCx8XFx8IikpCiAgICAgICAgICAgIGFubm08LWRhdGEuZnJhbWUobWF0cml4KGFubnMsbmNvbCA9IDE2LCBieXJvdyA9IFRSVUUpKQogICAgICAgICAgICBhbm5tPC1hbm5tWyxjKDIsMyw0LDUsOCldCiAgICAgICAgICAgIGNvbG5hbWVzKGFubm0pPC1jKCJFZmZlY3QiLCJQdXRhdGl2ZV9pbXBhY3QiLCJHZW5lX25hbWUiLCJHZW5lX0lEIiwiRmVhdHVyZSB0eXBlIikKICAgICAgICAgICAgYW5ubTwtYW5ubVshZHVwbGljYXRlZChhbm5tKSwgXQogICAgICAgICAgICBhbm5tJGNocjwtYW5vJFYxW2ldCiAgICAgICAgICAgIGFubm0kcG9zPC1hbm8kVjJbaV0KICAgICAgICAgICAgYW5ubSRBRjwtIGFubyRWM1tpXQogICAgICAgICAgICBhbm5vdGF0aW9uczwtcmJpbmQoYW5ub3RhdGlvbnMsIGFubm0pCiAgICAgICAgfSAgICAgCiAgICAgICAgYW5ub3RhdGlvbnM8LWFubm90YXRpb25zWyxjKDY6OCwxOjUpXQogICAgICAgIGFubm90YXRpb25zPC1hbm5vdGF0aW9uc1shZHVwbGljYXRlZChhbm5vdGF0aW9uc1ssMToyXSksXQogICAgICAgIHdyaXRlLmNzdihhbm5vdGF0aW9ucywgcGFzdGUwKCIuLi9PdXRwdXQvQ09WL0NPVnNjYW5fM3BvcC9HZW5lc18iLHBvcHNbcF0sIl9vdXRsaWVyc18xMDBrXyIsY3ZbY10sIi5jc3YiKSwgcm93Lm5hbWVzID0gRikKICAgIH0KICAgIH0KICAgIGlmIChjPT0zKXsKICAgICAgICBmb3IgKHAgaW4gMTozKXsKICAgICAgICBhbm88LXJlYWQudGFibGUocGFzdGUwKCIuLi9PdXRwdXQvQ09WL0NPVnNjYW5fM3BvcC8iLHBvcHNbcF0sIl9vdXRsaWVyc18iLGN2W2NdLCJfYW5ub3RhdGlvbiIpLCBoZWFkZXIgPSBGKQogICAgICAgIGFubm90YXRpb25zPC1kYXRhLmZyYW1lKCkKICAgICAgICBmb3IgKGkgaW4gMTogbnJvdyhhbm8pKXsKICAgICAgICAgICAgYW5uczwtdW5saXN0KHN0cnNwbGl0KGFubyRWNFtpXSwgIlxcLHxcXHwiKSkKICAgICAgICAgICAgYW5ubTwtZGF0YS5mcmFtZShtYXRyaXgoYW5ucyxuY29sID0gMTYsIGJ5cm93ID0gVFJVRSkpCiAgICAgICAgICAgIGFubm08LWFubm1bLGMoMiwzLDQsNSw4KV0KICAgICAgICAgICAgY29sbmFtZXMoYW5ubSk8LWMoIkVmZmVjdCIsIlB1dGF0aXZlX2ltcGFjdCIsIkdlbmVfbmFtZSIsIkdlbmVfSUQiLCJGZWF0dXJlIHR5cGUiKQogICAgICAgICAgICBhbm5tPC1hbm5tWyFkdXBsaWNhdGVkKGFubm0pLCBdCiAgICAgICAgICAgIGFubm0kY2hyPC1hbm8kVjFbaV0KICAgICAgICAgICAgYW5ubSRjaHI8LWFubyRWMVtpXQogICAgICAgICAgICBhbm5tJHBvczwtYW5vJFYyW2ldCiAgICAgICAgICAgIGFubm0kQUY8LSBhbm8kVjNbaV0KICAgICAgICAgICAgYW5ub3RhdGlvbnM8LXJiaW5kKGFubm90YXRpb25zLCBhbm5tKQogICAgICAgIH0gICAgIAogICAgICAgIGFubm90YXRpb25zPC1hbm5vdGF0aW9uc1ssYyg2OjgsMTo1KV0KICAgICAgICBhbm5vdGF0aW9uczwtYW5ub3RhdGlvbnNbIWR1cGxpY2F0ZWQoYW5ub3RhdGlvbnNbLDE6Ml0pLF0KICAgICAgICB3cml0ZS5jc3YoYW5ub3RhdGlvbnMsIHBhc3RlMCgiLi4vT3V0cHV0L0NPVi9DT1ZzY2FuXzNwb3AvR2VuZXNfIixwb3BzW3BdLCJfb3V0bGllcnNfMTAwa18iLGN2W2NdLCIuY3N2IiksIHJvdy5uYW1lcyA9IEYpCiAgICB9Cn0KICAKfQoKYGBgCgoKCiMjIEZpbmQgdGhlIG92ZXJsYXBwaW5nIGdlbmUgbmFtZXMgCgpgYGB7ciBldmFsPUZBTFNFLCBtZXNzYWdlPUZBTFNFLCB3YXJuaW5nPUZBTFNFfQoKZ25hbWVzZmlsZXM8LWxpc3QuZmlsZXMoIi4uL091dHB1dC9DT1YvQ09Wc2Nhbl8zcG9wLyIsIHBhdHRlcm49IkdlbmVzXy4rb3V0bGllcnNfMTAway4rXFxkLmNzdiQiKQoKZm9yIChpIGluIDE6bGVuZ3RoKGduYW1lc2ZpbGVzKSl7CiAgICBkZjwtcmVhZC5jc3YocGFzdGUwKCIuLi9PdXRwdXQvQ09WL0NPVnNjYW5fM3BvcC8iLGduYW1lc2ZpbGVzW2ldKSkKICAgIGRmPC1kZlssYygxLDY6NyldCiAgICBkZjwtZGZbIWR1cGxpY2F0ZWQoZGYpLF0KICAgIAogICAgZm5hbWU8LWdzdWIoIi5jc3YiLCIiLCBnbmFtZXNmaWxlc1tpXSkKICAgIGZuYW1lPC1nc3ViKCJHZW5lc18iLCIiLCBmbmFtZSkKICAgIAogICAgCiAgICAjYWRkIGdlbmUgbmFtZXMgZm9yIGZyb250IGFuZCBiYWNrIG9mIGludGVyZ2VuaWMgcmVnaW9ucwogICAgZGYyPC1kZltncmVwKCItIiwgZGYkR2VuZV9JRCksXQogICAgaz0xCiAgICBkZl9kaXY8LWRhdGEuZnJhbWUoKQogICAgb2RkbmFtZXM8LWRhdGEuZnJhbWUoKQogICAgZm9yIChqIGluIDE6bnJvdyhkZjIpKXsKICAgICAgICBuYW1lczwtdW5saXN0KHN0cnNwbGl0KGRmMiRHZW5lX25hbWVbal0sICItIikpCiAgICAgICAgaWRzPC11bmxpc3Qoc3Ryc3BsaXQoZGYyJEdlbmVfSURbal0sICItIikpCiAgICAgICAgCiAgICAgICAgaWYgKGxlbmd0aChuYW1lcyk9PTIpewogICAgICAgICAgICBkZl9kaXY8LXJiaW5kKGRmX2RpdiwgYyhkZjIkY2hyW2pdLG5hbWVzWzFdLGlkc1sxXSkpCiAgICAgICAgICAgIGs9aysxCiAgICAgICAgICAgIGRmX2RpdjwtcmJpbmQoZGZfZGl2LCBjKGRmMiRjaHJbal0sbmFtZXNbMl0saWRzWzJdKSkKICAgICAgICAgICAgaz1rKzEKICAgICAgICB9CiAgICAgICAKICAgICAgICBpZiAobGVuZ3RoKG5hbWVzKSE9Mil7CiAgICAgICAgICAgIG48LWdyZXAoInNpOiIsIG5hbWVzKQogICAgICAgICAgICBpZiAobGVuZ3RoKG4pPjApewogICAgICAgICAgICAgICAgaWYgKG49PTEpIG5ld25hbWVzPC1jKHBhc3RlMChuYW1lc1sxXSwiLSIsbmFtZXNbMl0pLCBuYW1lc1szXSkKICAgICAgICAgICAgICAgIGlmIChuPT0yKSBuZXduYW1lczwtYyhuYW1lc1sxXSxwYXN0ZTAobmFtZXNbMl0sIi0iLG5hbWVzWzNdKSkKICAgICAgICAgICAgICAgIGRmX2RpdjwtcmJpbmQoZGZfZGl2LCBjKGRmMiRjaHJbal0sbmV3bmFtZXNbMV0saWRzWzFdKSkKICAgICAgICAgICAgICAgIGs9aysxCiAgICAgICAgICAgICAgICBkZl9kaXY8LXJiaW5kKGRmX2RpdiwgYyhkZjIkY2hyW2pdLG5ld25hbWVzWzJdLGlkc1syXSkpCiAgICAgICAgICAgICAgICBrPWsrMQogICAgICAgICAgICB9CiAgICAgICAgICAgIAogICAgICAgICAgICBpZiAobGVuZ3RoKG4pPT0wKSB7CiAgICAgICAgICAgICAgICBvZGRuYW1lczwtcmJpbmQob2RkbmFtZXMsIGRmMltqLF0pCiAgICAgICAgICAgIH0KICAgICAgICB9CiAgICB9CiAgICBkZl9kaXY8LWRmX2RpdlshZHVwbGljYXRlZChkZl9kaXYpLF0KICAgIGRmX2RpdjwtZGZfZGl2W2RmX2RpdiRHZW5lX0lEIT0iQ0hSX0VORCIsXQogICAgZGZfZGl2PC1kZl9kaXZbZGZfZGl2JEdlbmVfSUQhPSJDSFJfU1RBUlQiLF0KICAgIAogICAgcmVtb3ZlPC1ncmVwKCItIiwgZGYkR2VuZV9JRCkKICAgIGRmPC1kZlstcmVtb3ZlLF0KICAgIGRmPC1yYmluZChkZiwgZGZfZGl2KQogICAgZGY8LWRmWyFkdXBsaWNhdGVkKGRmKSxdCiAgICAKICAgIGlmIChucm93KG9kZG5hbWVzKSE9MCl7CiAgICAgICAgd3JpdGUuY3N2KGRmLCBwYXN0ZTAoIi4uL091dHB1dC9DT1YvQ09Wc2Nhbl8zcG9wLyIsZm5hbWUsIkdlbmVMaXN0X3dpdGhJbnRlcmdlbmljR2VuZXMuY3N2IiApLCByb3cubmFtZXMgPSBGKQogICAgICAgIHdyaXRlLmNzdihvZGRuYW1lcywgcGFzdGUwKCIuLi9PdXRwdXQvQ09WL0NPVnNjYW5fM3BvcC9PZGRuYW1lc18iLCBmbmFtZSwiLmNzdiIpKQogICAgfQogICAgaWYgKG5yb3cob2RkbmFtZXMpPT0wKXsKICAgICAgICB3cml0ZS5jc3YoZGYsIHBhc3RlMCgiLi4vT3V0cHV0L0NPVi9DT1ZzY2FuXzNwb3AvIixmbmFtZSwiR2VuZUxpc3Rfd2l0aEludGVyZ2VuaWNHZW5lc19uZXcuY3N2IiApLCByb3cubmFtZXMgPSBGKQogICAgIH0KfQogICAKCiMjICEhICMjCiMjIE1hbnVhbGx5IGNoYW5nZSB0aGUgb2RkbmFtZXMgYW5kIGFkZCB0aGVtIGV0ZW8gdGhlIEdlbmVMaXN0IGZpbGVzIAojKHVwZGF0ZWQgZmlsZSBuYW1lcyBoYXMgIl9uZXciIGF0IHRoZSBlbmQpCgojYWdncmVnYXRlIGFsbCBnZW5lIG5hbWVzCmduZXc8LWxpc3QuZmlsZXMoIi4uL091dHB1dC9DT1YvQ09Wc2Nhbl8zcG9wLyIsIHBhdHRlcm49IkdlbmVMaXN0X3dpdGhJbnRlcmdlbmljR2VuZXNfbmV3LmNzdiQiKQpHZW5lczwtZGF0YS5mcmFtZSgpCkdlbmVMaXN0PC1saXN0KCkKZm9yIChpIGluIDE6bGVuZ3RoKGduZXcpKXsKICAgIGRmPC1yZWFkLmNzdihwYXN0ZTAoIi4uL091dHB1dC9DT1YvQ09Wc2Nhbl8zcG9wLyIsIGduZXdbaV0pKQogICAgR2VuZUxpc3RbW2ldXTwtZGYKICAgIGZuYW1lPC1nc3ViKCJHZW5lTGlzdF93aXRoSW50ZXJnZW5pY0dlbmVzX25ldy5jc3YiLCcnLGduZXdbaV0pCiAgICBuYW1lcyhHZW5lTGlzdClbaV08LWZuYW1lCiAgICBkdXA8LWRmW2R1cGxpY2F0ZWQoZGYpLF0KICAgIGRmPC1kZlshZHVwbGljYXRlZChkZiksXQogICAgR2VuZXM8LXJiaW5kKEdlbmVzLCBkZikKICAgIEdlbmVzPC1HZW5lc1shZHVwbGljYXRlZChHZW5lcyksXQogICAgCn0KCgojMS4gQmV0d2VlbiBwb3B1bGF0aW9ucwp0aW1lczwtYygiY292MTIiLCJjb3YxMyIsImNvdjIzIikKY29tbW9uPC1saXN0KCkKY29tbW9uX3N1bW1hcnk8LWRhdGEuZnJhbWUodGltZT10aW1lcykKZm9yIChpIGluIDE6Myl7CiAgICB0bGlzdDwtR2VuZUxpc3RbZ3JlcCh0aW1lc1tpXSwgbmFtZXMoR2VuZUxpc3QpKV0KICAgIGlmIChpICE9Myl7CiAgICAgICAgY29tbW9uX2dlbmVzPC1pbnRlcnNlY3QodGxpc3RbWzFdXVsiR2VuZV9uYW1lIl0sIHRsaXN0W1syXV1bIkdlbmVfbmFtZSJdKQogICAgICAgIGNvbW1vbltbaV1dPC1jb21tb25fZ2VuZXMKICAgICAgICBuYW1lcyhjb21tb24pW1tpXV08LXRpbWVzW2ldCiAgICAgICAgY29tbW9uX3N1bW1hcnkkUFdTW2ldPC1ucm93KHRsaXN0W1tncmVwKCJQV1MiLCBuYW1lcyh0bGlzdCkpXV0pCiAgICAgICAgY29tbW9uX3N1bW1hcnkkVEJbaV08LW5yb3codGxpc3RbW2dyZXAoIlRCIiwgbmFtZXModGxpc3QpKV1dKQogICAgICAgIGNvbW1vbl9zdW1tYXJ5JFNTW2ldPC1OQQogICAgICAgIGNvbW1vbl9zdW1tYXJ5JGNvbW1vbl9QV1MuVEJbaV08LW5yb3coY29tbW9uX2dlbmVzKQogICAgICAgIAogICAgICAgIHB3czwtdGxpc3RbWzFdXVsiR2VuZV9uYW1lIl0KICAgICAgICB0YjwtdGxpc3RbWzJdXVsiR2VuZV9uYW1lIl0KICAgICAgICB4PC1saXN0KFBXUz1wd3MkR2VuZV9uYW1lLFRCPXRiJEdlbmVfbmFtZSkKICAgICAgICBnZ3Zlbm4oeCwgZmlsbF9jb2xvciA9IGNvbHNbYygyLDEpXSwgc3Ryb2tlX3NpemUgPSAwLjUsIHNldF9uYW1lX3NpemUgPSA0LHRleHRfc2l6ZT0zKStnZ3RpdGxlKHRpbWVzW2ldKQogICAgICAgIGdnc2F2ZShwYXN0ZTAoIi4uL091dHB1dC9DT1YvQ09Wc2Nhbl8zcG9wL1Zlbm5fIix0aW1lc1tpXSwiLnBuZyIpLCB3aWR0aCA9IDMsIGhlaWdodD0zLCBkcGk9MzAwKQogICAgfQogICAgaWYgKGk9PTMpewogICAgICAgIGNvbW1vbl9zdW1tYXJ5JFBXU1tpXTwtbnJvdyh0bGlzdFtbZ3JlcCgiUFdTIiwgbmFtZXModGxpc3QpKV1dKQogICAgICAgIGNvbW1vbl9zdW1tYXJ5JFRCW2ldPC0gbnJvdyh0bGlzdFtbZ3JlcCgiVEIiLCBuYW1lcyh0bGlzdCkpXV0pCiAgICAgICAgY29tbW9uX3N1bW1hcnkkU1NbaV08LSBucm93KHRsaXN0W1tncmVwKCJTUyIsIG5hbWVzKHRsaXN0KSldXSkKICAgICAgICAKICAgICAgICBnZW5lczE8LWludGVyc2VjdCh0bGlzdFtbMV1dWyJHZW5lX25hbWUiXSwgdGxpc3RbWzNdXVsiR2VuZV9uYW1lIl0pCiAgICAgICAgZ2VuZXMyPC1pbnRlcnNlY3QodGxpc3RbWzFdXVsiR2VuZV9uYW1lIl0sIHRsaXN0W1syXV1bIkdlbmVfbmFtZSJdKQogICAgICAgIGdlbmVzMzwtaW50ZXJzZWN0KHRsaXN0W1syXV1bIkdlbmVfbmFtZSJdLCB0bGlzdFtbM11dWyJHZW5lX25hbWUiXSkKICAgICAgICBnZW5lczQ8LWludGVyc2VjdCh0bGlzdFtbMV1dWyJHZW5lX25hbWUiXSxpbnRlcnNlY3QodGxpc3RbWzJdXVsiR2VuZV9uYW1lIl0sIHRsaXN0W1szXV1bIkdlbmVfbmFtZSJdKSkKICAgICAgICBjb21tb25fc3VtbWFyeSRjb21tb25fUFdTLlRCW2ldPC1ucm93KGdlbmVzMSkKICAgICAgICBjb21tb25fc3VtbWFyeSRjb21tb25fUFdTLlNTW2ldPC1ucm93KGdlbmVzMikKICAgICAgICBjb21tb25fc3VtbWFyeSRjb21tb25fU1MuVEJbaV08LW5yb3coZ2VuZXMzKQogICAgICAgIGNvbW1vbl9zdW1tYXJ5JGNvbW1vbjNbaV08LW5yb3coZ2VuZXM0KQogICAgICAgIGs9aQogICAgICAgIGNvbW1vbltba11dPC1nZW5lczIKICAgICAgICBuYW1lcyhjb21tb24pW1trXV08LXBhc3RlMCh0aW1lc1tpXSwiX1BXUy5TUyIpCiAgICAgICAgaz1rKzEKICAgICAgICBjb21tb25bW2tdXTwtZ2VuZXMxCiAgICAgICAgbmFtZXMoY29tbW9uKVtba11dPC1wYXN0ZTAodGltZXNbaV0sIl9QV1MuVEIiKQogICAgICAgIGs9aysxCiAgICAgICAgY29tbW9uW1trXV08LWdlbmVzMwogICAgICAgIG5hbWVzKGNvbW1vbilbW2tdXTwtcGFzdGUwKHRpbWVzW2ldLCJfU1MuVEIiKQogICAgICAgIGs9aysxCiAgICAgICAgY29tbW9uW1trXV08LWdlbmVzNAogICAgICAgIG5hbWVzKGNvbW1vbilbW2tdXTwtcGFzdGUwKHRpbWVzW2ldLCJfM3BvcHMiKQogICAgICAgIAogICAgICAgIHB3czwtdGxpc3RbWzFdXVsiR2VuZV9uYW1lIl0KICAgICAgICB0YjwtdGxpc3RbWzNdXVsiR2VuZV9uYW1lIl0KICAgICAgICBzczwtdGxpc3RbWzJdXVsiR2VuZV9uYW1lIl0KICAgICAgICB4PC1saXN0KFBXUz1wd3MkR2VuZV9uYW1lLFRCPXRiJEdlbmVfbmFtZSwgU1M9c3MkR2VuZV9uYW1lKQogICAgICAgIGdndmVubih4LCBmaWxsX2NvbG9yID0gY29sc1tjKDIsMSwzKV0sIHN0cm9rZV9zaXplID0gMC41LCBzZXRfbmFtZV9zaXplID0gNCx0ZXh0X3NpemU9MykrZ2d0aXRsZSh0aW1lc1tpXSkKICAgICAgICBnZ3NhdmUocGFzdGUwKCIuLi9PdXRwdXQvQ09WL0NPVnNjYW5fM3BvcC9WZW5uXyIsdGltZXNbaV0sIi5wbmciKSwgd2lkdGggPSA0LCBoZWlnaHQ9NCwgZHBpPTMwMCkKICAgICAgICAKICAgICAgICAgeDE8LWxpc3QoUFdTPXB3cyRHZW5lX25hbWUsVEI9dGIkR2VuZV9uYW1lKQogICAgICAgIGdndmVubih4MSwgZmlsbF9jb2xvciA9IGNvbHNbYygyLDEpXSwgc3Ryb2tlX3NpemUgPSAwLjUsIHNldF9uYW1lX3NpemUgPSA0LHRleHRfc2l6ZT0zKStnZ3RpdGxlKHRpbWVzW2ldKQogICAgICAgIGdnc2F2ZShwYXN0ZTAoIi4uL091dHB1dC9DT1YvQ09Wc2Nhbl8zcG9wL1Zlbm5fUFdTX1RCXyIsdGltZXNbaV0sIi5wbmciKSwgd2lkdGggPSAzLCBoZWlnaHQ9MywgZHBpPTMwMCkKICAgICAgICAgeDI8LWxpc3QoUFdTPXB3cyRHZW5lX25hbWUsU1M9c3MkR2VuZV9uYW1lKQogICAgICAgIGdndmVubih4MiwgZmlsbF9jb2xvciA9IGNvbHNbYygyLDMpXSwgc3Ryb2tlX3NpemUgPSAwLjUsIHNldF9uYW1lX3NpemUgPSA0LHRleHRfc2l6ZT0zKStnZ3RpdGxlKHRpbWVzW2ldKQogICAgICAgIGdnc2F2ZShwYXN0ZTAoIi4uL091dHB1dC9DT1YvQ09Wc2Nhbl8zcG9wL1Zlbm5fUFdTX1NTXyIsdGltZXNbaV0sIi5wbmciKSwgd2lkdGggPSAzLCBoZWlnaHQ9MywgZHBpPTMwMCkKICAgICAgICAgIHgzPC1saXN0KFNTPXNzJEdlbmVfbmFtZSwgVEI9dGIkR2VuZV9uYW1lKQogICAgICAgIGdndmVubih4MywgZmlsbF9jb2xvciA9IGNvbHNbYygzLDEpXSwgc3Ryb2tlX3NpemUgPSAwLjUsIHNldF9uYW1lX3NpemUgPSA0LHRleHRfc2l6ZT0zKStnZ3RpdGxlKHRpbWVzW2ldKQogICAgICAgIGdnc2F2ZShwYXN0ZTAoIi4uL091dHB1dC9DT1YvQ09Wc2Nhbl8zcG9wL1Zlbm5fU1NfVEJfIix0aW1lc1tpXSwiLnBuZyIpLCB3aWR0aCA9IDMsIGhlaWdodD0zLCBkcGk9MzAwKQogICAgICAgIAogICAgICAgIAogICAgICAgIH0KfQp3cml0ZS5jc3YoY29tbW9uX3N1bW1hcnksICIuLi9PdXRwdXQvQ09WL0NPVnNjYW5fM3BvcC9Db21tb25fZ2VuZXNfd2l0aEludGVyZ2VuZXNfM3BvcHMuY3N2IikKCgojV2hhdCBhcmUgdGhlIG92ZXJsYXBwaW5nIGdlbmUgbmFtZXMgYmV0d2VlbiBwb3B1bGF0aW9ucwpjb21tb25fdGltZXM8LWxpc3QoKQpmb3IgKGkgaW4gMTogbGVuZ3RoKGNvbW1vbikpewogICAgZ2lkczwtY29tbW9uW1tpXV0KICAgIGRmPC1kYXRhLmZyYW1lKEdlbmVfbmFtZT1naWRzKQogICAgCiAgICBkZjI8LW1lcmdlKGRmLCBHZW5lcywgYnk9IkdlbmVfbmFtZSIpCiAgICB3cml0ZS5jc3YoZGYyLCBwYXN0ZTAoIi4uL091dHB1dC9DT1YvQ09Wc2Nhbl8zcG9wL0NvbW1vbl9nZW5lc18iLCBuYW1lcyhjb21tb24pW2ldLCIuY3N2IiksIHJvdy5uYW1lcyA9IEYpCiAgICBjb21tb25fdGltZXNbW2ldXTwtZGYyCiAgICBuYW1lcyhjb21tb25fdGltZXMpW2ldPC0gbmFtZXMoY29tbW9uKVtpXQp9Cgojbm9uX292ZXJsYXBwaW5nIGdlbmVzIENPVjIzCgp0bGlzdDwtR2VuZUxpc3RbZ3JlcCh0aW1lc1szXSwgbmFtZXMoR2VuZUxpc3QpKV0KcHdzMjM8LXRsaXN0W1sxXV1bIkdlbmVfSUQiXQpzczIzPC10bGlzdFtbMl1dWyJHZW5lX0lEIl0KdGIyMzwtdGxpc3RbWzNdXVsiR2VuZV9JRCJdCmdlbmVzMTwtaW50ZXJzZWN0KHRsaXN0W1sxXV1bIkdlbmVfSUQiXSwgdGxpc3RbWzNdXVsiR2VuZV9JRCJdKQpnZW5lczI8LWludGVyc2VjdCh0bGlzdFtbMV1dWyJHZW5lX0lEIl0sIHRsaXN0W1syXV1bIkdlbmVfSUQiXSkKZ2VuZXMzPC1pbnRlcnNlY3QodGxpc3RbWzJdXVsiR2VuZV9JRCJdLCB0bGlzdFtbM11dWyJHZW5lX0lEIl0pCmdlbmVzNDwtaW50ZXJzZWN0KHRsaXN0W1sxXV1bIkdlbmVfSUQiXSxpbnRlcnNlY3QodGxpc3RbWzJdXVsiR2VuZV9JRCJdLCB0bGlzdFtbM11dWyJHZW5lX0lEIl0pKQoKb3ZlcmxhcHM8LXJiaW5kKGdlbmVzMSwgZ2VuZXMyLCBnZW5lczMsIGdlbmVzNCkKb3ZlcmxhcHM8LW92ZXJsYXBzWyFkdXBsaWNhdGVkKG92ZXJsYXBzJEdlbmVfSUQpLF0KcHdzMjNfb25seTwtZGF0YS5mcmFtZShHZW5lX0lEPXB3czIzWyEocHdzMjMkR2VuZV9JRCAlaW4lIG92ZXJsYXBzKSwgXSkKZzE8LXVuaXF1ZShwd3MyM19vbmx5JEdlbmVfSUQpCnNpbmsoIi4uL091dHB1dC9DT1YvQ09Wc2Nhbl8zcG9wL2dlbmVzL0dlbmVJRF9QV1Nvbmx5X0NPVjIzLnR4dCIpCmNhdChwYXN0ZTAoZzEsICI7IikpCnNpbmsoTlVMTCkKICAgICAgICAKdGIyM19vbmx5PC10YjIzWyEodGIyMyRHZW5lX0lEICVpbiUgb3ZlcmxhcHMpLCBdCmcxPC11bmlxdWUodGIyM19vbmx5KQpzaW5rKCIuLi9PdXRwdXQvQ09WL0NPVnNjYW5fM3BvcC9nZW5lcy9HZW5lSURfVEJvbmx5X0NPVjIzLnR4dCIpCmNhdChwYXN0ZTAoZzEsICI7IikpCnNpbmsoTlVMTCkKCnNzMjNfb25seTwtc3MyM1shKHNzMjMkR2VuZV9JRCAlaW4lIG92ZXJsYXBzKSwgXQpnMTwtdW5pcXVlKHNzMjNfb25seSkKc2luaygiLi4vT3V0cHV0L0NPVi9DT1ZzY2FuXzNwb3AvZ2VuZXMvR2VuZUlEX1NTb25seV9DT1YyMy50eHQiKQpjYXQocGFzdGUwKGcxLCAiOyIpKQpzaW5rKE5VTEwpCgoKCgojMi4gQmV0d2VlbiBUaW1lLVBvaW50cyB3aXRoaW4gYSBwb3B1bGF0aW9uCgp0aW1lczwtYygiY292MTIiLCJjb3YxMyIsImNvdjIzIikKcG9wczwtYygiUFdTIiwiVEIiKQpjb21tb24yPC1saXN0KCkKY29tbW9uX3N1bW1hcnkyPC1kYXRhLmZyYW1lKHBvcD1yZXAocG9wc1sxOjJdLCBlYWNoPTQpKQpmb3IgKGkgaW4gMTpsZW5ndGgocG9wcykpewogICAgcGxpc3Q8LUdlbmVMaXN0W2dyZXAocG9wc1tpXSwgbmFtZXMoR2VuZUxpc3QpKV0KICAgIGs9NCppLTMKICAgICNjb21tb24gZ2VuZXMgYmV0d2VlbiBDT1YxMiBhbmQgQ09WMTMKICAgIGNvbW1vbl9nZW5lczE8LWludGVyc2VjdChwbGlzdFtbMV1dWyJHZW5lX25hbWUiXSwgcGxpc3RbWzJdXVsiR2VuZV9uYW1lIl0pCiAgICBjb21tb24yW1trXV08LWNvbW1vbl9nZW5lczEKICAgIG5hbWVzKGNvbW1vbjIpW1trXV08LXBhc3RlMChwb3BzW2ldLCIuIiwgdGltZXNbMV0sIl8iLHRpbWVzWzJdKQogICAgY29tbW9uX3N1bW1hcnkyJFRpbWVba108LXBhc3RlMCh0aW1lc1sxXSwiXyIsdGltZXNbMl0pCiAgICBjb21tb25fc3VtbWFyeTIkbm8ub2YuZ2VuZXNba108LW5yb3coY29tbW9uX2dlbmVzMSkgCiAgICAKICAgIGMxMjwtcGxpc3RbWzFdXVsiR2VuZV9uYW1lIl0KICAgIGMxMzwtcGxpc3RbWzJdXVsiR2VuZV9uYW1lIl0KICAgIGMyMzwtcGxpc3RbWzNdXVsiR2VuZV9uYW1lIl0KICAgIHg8LWxpc3QoQ09WMTI9YzEyJEdlbmVfbmFtZSxDT1YxMz1jMTMkR2VuZV9uYW1lLCBDT1YyMz1jMjMkR2VuZV9uYW1lKQogICAgZ2d2ZW5uKHgsIGZpbGxfY29sb3IgPSBjb2xzW2MoMSw1LDcpXSwgc3Ryb2tlX3NpemUgPSAwLjUsIHNldF9uYW1lX3NpemUgPSA0LHRleHRfc2l6ZT0zKStnZ3RpdGxlKHBvcHNbaV0pCiAgICBnZ3NhdmUocGFzdGUwKCIuLi9PdXRwdXQvQ09WL0NPVnNjYW5fM3BvcC9WZW5uXyIscG9wc1tpXSwiLnBuZyIpLCB3aWR0aCA9IDQsIGhlaWdodD00LCBkcGk9MzAwKQoKICAgIAogICAgaz1rKzEKICAgICNjb21tb24gZ2VuZXMgYmV0d2VlbiBDT1YxMiBhbmQgQ09WMjMKICAgIGNvbW1vbl9nZW5lczI8LWludGVyc2VjdChwbGlzdFtbMV1dWyJHZW5lX25hbWUiXSwgcGxpc3RbWzNdXVsiR2VuZV9uYW1lIl0pCiAgICBjb21tb24yW1trXV08LWNvbW1vbl9nZW5lczIKICAgIG5hbWVzKGNvbW1vbjIpW1trXV08LXBhc3RlMChwb3BzW2ldLCIuIiwgdGltZXNbMV0sIl8iLHRpbWVzWzNdKQogICAgY29tbW9uX3N1bW1hcnkyJFRpbWVba108LXBhc3RlMCh0aW1lc1sxXSwiXyIsdGltZXNbM10pCiAgICBjb21tb25fc3VtbWFyeTIkbm8ub2YuZ2VuZXNba108LW5yb3coY29tbW9uX2dlbmVzMikgCiAKICAgIGs9aysxCiAgICAjY29tbW9uIGdlbmVzIGJldHdlZW4gQ09WMTMgYW5kIENPVjIzCiAgICBjb21tb25fZ2VuZXMzPC1pbnRlcnNlY3QocGxpc3RbWzJdXVsiR2VuZV9uYW1lIl0sIHBsaXN0W1szXV1bIkdlbmVfbmFtZSJdKQogICAgY29tbW9uMltba11dPC1jb21tb25fZ2VuZXMzCiAgICBuYW1lcyhjb21tb24yKVtba11dPC1wYXN0ZTAocG9wc1tpXSwiLiIsIHRpbWVzWzJdLCJfIix0aW1lc1szXSkKICAgIGNvbW1vbl9zdW1tYXJ5MiRUaW1lW2tdPC1wYXN0ZTAodGltZXNbMl0sIl8iLHRpbWVzWzNdKQogICAgY29tbW9uX3N1bW1hcnkyJG5vLm9mLmdlbmVzW2tdPC1ucm93KGNvbW1vbl9nZW5lczMpIAogCiAgICBrPWsrMQogICAgI2NvbW1vbiBnZW5lcyBhbW9uZyBhbGwgdGltZSBwZXJpb2RzCiAgICBjb21tb25fZ2VuZXM0PC1pbnRlcnNlY3QocGxpc3RbWzFdXVsiR2VuZV9uYW1lIl0sIChpbnRlcnNlY3QocGxpc3RbWzJdXVsiR2VuZV9uYW1lIl0sIHBsaXN0W1szXV1bIkdlbmVfbmFtZSJdKSkpCiAgICBjb21tb24yW1trXV08LWNvbW1vbl9nZW5lczQKICAgIG5hbWVzKGNvbW1vbjIpW1trXV08LXBhc3RlMChwb3BzW2ldLCIuYWxsIikKICAgIGNvbW1vbl9zdW1tYXJ5MiRUaW1lW2tdPC0iQWxsIgogICAgY29tbW9uX3N1bW1hcnkyJG5vLm9mLmdlbmVzW2tdPC1ucm93KGNvbW1vbl9nZW5lczQpIAp9CndyaXRlLmNzdihjb21tb25fc3VtbWFyeTIsICIuLi9PdXRwdXQvQ09WL0NPVnNjYW5fM3BvcC9Db21tb25fZ2VuZXNfYmV0d2VlblRpbWVQb2ludHMuY3N2IikKCgojQ29tbW9uIGdlbmUgbmFtZXMgYmV0d2VlbiB0aW1lIHBvaW50cwoKZm9yIChpIGluIDE6Mil7CiAgICBDb21tb25HZW5lczwtZGF0YS5mcmFtZSgpCiAgICBnbGlzdDwtY29tbW9uMltncmVwKHBvcHNbaV0sIG5hbWVzKGNvbW1vbjIpKV0KICAgIGZvcihqIGluIDE6bGVuZ3RoKGdsaXN0KSl7CiAgICAgICAgZ2lkczwtZ2xpc3RbW2pdXQogICAgICAgIGRmPC1kYXRhLmZyYW1lKEdlbmVfbmFtZT1naWRzKQogICAgICAgIGRmMjwtbWVyZ2UoZGYsIEdlbmVzLCBieT0iR2VuZV9uYW1lIiwgYWxsLng9VCkKICAgICAgICB3cml0ZS5jc3YoZGYyLCBwYXN0ZTAoIi4uL091dHB1dC9DT1YvQ09Wc2Nhbl8zcG9wL0NvbW1vbl9nZW5lc18iLCBuYW1lcyhnbGlzdClbal0sIi5jc3YiKSwgcm93Lm5hbWVzID0gRikKICAgICAgICBkZjIkVGltZTwtbmFtZXMoZ2xpc3QpW2pdCiAgICAgICAgQ29tbW9uR2VuZXM8LXJiaW5kKENvbW1vbkdlbmVzLCBkZjIpCiAgICB9CiAgICB3cml0ZS5jc3YoQ29tbW9uR2VuZXMsIHBhc3RlMCgiLi4vT3V0cHV0L0NPVi9DT1ZzY2FuXzNwb3AvQ29tbW9uX2dlbmVzXyIscG9wc1tpXSAsIi5jc3YiKSwgcm93Lm5hbWVzID0gRikKfQoKCgpgYGAKCgoKIyMjIyBPdmVybGFwcGluZyBnZW5lIG51bWJlcnMgICAKYGBge3IgZWNobz1GQUxTRSwgbWVzc2FnZT1GQUxTRSwgd2FybmluZz1GQUxTRX0KCiMgU3VtbWFyeSB0YWJsZQpjb21tb25fZ2VuZXM8LXJlYWQuY3N2KCIuLi9PdXRwdXQvQ09WL0NPVnNjYW5fM3BvcC9Db21tb25fZ2VuZXNfM3BvcHMuY3N2Iiwgcm93Lm5hbWVzID0gMSkKa25pdHI6OmthYmxlKGNvbW1vbl9nZW5lcykKCmBgYAoKCiMjIFdoYXQgYXJlIHRoZSBnZW5lcyBvdmVybGFwcGluZyBhY3Jvc3MgZGlmZmVyZW50IHRpbWUgcG9pbnRzIGJldHdlZW4gcG9wdWxhdGlvbnM/IAoKYGBge3IgZXZhbD1GQUxTRSwgbWVzc2FnZT1GQUxTRSwgd2FybmluZz1GQUxTRX0KIzEuIGJldHdlZW4gUFdTIGFuZCBUQgpwd3MudGI8LWNvbW1vbl90aW1lc1tjKDEsMiw0KV0KCmdlbmVzMTIxMzwtaW50ZXJzZWN0KHB3cy50YltbMV1dWyJHZW5lX25hbWUiXSwgcHdzLnRiW1syXV1bIkdlbmVfbmFtZSJdKQpnZW5lczEyMTM8LW1lcmdlKGdlbmVzMTIxMywgR2VuZXMsIGJ5PSJHZW5lX25hbWUiKQp3cml0ZS5jc3YoZ2VuZXMxMjIzLCAiLi4vT3V0cHV0L0NPVi9DT1ZzY2FuXzNwb3AvQ29tbW9uX2dlbmVzX1BXUy5UQi5jb3YxMi1jb3YyMy5jc3YiKQojICAgICAgICAgICBHZW5lX25hbWUgICBjaHIgICAgICAgICAgICBHZW5lX0lECiMxIEVOU0NIQUcwMDAwMDAwMTY4NyBjaHIxMyBFTlNDSEFHMDAwMDAwMDE2ODcKIzIgICAgICAgICAgICAgbmRzdDJhIGNocjEzIEVOU0NIQUcwMDAwMDAwMjY0OQojMyAgICAgICAgICAgICB6c3dpbTggY2hyMTMgRU5TQ0hBRzAwMDAwMDA1OTU2CgojIyBDb21tb24gZ2VuZXMgYmV0d2VlbiBwb3B1bGF0aW9uIGFjcm9zcyB0aW1lIHBvaW50cyAoQ09WMTIgLSBDT1YxMykKcDEyMTM8LXB3cy50YltbMV1dWyJHZW5lX25hbWUiXQp0MTIxMzwtcHdzLnRiW1syXV1bIkdlbmVfbmFtZSJdCng8LWxpc3QoUFdTPXAxMjEzJEdlbmVfbmFtZSxUQj10MTIxMyRHZW5lX25hbWUpCmdndmVubih4LCBmaWxsX2NvbG9yID0gY29sc1tjKDIsMSldLCBzdHJva2Vfc2l6ZSA9IDAuNSwgc2V0X25hbWVfc2l6ZSA9IDQsdGV4dF9zaXplPTMpK2dndGl0bGUoIkNPVjEyLUNPVjEzIGluIFBXUyAmIFRCIikKZ2dzYXZlKHBhc3RlMCgiLi4vT3V0cHV0L0NPVi9DT1ZzY2FuXzNwb3AvVmVubl9QV1NfVEJfQ09WMTItQ09WMTMucG5nIiksIHdpZHRoID0gMywgaGVpZ2h0PTMsIGRwaT0zMDApCiAgICAgICAgCmdlbmVzMTIyMzwtaW50ZXJzZWN0KHB3cy50YltbMV1dWyJHZW5lX25hbWUiXSwgcHdzLnRiW1szXV1bIkdlbmVfbmFtZSJdKQpnZW5lczEyMjM8LW1lcmdlKGdlbmVzMTIyMywgR2VuZXMsIGJ5PSJHZW5lX25hbWUiKQp3cml0ZS5jc3YoZ2VuZXMxMjIzLCAiLi4vT3V0cHV0L0NPVi9DT1ZzY2FuXzNwb3AvQ29tbW9uX2dlbmVzX1BXUy5UQi5jb3YxMi1jb3YxMy5jc3YiKQojICAgICAgICAgICBHZW5lX25hbWUgICBjaHIgICAgICAgICAgICBHZW5lX0lECiMxIEVOU0NIQUcwMDAwMDAwMTY4NyBjaHIxMyBFTlNDSEFHMDAwMDAwMDE2ODcKIzIgRU5TQ0hBRzAwMDAwMDIyNzA5IGNocjIwIEVOU0NIQUcwMDAwMDAyMjcwOQojMyBFTlNDSEFHMDAwMDAwMjI4MTUgY2hyMjAgRU5TQ0hBRzAwMDAwMDIyODE1CiM0ICAgICAgICAgICAgIG5kc3QyYSBjaHIxMyBFTlNDSEFHMDAwMDAwMDI2NDkKCnAxMjIzPC1wd3MudGJbWzFdXVsiR2VuZV9uYW1lIl0KdDEyMjM8LXB3cy50YltbM11dWyJHZW5lX25hbWUiXQp4PC1saXN0KFBXUz1wMTIyMyRHZW5lX25hbWUsVEI9dDEyMjMkR2VuZV9uYW1lKQpnZ3Zlbm4oeCwgZmlsbF9jb2xvciA9IGNvbHNbYygyLDEpXSwgc3Ryb2tlX3NpemUgPSAwLjUsIHNldF9uYW1lX3NpemUgPSA0LHRleHRfc2l6ZT0zKStnZ3RpdGxlKCJDT1YxMi1DT1YyMyBpbiBQV1MgJiBUQiIpCmdnc2F2ZShwYXN0ZTAoIi4uL091dHB1dC9DT1YvQ09Wc2Nhbl8zcG9wL1Zlbm5fUFdTX1RCX0NPVjEyLUNPVjIzLnBuZyIpLCB3aWR0aCA9IDMsIGhlaWdodD0zLCBkcGk9MzAwKQoKCmdlbmVzMTMyMzwtaW50ZXJzZWN0KHB3cy50YltbMl1dWyJHZW5lX25hbWUiXSwgcHdzLnRiW1szXV1bIkdlbmVfbmFtZSJdKQpnZW5lczEzMjM8LW1lcmdlKGdlbmVzMTMyMywgR2VuZXMsIGJ5PSJHZW5lX25hbWUiKQp3cml0ZS5jc3YoZ2VuZXMxMzIzLCAiLi4vT3V0cHV0L0NPVi9DT1ZzY2FuXzNwb3AvQ29tbW9uX2dlbmVzX1BXUy5UQi5jb3YxMy1jb3YyMy5jc3YiKQojICAgICAgICAgICBHZW5lX25hbWUgICBjaHIgICAgICAgICAgICBHZW5lX0lECiMxIEVOU0NIQUcwMDAwMDAwMTY4NyBjaHIxMyBFTlNDSEFHMDAwMDAwMDE2ODcKIzIgICAgICAgICAgICAgbmRzdDJhIGNocjEzIEVOU0NIQUcwMDAwMDAwMjY0OQoKcDEzMjM8LXB3cy50YltbMl1dWyJHZW5lX25hbWUiXQp0MTMyMzwtcHdzLnRiW1szXV1bIkdlbmVfbmFtZSJdCng8LWxpc3QoUFdTPXAxMzIzJEdlbmVfbmFtZSxUQj10MTMyMyRHZW5lX25hbWUpCmdndmVubih4LCBmaWxsX2NvbG9yID0gY29sc1tjKDIsMSldLCBzdHJva2Vfc2l6ZSA9IDAuNSwgc2V0X25hbWVfc2l6ZSA9IDQsdGV4dF9zaXplPTMpK2dndGl0bGUoIkNPVjEzLUNPVjIzIGluIFBXUyAmIFRCIikKZ2dzYXZlKHBhc3RlMCgiLi4vT3V0cHV0L0NPVi9DT1ZzY2FuXzNwb3AvVmVubl9QV1NfVEJfQ09WMTMtQ09WMjMucG5nIiksIHdpZHRoID0gMywgaGVpZ2h0PTMsIGRwaT0zMDApCgoKCgpgYGAKCgoKIyBDb21wYXJlIHJlc3VsdHMgZnJvbSBQSF9NRDcwMDAgKGFsbCBwb3BzIHRvZ2V0aGVyKSBhbmQgM1BvcHMtTlMwLjUgKGZpbHRlciBmb3IgMyBwb3BzKSBWQ0YgZmlsZXMgIAoKYGBge3IgZXZhbD1GQUxTRSwgbWVzc2FnZT1GQUxTRSwgd2FybmluZz1GQUxTRX0KcHdzMTwtcmVhZC5jc3YoIi4uL091dHB1dC9DT1YvQ09Wc2Nhbl8zcG9wLzNwb3BzX3RvcDFwZXJjZW50X291dGxpZXJfcmVnaW9ucy5jb3YxMl9uZXcuY3N2IikKcHdzMTwtcHdzMVtwd3MxJHBvcD09IlBXUyIgJiBwd3MxJHdpbmRvdz09IjEwMGsiLF0KcHdzMjwtcmVhZC5jc3YoIi4uL091dHB1dC9DT1YvM3BvcHNfdG9wMXBlcmNlbnRfb3V0bGllcl9yZWdpb25zLmNvdjEyLmNzdiIpCnB3czI8LXB3czJbcHdzMiRwb3A9PSJQV1MiJnB3czIkd2luZG93PT0iMTAwayIsXQoKcHdzMSR2Y2Y8LSIzUG9wcyIKcHdzMiR2Y2Y8LSJQSCIKcHdzPC1yYmluZChwd3MxWyxjKCJjaHJvbSIsInN0YXJ0IiwiZW5kIiwiY292MTIiLCJ2Y2YiKV0scHdzMlssYygiY2hyb20iLCJzdGFydCIsImVuZCIsImNvdjEyIiwidmNmIildKQpwd3MkY2hyb208LWZhY3Rvcihwd3MkY2hyb20sIGxldmVscz1wYXN0ZTAoImNociIsIDE6MjYpKQpnZ3Bsb3QoZGF0YT1wd3MsYWVzKHg9c3RhcnQvMTAwMDAwMCwgeT1jb3YxMiwgY29sb3I9dmNmLCBmaWxsPXZjZikpKwogICAgZmFjZXRfd3JhcCh+Y2hyb20pKwogICAgZ2VvbV9wb2ludChwb3NpdGlvbj1wb3NpdGlvbl9kb2RnZSh3aWR0aCA9IDAuNyksIGFscGhhPTAuNSkreGxhYigiUG9zaXRpb24gKE1iKSIpKwogICAgZ2d0aXRsZSgiUFdTIENPVjEyIikKZ2dzYXZlKCIuLi9PdXRwdXQvQ09WL0NPVnNjYW5fM3BvcC9QV1NfUEgudnMuM1BvcHNfb3V0bGllcl9vdmVybGFwX2NvdjEyLnBuZyIsIHdpZHRoID0gMTAsIGhlaWdodCA9IDgsIGRwaT0zMDApCmBgYAohW10oLi4vT3V0cHV0L0NPVi9DT1ZzY2FuXzNwb3AvUFdTX1BILnZzLjNQb3BzX291dGxpZXJfb3ZlcmxhcF9jb3YxMi5wbmcpCgpgYGB7ciBldmFsPUZBTFNFLCBtZXNzYWdlPUZBTFNFLCB3YXJuaW5nPUZBTFNFfQoKcHdzMTwtcmVhZC5jc3YoIi4uL091dHB1dC9DT1YvQ09Wc2Nhbl8zcG9wLzNwb3BzX3RvcDFwZXJjZW50X291dGxpZXJfcmVnaW9ucy5jb3YyM19uZXcuY3N2IikKcHdzMTwtcHdzMVtwd3MxJHBvcD09IlBXUyIgJiBwd3MxJHdpbmRvdz09IjEwMGsiLF0KcHdzMjwtcmVhZC5jc3YoIi4uL091dHB1dC9DT1YvM3BvcHNfdG9wMXBlcmNlbnRfb3V0bGllcl9yZWdpb25zLmNvdjIzLmNzdiIpCnB3czI8LXB3czJbcHdzMiRwb3A9PSJQV1MiJnB3czIkd2luZG93PT0iMTAwayIsXQoKcHdzMSR2Y2Y8LSJQSCIKcHdzMiR2Y2Y8LSIzUG9wcyIKcHdzPC1yYmluZChwd3MxWyxjKCJjaHJvbSIsInN0YXJ0IiwiZW5kIiwiY292MjMiLCJ2Y2YiKV0scHdzMlssYygiY2hyb20iLCJzdGFydCIsImVuZCIsImNvdjIzIiwidmNmIildKQpwd3MkY2hyb208LWZhY3Rvcihwd3MkY2hyb20sIGxldmVscz1wYXN0ZTAoImNociIsIDE6MjYpKQpnZ3Bsb3QoZGF0YT1wd3MsYWVzKHg9c3RhcnQvMTAwMDAwMCwgeT1jb3YyMywgY29sb3I9dmNmLCBmaWxsPXZjZikpKwogICAgZmFjZXRfd3JhcCh+Y2hyb20pKwogICAgZ2VvbV9wb2ludChwb3NpdGlvbj1wb3NpdGlvbl9kb2RnZSh3aWR0aCA9IDAuNyksIGFscGhhPTAuNSkreGxhYigiUG9zaXRpb24gKE1iKSIpKwogICAgZ2d0aXRsZSgiUFdTIENPVjIzIikKZ2dzYXZlKCIuLi9PdXRwdXQvQ09WL0NPVnNjYW5fM3BvcC9QV1NfUEgudnMuM1BvcHNfb3V0bGllcl9vdmVybGFwX2NvdjIzLnBuZyIsIHdpZHRoID0gMTAsIGhlaWdodCA9IDgsIGRwaT0zMDApCmBgYCAgICAKIVtdKC4uL091dHB1dC9DT1YvQ09Wc2Nhbl8zcG9wL1BXU19QSC52cy4zUG9wc19vdXRsaWVyX292ZXJsYXBfY292MjMucG5nKQoKCmBgYHtyIGV2YWw9RkFMU0UsIG1lc3NhZ2U9RkFMU0UsIHdhcm5pbmc9RkFMU0V9Cgpwd3MxPC1yZWFkLmNzdigiLi4vT3V0cHV0L0NPVi9DT1ZzY2FuXzNwb3AvM3BvcHNfdG9wMXBlcmNlbnRfb3V0bGllcl9yZWdpb25zLmNvdjEzX25ldy5jc3YiKQpwd3MxPC1wd3MxW3B3czEkcG9wPT0iUFdTIiAmIHB3czEkd2luZG93PT0iMTAwayIsXQpwd3MyPC1yZWFkLmNzdigiLi4vT3V0cHV0L0NPVi8zcG9wc190b3AxcGVyY2VudF9vdXRsaWVyX3JlZ2lvbnMuY292MTMuY3N2IikKcHdzMjwtcHdzMltwd3MyJHBvcD09IlBXUyImcHdzMiR3aW5kb3c9PSIxMDBrIixdCgpwd3MxJHZjZjwtIlBIIgpwd3MyJHZjZjwtIjNQb3BzIgpwd3M8LXJiaW5kKHB3czFbLGMoImNocm9tIiwic3RhcnQiLCJlbmQiLCJjb3YxMyIsInZjZiIpXSxwd3MyWyxjKCJjaHJvbSIsInN0YXJ0IiwiZW5kIiwiY292MTMiLCJ2Y2YiKV0pCnB3cyRjaHJvbTwtZmFjdG9yKHB3cyRjaHJvbSwgbGV2ZWxzPXBhc3RlMCgiY2hyIiwgMToyNikpCmdncGxvdChkYXRhPXB3cyxhZXMoeD1zdGFydC8xMDAwMDAwLCB5PWNvdjEzLCBjb2xvcj12Y2YsIGZpbGw9dmNmKSkrCiAgICBmYWNldF93cmFwKH5jaHJvbSkrCiAgICBnZW9tX3BvaW50KHBvc2l0aW9uPXBvc2l0aW9uX2RvZGdlKHdpZHRoID0gMC43KSwgYWxwaGE9MC41KSt4bGFiKCJQb3NpdGlvbiAoTWIpIikrCiAgICBnZ3RpdGxlKCJQV1MgQ09WMTMiKQpnZ3NhdmUoIi4uL091dHB1dC9DT1YvQ09Wc2Nhbl8zcG9wL1BXU19QSC52cy4zUG9wc19vdXRsaWVyX292ZXJsYXBfY292MTMucG5nIiwgd2lkdGggPSAxMCwgaGVpZ2h0ID0gOCwgZHBpPTMwMCkKYGBgICAgIAoKIVtdKC4uL091dHB1dC9DT1YvQ09Wc2Nhbl8zcG9wL1BXU19QSC52cy4zUG9wc19vdXRsaWVyX292ZXJsYXBfY292MTMucG5nKQoKCgoKCgojIEludGVycG9wdWxhdGlvbiBjb21wYXJpc29uIHBlciB0aW1lIHBlcmlvZCAKYGBge3IgZXZhbD1GQUxTRSwgbWVzc2FnZT1GQUxTRSwgd2FybmluZz1GQUxTRX0KIyMjIEludGVycG9wdWxhdGlvbiBjb21wYXJpc29ucwojZGVjb2RlIHRoZSBzYW1wbGVzIHRvIGNyZWF0ZSB0aGUgcmlnaHQgbWF0cml4CmN2PC1yZWFkLmNzdigifi9Qcm9qZWN0cy9QYWNoZXJyaW5nX1ZpbmNlbnQvTUQ3MDAwL0dXX0NvdnNfaW50ZXJQb3B1bGF0aW9uc18xMDBrXzNwb3BzLmNzdiIsIGhlYWRlciA9IEYpCmxhYnM8LXJlYWQuY3N2KCJ+L1Byb2plY3RzL1BhY2hlcnJpbmdfVmluY2VudC9NRDcwMDAvR1dfQ292c19pbnRlclBvcHVsYXRpb25zXzEwMGtfbGFiZWxzXzNwb3BzLmNzdiIgKQpsYWJzPC1sYWJzWywtMV0KY3ZtPC1kYXRhLmZyYW1lKGxhYmVsPWFzLnZlY3Rvcih0KGxhYnMpKSwgY292PWFzLnZlY3Rvcih0KGN2KSkpCgojcmVhcnJhbmdlIGJhc2VkIG9uIGNvbXBhcmlvbnM6IGNvdmFyaWFuY2UgYmV0d2VlbiBwb3B1bGF0aW9ucyB3aXRoaW4gdGhlIHNhbWUgcGVyaW9kCiNQb3BZciBTeW1ib2xzCiMgUEggMSAnUFdTJywgMTk5MQojIFBIIDIgJ1BXUycsIDE5OTYKIyBQSCAzICdQV1MnLCAyMDA2CiMgUEggNCAnUFdTJywgMjAxNwojIFBIIDUgJ1NTJywgIDE5OTEKIyBQSCA2ICdTUycsICAxOTk2CiMgUEggNyAnU1MnLCAgMjAwNgojIFBIIDggJ1NTJywgIDIwMTcKIyBQSCA5ICdUQicsICAxOTkxCiMgUEggMTAnVEInLCAgMTk5NgojIFBIIDExJ1RCJywgIDIwMDYKIyBQSCAxMidUQicsICAyMDE3CgpDb3ZzPC1kYXRhLmZyYW1lKHBvcHM9cmVwKGMoIlBXUy52cy5TUyIsICJQV1MudnMuVEIiLCAgIlNTLnZzLlRCIiksIHRpbWVzPTYpLAogICAgICAgICAgICAgICAgIHBlcmlvZD1jKHJlcCgiMTk5MS0xOTk2IiwgdGltZXM9MykscmVwKCIxOTk2LTIwMDYiLCB0aW1lcz0zKSwgcmVwKCIyMDA2LTIwMTciLCB0aW1lcz0zKSkpCgpDb3ZzJGNvdjwtYyhOQSwgY3ZtJGNvdltjdm0kbGFiZWw9PSJjb3YoUEg6IDItMSwgUEg6IDEwLTkpIl0sTkEsCiAgICAgICAgICAgIGN2bSRjb3ZbY3ZtJGxhYmVsPT0iY292KFBIOiAzLTIsIFBIOiA3LTYpIl0sY3ZtJGNvdltjdm0kbGFiZWw9PSJjb3YoUEg6IDMtMiwgUEg6IDExLTEwKSJdLCAKICAgICAgICAgICAgY3ZtJGNvdltjdm0kbGFiZWw9PSJjb3YoUEg6IDctNiwgUEg6IDExLTEwKSJdLAogICAgICAgICAgICBjdm0kY292W2N2bSRsYWJlbD09ImNvdihQSDogNC0zLCBQSDogOC03KSJdLGN2bSRjb3ZbY3ZtJGxhYmVsPT0iY292KFBIOiA0LTMsIFBIOiAxMi0xMSkiXSxjdm0kY292W2N2bSRsYWJlbD09ImNvdihQSDogOC03LCBQSDogMTItMTEpIl0pCgoKI0MuSS4KY2lzPC1yZWFkLmNzdigifi9Qcm9qZWN0cy9QYWNoZXJyaW5nX1ZpbmNlbnQvTUQ3MDAwL0dXX0NPVl9JbnRlcnBvcF9jb21wYXJpc29uX0NJcy5jc3YiKQpjaXM8LWNpc1ssLTFdCmNpbTwtZGF0YS5mcmFtZShsYWJlbD1hcy52ZWN0b3IodChsYWJzKSksIGNpX2w9YXMudmVjdG9yKHQoY2lzWzE6MTEsXSkpKQpjaW0kY2lfaDwtYXMudmVjdG9yKHQoY2lzWzEyOjIyLF0pKQoKQ292cyRjaV9sPC1hcy5udW1lcmljKGMoTkEsY2ltJGNpX2xbY2ltJGxhYmVsPT0iY292KFBIOiAyLTEsIFBIOiAxMC05KSJdLE5BLAogICAgICAgICAgICAgICAgICAgICAgY2ltJGNpX2xbY2ltJGxhYmVsPT0iY292KFBIOiAzLTIsIFBIOiA3LTYpIl0sY2ltJGNpX2xbY2ltJGxhYmVsPT0iY292KFBIOiAzLTIsIFBIOiAxMS0xMCkiXSwgY2ltJGNpX2xbY2ltJGxhYmVsPT0iY292KFBIOiA3LTYsIFBIOiAxMS0xMCkiXSwKICAgICAgICAgICAgICAgICAgICAgIGNpbSRjaV9sW2NpbSRsYWJlbD09ImNvdihQSDogNC0zLCBQSDogOC03KSJdLGNpbSRjaV9sW2NpbSRsYWJlbD09ImNvdihQSDogNC0zLCBQSDogMTItMTEpIl0sIGNpbSRjaV9sW2NpbSRsYWJlbD09ImNvdihQSDogOC03LCBQSDogMTItMTEpIl0pKQoKQ292cyRjaV9oPC1hcy5udW1lcmljKGMoTkEsIGNpbSRjaV9oW2NpbSRsYWJlbD09ImNvdihQSDogMi0xLCBQSDogMTAtOSkiXSxOQSwKICAgICAgICAgICAgICAgICAgICAgIGNpbSRjaV9oW2NpbSRsYWJlbD09ImNvdihQSDogMy0yLCBQSDogNy02KSJdLGNpbSRjaV9oW2NpbSRsYWJlbD09ImNvdihQSDogMy0yLCBQSDogMTEtMTApIl0sIGNpbSRjaV9oW2NpbSRsYWJlbD09ImNvdihQSDogNy02LCBQSDogMTEtMTApIl0sCiAgICAgICAgICAgICAgICAgICAgICBjaW0kY2lfaFtjaW0kbGFiZWw9PSJjb3YoUEg6IDQtMywgUEg6IDgtNykiXSxjaW0kY2lfaFtjaW0kbGFiZWw9PSJjb3YoUEg6IDQtMywgUEg6IDEyLTExKSJdLCBjaW0kY2lfaFtjaW0kbGFiZWw9PSJjb3YoUEg6IDgtNywgUEg6IDEyLTExKSJdKSkKCiNCYXJwbG90CmdncGxvdChDb3ZzLCBhZXMoeD1wZXJpb2QsIHk9Y292LCBmaWxsPXBvcHMpKSsKICAgIGdlb21fYmFyKHN0YXQ9ImlkZW50aXR5Iixwb3NpdGlvbj1wb3NpdGlvbl9kb2RnZSh3aWR0aCA9IDAuNyksIHdpZHRoPTAuOCkrCiAgICB5bGFiKCJDb3ZhcmlhbmNlIikreGxhYignJykrdGhlbWVfY2xhc3NpYygpKwogICAgZ2VvbV9obGluZSh5aW50ZXJjZXB0ID0gMCxjb2xvcj0iZ3JheTcwIiwgc2l6ZT0wLjMpKwogICAgc2NhbGVfZmlsbF9tYW51YWwodmFsdWVzPWNvbHNbYygyLDEsMyldKSsKICAgIHRoZW1lKGxlZ2VuZC50aXRsZSA9IGVsZW1lbnRfYmxhbmsoKSkrCiAgICBzY2FsZV95X2NvbnRpbnVvdXMobGFiZWxzID0gY29tbWEpKwogICAgeWxpbSgtMC4wMDEzLCAwLjAwMikrCiAgICBnZW9tX2Vycm9yYmFyKGFlcyh5bWluPWNpX2wsIHltYXg9Y2lfaCksIHdpZHRoPS4yLCBzaXplPS4yLCBwb3NpdGlvbj1wb3NpdGlvbl9kb2RnZSh3aWR0aCA9IDAuNykpCmdnc2F2ZSgiLi4vT3V0cHV0L0NPVi9JbnRlcnBvcF9jb3ZfY29tcGFyaXNvbl8zUG9wc19uZXcucG5nIix3aWR0aCA9IDQuNywgaGVpZ2h0ID0gMywgZHBpPTMwMCkKCiNQb2ludCBwbG90CmdncGxvdChDb3ZzLCBhZXMoeD1wZXJpb2QsIHk9Y292LCBjb2xvcj1wb3BzKSkrCiAgICBnZW9tX3BvaW50KHBvc2l0aW9uPXBvc2l0aW9uX2RvZGdlKHdpZHRoID0gMC43KSwgc2l6ZT00KSsKICAgIHlsYWIoIkNvdmFyaWFuY2UiKSt4bGFiKCcnKSt0aGVtZV9jbGFzc2ljKCkrCiAgICBnZW9tX2hsaW5lKHlpbnRlcmNlcHQgPSAwLGNvbG9yPSJncmF5NzAiLCBzaXplPTAuMykrCiAgICBzY2FsZV9jb2xvcl9tYW51YWwodmFsdWVzPWNvbHNbYygyLDEsMyldKSsKICAgIHRoZW1lKGxlZ2VuZC50aXRsZSA9IGVsZW1lbnRfYmxhbmsoKSkrCiAgICBzY2FsZV95X2NvbnRpbnVvdXMobGFiZWxzID0gY29tbWEpKwogICAgZ2VvbV9lcnJvcmJhcihhZXMoeW1pbj1jaV9sLCB5bWF4PWNpX2gpLCB3aWR0aD0uMiwgc2l6ZT0uMiwgcG9zaXRpb249cG9zaXRpb25fZG9kZ2Uod2lkdGggPSAwLjcpKSsKICAgIHlsaW0oLTAuMDAxMywgMC4wMDIpKwogICAgIGdlb21fdmxpbmUoeGludGVyY2VwdCA9IGMoMS41LDIuNSksIGNvbG9yPSJncmF5Iiwgc2l6ZT0wLjMpCmdnc2F2ZSgiLi4vT3V0cHV0L0NPVi9JbnRlcnBvcF9jb3ZfY29tcGFyaXNvbl8zUG9wc19uZXdfUG9pbnRQbG90LnBuZyIsd2lkdGggPSA0LjcsIGhlaWdodCA9IDMsIGRwaT0zMDApCgojbGluZSBwbG90CkNvdnMkdGltZTwtMQpDb3ZzJHRpbWVbQ292cyRwZXJpb2Q9PSIxOTk2LTIwMDYiXTwtMgpDb3ZzJHRpbWVbQ292cyRwZXJpb2Q9PSIyMDA2LTIwMTciXTwtMwpDb3ZzPC1Db3ZzW29yZGVyKENvdnMkdGltZSksXQpnZ3Bsb3QoQ292cywgYWVzKHg9dGltZSwgeT1jb3YsIGNvbG9yPXBvcHMsIGdyb3VwPXBvcHMpKSsKICAgIGdlb21fcG9pbnQocG9zaXRpb249cG9zaXRpb25fZG9kZ2Uod2lkdGggPSAwLjcpLCBzaXplPTQpKwogICAgZ2VvbV9wYXRoKHBvc2l0aW9uPXBvc2l0aW9uX2RvZGdlKHdpZHRoID0gMC43KSkrCiAgICB5bGFiKCJDb3ZhcmlhbmNlIikreGxhYignJykrdGhlbWVfY2xhc3NpYygpKwogICAgZ2VvbV9obGluZSh5aW50ZXJjZXB0ID0gMCxjb2xvcj0iZ3JheTcwIiwgc2l6ZT0wLjMpKwogICAgc2NhbGVfY29sb3JfbWFudWFsKHZhbHVlcz1jb2xzW2MoMiwxLDMpXSkrCiAgICB0aGVtZShsZWdlbmQudGl0bGUgPSBlbGVtZW50X2JsYW5rKCkpKwogICAgc2NhbGVfeV9jb250aW51b3VzKGxhYmVscyA9IGNvbW1hKSsKICAgIGdlb21fZXJyb3JiYXIoYWVzKHltaW49Y2lfbCwgeW1heD1jaV9oKSwgd2lkdGg9LjIsIHNpemU9LjIsIHBvc2l0aW9uPXBvc2l0aW9uX2RvZGdlKHdpZHRoID0gMC43KSkrCiAgICB5bGltKC0wLjAwMTMsIDAuMDAyKSsKICAgICBnZW9tX3ZsaW5lKHhpbnRlcmNlcHQgPSBjKDEuNSwyLjUpLCBjb2xvcj0iZ3JheSIsIHNpemU9MC4zKSsKICAgIHNjYWxlX3hfY29udGludW91cyhicmVha3M9YygxLDIsMyksIGxhYmVscyA9IGMoIjE5OTEtMTk5NiIsIjE5OTYtMjAwNiIsIjIwMDYtMjAxNyIpKQpnZ3NhdmUoIi4uL091dHB1dC9DT1YvSW50ZXJwb3BfY292X2NvbXBhcmlzb25fM1BvcHNfbmV3X0xpbmVQbG90LnBuZyIsd2lkdGggPSA0LjcsIGhlaWdodCA9IDMsIGRwaT0zMDApCgpgYGAKCiFbXSguLi9PdXRwdXQvQ09WL0ludGVycG9wX2Nvdl9jb21wYXJpc29uXzNQb3BzX25ld19Qb2ludFBsb3QucG5nKQoKIyMgTG9uZ2VyIHRpbWUgcGVyaW9kCmBgYHtyIGV2YWw9RkFMU0UsIG1lc3NhZ2U9RkFMU0UsIHdhcm5pbmc9RkFMU0V9CiMjIExvbmdlciB0aW1lLXBlcmlvZApDb3ZzMjwtZGF0YS5mcmFtZShwb3BzPXJlcChjKCJQV1MudnMuU1MiLCAiUFdTLnZzLlRCIiwgICJTUy52cy5UQiIpLCB0aW1lcz0zKSwKICAgICAgICAgICAgICAgICBwZXJpb2Q9YyhyZXAoIjE5OTEtMjAwNiIsIHRpbWVzPTMpLHJlcCgiMTk5MS0yMDE3IiwgdGltZXM9MykscmVwKCIxOTk2LTIwMTciLCB0aW1lcz0zKSkpCgpjdjE8LXJlYWQuY3N2KCJ+L1Byb2plY3RzL1BhY2hlcnJpbmdfVmluY2VudC9NRDcwMDAvR1dfQ292c19pbnRlclBvcHVsYXRpb25zXzEwMGtfMTk5MS0yMDA2LmNzdiIsIGhlYWRlciA9IEYpCmxhYnMxPC1yZWFkLmNzdigifi9Qcm9qZWN0cy9QYWNoZXJyaW5nX1ZpbmNlbnQvTUQ3MDAwL0dXX0NvdnNfaW50ZXJQb3B1bGF0aW9uc18xMDBrX2xhYmVsc18xOTkxLTIwMDYuY3N2IiApCmxhYnMxPC1sYWJzMVssLTFdCmN2bTE8LWRhdGEuZnJhbWUobGFiZWw9YXMudmVjdG9yKHQobGFiczEpKSwgY292PWFzLnZlY3Rvcih0KGN2MSkpKQoKY3YyPC1yZWFkLmNzdigifi9Qcm9qZWN0cy9QYWNoZXJyaW5nX1ZpbmNlbnQvTUQ3MDAwL0dXX0NvdnNfaW50ZXJQb3B1bGF0aW9uc18xMDBrXzE5OTEtMjAxNy5jc3YiLCBoZWFkZXIgPSBGKQpsYWJzMjwtcmVhZC5jc3YoIn4vUHJvamVjdHMvUGFjaGVycmluZ19WaW5jZW50L01ENzAwMC9HV19Db3ZzX2ludGVyUG9wdWxhdGlvbnNfMTAwa19sYWJlbHNfMTk5MS0yMDE3LmNzdiIgKQpsYWJzMjwtbGFiczJbLC0xXQpjdm0yPC1kYXRhLmZyYW1lKGxhYmVsPWFzLnZlY3Rvcih0KGxhYnMyKSksIGNvdj1hcy52ZWN0b3IodChjdjIpKSkKCmN2MzwtcmVhZC5jc3YoIn4vUHJvamVjdHMvUGFjaGVycmluZ19WaW5jZW50L01ENzAwMC9HV19Db3ZzX2ludGVyUG9wdWxhdGlvbnNfMTAwa18xOTk2LTIwMTcuY3N2IiwgaGVhZGVyID0gRikKbGFiczM8LXJlYWQuY3N2KCJ+L1Byb2plY3RzL1BhY2hlcnJpbmdfVmluY2VudC9NRDcwMDAvR1dfQ292c19pbnRlclBvcHVsYXRpb25zXzEwMGtfbGFiZWxzXzE5OTYtMjAxNy5jc3YiICkKbGFiczM8LWxhYnMzWywtMV0KY3ZtMzwtZGF0YS5mcmFtZShsYWJlbD1hcy52ZWN0b3IodChsYWJzMykpLCBjb3Y9YXMudmVjdG9yKHQoY3YzKSkpCgpDb3ZzMiRjb3Y8LWMoTkEsIGN2bTEkY292W2N2bTEkbGFiZWw9PSJjb3YoUEg6IDItMSwgUEg6IDQtMykiXSwgTkEsCiAgICAgICAgICAgICBOQSwgY3ZtMiRjb3ZbY3ZtMiRsYWJlbD09ImNvdihQSDogMi0xLCBQSDogNC0zKSJdLCBOQSwKICAgICAgICAgICAgIGN2bTMkY292W2N2bTMkbGFiZWw9PSJjb3YoUEg6IDItMSwgUEg6IDQtMykiXSwgY3ZtMyRjb3ZbY3ZtMyRsYWJlbD09ImNvdihQSDogMi0xLCBQSDogNi01KSJdLCBjdm0zJGNvdltjdm0zJGxhYmVsPT0iY292KFBIOiA0LTMsIFBIOiA2LTUpIl0pCgojQy5JLgpjaXMxPC1yZWFkLmNzdigifi9Qcm9qZWN0cy9QYWNoZXJyaW5nX1ZpbmNlbnQvTUQ3MDAwL0dXX0NPVl9JbnRlcnBvcF9jb21wYXJpc29uX0NJc18xOTkxLTIwMDYuY3N2IikKY2lzMTwtY2lzMVssLTFdCmNpczI8LXJlYWQuY3N2KCJ+L1Byb2plY3RzL1BhY2hlcnJpbmdfVmluY2VudC9NRDcwMDAvR1dfQ09WX0ludGVycG9wX2NvbXBhcmlzb25fQ0lzXzE5OTEtMjAxNy5jc3YiKQpjaXMyPC1jaXMyWywtMV0KY2lzMzwtcmVhZC5jc3YoIn4vUHJvamVjdHMvUGFjaGVycmluZ19WaW5jZW50L01ENzAwMC9HV19DT1ZfSW50ZXJwb3BfY29tcGFyaXNvbl9DSXNfMTk5Ni0yMDE3LmNzdiIpCmNpczM8LWNpczNbLC0xXQoKI2NpbTwtZGF0YS5mcmFtZShsYWJlbD1hcy52ZWN0b3IodChsYWJzKSksIGNpX2w9YXMudmVjdG9yKHQoY2lzMVsxOjQsXSkpKQojY2ltJGNpX2g8LWFzLnZlY3Rvcih0KGNpc1sxMjoyMixdKSkKCkNvdnMyJGNpX2w8LWFzLm51bWVyaWMoYyhOQSxjaXMxWzEsM10sTkEsCiAgICAgICAgICAgICAgICAgICAgICAgIE5BLGNpczJbMSwzXSxOQSwKICAgICAgICAgICAgICAgICAgICAgIGNpczNbMSwzXSxjaXMzWzEsNV0sY2lzM1szLDVdKSkKCkNvdnMyJGNpX2g8LWFzLm51bWVyaWMoYyhOQSxjaXMxWzQsM10sTkEsCiAgICAgICAgICAgICAgICAgICAgICAgIE5BLGNpczJbNCwzXSxOQSwKICAgICAgICAgICAgICAgICAgICAgIGNpczNbNiwzXSxjaXMzWzYsNV0sY2lzM1s4LDVdKSkKCgpnZ3Bsb3QoQ292czIsIGFlcyh4PXBlcmlvZCwgeT1jb3YsIGZpbGw9cG9wcykpKwogICAgZ2VvbV9iYXIoc3RhdD0iaWRlbnRpdHkiLHBvc2l0aW9uPXBvc2l0aW9uX2RvZGdlKHdpZHRoID0gMC43KSwgd2lkdGg9MC44KSsKICAgIHlsYWIoIkNvdmFyaWFuY2UiKSt4bGFiKCcnKSt0aGVtZV9jbGFzc2ljKCkrCiAgICBnZW9tX2hsaW5lKHlpbnRlcmNlcHQgPSAwLGNvbG9yPSJncmF5NzAiLCBzaXplPTAuMykrCiAgICBzY2FsZV9maWxsX21hbnVhbCh2YWx1ZXM9Y29sc1tjKDIsMSwzKV0pKwogICAgdGhlbWUobGVnZW5kLnRpdGxlID0gZWxlbWVudF9ibGFuaygpKStzY2FsZV95X2NvbnRpbnVvdXMobGFiZWxzID0gY29tbWEpKwogICAgeWxpbSgtMC4wMDEzLCAwLjAwMikrCiAgICBnZW9tX2Vycm9yYmFyKGFlcyh5bWluPWNpX2wsIHltYXg9Y2lfaCksIHdpZHRoPS4yLCBzaXplPS4yLCBwb3NpdGlvbj1wb3NpdGlvbl9kb2RnZSh3aWR0aCA9IDAuNykpCmdnc2F2ZSgiLi4vT3V0cHV0L0NPVi9JbnRlcnBvcF9jb3ZfY29tcGFyaXNvbl8zUG9wc19Mb25vZ2VyUGVyaW9kLnBuZyIsd2lkdGggPSA0LjcsIGhlaWdodCA9IDMsIGRwaT0zMDApCgpnZ3Bsb3QoQ292czIsIGFlcyh4PXBlcmlvZCwgeT1jb3YsIGNvbG9yPXBvcHMpKSsKICAgIGdlb21fcG9pbnQocG9zaXRpb249cG9zaXRpb25fZG9kZ2Uod2lkdGggPSAwLjcpLCBzaXplPTQpKwogICAgeWxhYigiQ292YXJpYW5jZSIpK3hsYWIoJycpK3RoZW1lX2NsYXNzaWMoKSsKICAgIGdlb21faGxpbmUoeWludGVyY2VwdCA9IDAsY29sb3I9ImdyYXk3MCIsIHNpemU9MC4zKSsKICAgIHNjYWxlX2NvbG9yX21hbnVhbCh2YWx1ZXM9Y29sc1tjKDIsMSwzKV0pKwogICAgdGhlbWUobGVnZW5kLnRpdGxlID0gZWxlbWVudF9ibGFuaygpKSsKICAgIHNjYWxlX3lfY29udGludW91cyhsYWJlbHMgPSBjb21tYSkrCiAgICAgeWxpbSgtMC4wMDEzLCAwLjAwMikrCiAgICBnZW9tX2Vycm9yYmFyKGFlcyh5bWluPWNpX2wsIHltYXg9Y2lfaCksIHdpZHRoPS4yLCBzaXplPS4yLCBwb3NpdGlvbj1wb3NpdGlvbl9kb2RnZSh3aWR0aCA9IDAuNykpKwogICAgZ2VvbV92bGluZSh4aW50ZXJjZXB0ID0gYygxLjUsMi41KSwgY29sb3I9ImdyYXkiLCBzaXplPTAuMykKZ2dzYXZlKCIuLi9PdXRwdXQvQ09WL0ludGVycG9wX2Nvdl9jb21wYXJpc29uXzNQb3BzTG9ub2dlclBlcmlvZF9Qb2ludFBsb3QucG5nIix3aWR0aCA9IDQuNywgaGVpZ2h0ID0gMywgZHBpPTMwMCkKCgpgYGAKIVtdKC4uL091dHB1dC9DT1YvSW50ZXJwb3BfY292X2NvbXBhcmlzb25fM1BvcHNMb25vZ2VyUGVyaW9kX1BvaW50UGxvdC5wbmcpCgojIEZvY3VzZWQgZnJlcSBhbmFseXNpcwojIyByZWwgZ2VuZSAoY2hyMTM6IDIzMDcwMDAwIC0gMjMwODAwMDApCgpgYGB7ciBldmFsPUZBTFNFLCBtZXNzYWdlPUZBTFNFLCB3YXJuaW5nPUZBTFNFfQpwb3BzPC1jKCJQV1M5MSIsIlBXUzk2IiwiUFdTMDciLCJQV1MxNyIpCnlyPC1jKDE5OTEsMTk5NiwyMDA3LDIwMTcpCm1hZjwtZGF0YS5mcmFtZSgpCmZvciAoaSBpbiAxOjQpewogICAgYWY8LXJlYWQudGFibGUocGFzdGUwKCIuLi9EYXRhL25ld192Y2YvQUYvIixwb3BzW2ldLCIubWFmcyIpLHNlcD0iXHQiLCBoZWFkZXIgPSBUKQogICAgYWY8LWFmW2FmJGNocm9tbz09ImNocjEzIiZhZiRwb3NpdGlvbj49MjMwNzAwMDAmYWYkcG9zaXRpb248PTIzMDgwMDAwLF0KICAgIGFmJHllYXI8LXlyW2ldCiAgICBtYWY8LXJiaW5kKG1hZixhZikKfQojd3JpdGUuY3N2KG1hZiwiLi4vT3V0cHV0L0NPVi9DT1ZzY2FuXzNwb3AvQUZfbWFmX2NocjEzXzIzTWIuY3N2IikKCmdncGxvdChtYWYsIGFlcyh4PXllYXIsIHk9a25vd25FTSwgY29sb3I9ZmFjdG9yKHBvc2l0aW9uKSkpKwogICAgZ2VvbV9wb2ludCgpKwogICAgZ2VvbV9wYXRoKCkrZ2d0aXRsZSgiTUFGIChBTkdTRCkiKSsKICAgIHlsYWIoIm1hZiIpKwogICAgdGhlbWUobGVnZW5kLnRpdGxlPWVsZW1lbnRfYmxhbmsoKSkKZ2dzYXZlKCIuLi9PdXRwdXQvQ09WL0NPVnNjYW4vQUZfY2gxM18yMy4wNy0yMy4wOF9hbmdzZC5wbmciLCB3aWR0aCA9IDUsIGhlaWdodD0zLCBkcGk9MzAwKQoKCkFGPC1kYXRhLmZyYW1lKCkKZm9yIChpIGluIDE6NCl7CiAgICBhZjwtcmVhZC50YWJsZShwYXN0ZTAoIi4uL0RhdGEvbmV3X3ZjZi9BRi8iLHBvcHNbaV0sIl9tYWYwNV9hZi5mcnEiKSxoZWFkZXIgPSBGQUxTRSwgc2tpcD0xLCBjb2wubmFtZXMgPSBjKCJjaHIiLCJwb3MiLCJuX2FsbGVsZSIsIm5fc2FtcGxlIiwiTWFqb3JBRiIsIk1BRiIpKQogICAgYWY8LWFmW2FmJGNocj09ImNocjEzIiZhZiRwb3M+PTIzMDcwMDAwJmFmJHBvczw9MjMwODAwMDAsXQogICAgYWYkeWVhcjwteXJbaV0KICAgIGFmJG1hZjwtc3Vic3RyKGFmJE1BRiwgMywxMCkKICAgIGFmJG1hZjwtYXMubnVtZXJpYyhhZiRtYWYpCiAgICBhZjwtYWZbLGMoImNociIsInBvcyIsIm1hZiIsInllYXIiKV0KICAgIEFGPC1yYmluZChBRixhZikKfQpnZ3Bsb3QoQUYsIGFlcyh4PXllYXIsIHk9bWFmLCBjb2xvcj1mYWN0b3IocG9zKSkpKwogICAgZ2VvbV9wb2ludCgpKwogICAgZ2VvbV9wYXRoKCkrZ2d0aXRsZSgiTUFGICh2Y2Z0b29scykiKSsKICAgIHRoZW1lKGxlZ2VuZC50aXRsZT1lbGVtZW50X2JsYW5rKCkpCmdnc2F2ZSgiLi4vT3V0cHV0L0NPVi9DT1ZzY2FuL0FGX2NoMTNfMjMuMDctMjMuMDgucG5nIiwgd2lkdGggPSA1LCBoZWlnaHQ9MywgZHBpPTMwMCkKCgojIyNUQgpwb3BzPC1jKCJUQjkxIiwiVEI5NiIsIlRCMDYiLCJUQjE3IiwiUFdTOTEiLCJQV1M5NiIsIlBXUzA3IiwiUFdTMTciLCJTUzk2IiwiU1MwNiIsIlNTMTciKQp5cjwtYygxOTkxLDE5OTYsMjAwNiwyMDE3LDE5OTEsMTk5NiwyMDA3LDIwMTcsMTk5NiwyMDA2LDIwMTcpCm1hZjwtZGF0YS5mcmFtZSgpCmZvciAoaSBpbiAxOmxlbmd0aChwb3BzKSl7CiAgICBhZjwtcmVhZC50YWJsZShwYXN0ZTAoIi4uL0RhdGEvbmV3X3ZjZi9BRi8iLHBvcHNbaV0sIi5tYWZzIiksc2VwPSJcdCIsIGhlYWRlciA9IFQpCiAgICBhZjwtYWZbYWYkY2hyb21vPT0iY2hyMTMiJmFmJHBvc2l0aW9uPj0yMzA3MDAwMCZhZiRwb3NpdGlvbjw9MjMwODAwMDAsXQogICAgYWYkeWVhcjwteXJbaV0KICAgIGFmJHBvcDwtc3ViKCJcXGRcXGQiLCIiLCBwb3BzW2ldKQogICAgbWFmPC1yYmluZChtYWYsYWYpCn0KI3dyaXRlLmNzdihtYWYsIi4uL091dHB1dC9DT1YvQ09Wc2Nhbl8zcG9wL0FGX21hZl9jaHIxM18yM01iLmNzdiIpCmdncGxvdChtYWYsIGFlcyh4PXllYXIsIHk9a25vd25FTSwgY29sb3I9cG9wKSkrCiAgICBmYWNldF93cmFwKH5mYWN0b3IocG9zaXRpb24pKSsKICAgIGdlb21fcG9pbnQoKSsKICAgIGdlb21fcGF0aCgpK2dndGl0bGUoIkNocjEzIChyZWwgZ2VuZSkiKSsKICAgIHlsYWIoIm1hZiIpKwogICAgdGhlbWUobGVnZW5kLnRpdGxlPWVsZW1lbnRfYmxhbmsoKSkKZ2dzYXZlKCIuLi9PdXRwdXQvQ09WL0NPVnNjYW4vQUZfY2gxM18yMy4wNy0yMy4wOF9hbmdzZC5wbmciLCB3aWR0aCA9IDUsIGhlaWdodD0zLCBkcGk9MzAwKQoKCgoKYGBgCiFbXSguLi9PdXRwdXQvQ09WL0NPVnNjYW4vQUZfY2gxM18yMy4wNy0yMy4wOF9hbmdzZC5wbmcpCgoKCgoKYGBge3IgZXZhbD1GQUxTRSwgbWVzc2FnZT1GQUxTRSwgd2FybmluZz1GQUxTRX0KCgoKCmBgYAoKYGBge3IgZXZhbD1GQUxTRSwgbWVzc2FnZT1GQUxTRSwgd2FybmluZz1GQUxTRX0KYGBgCgo=